&star_job

      ! load initial model from inlist_150M
      load_saved_model = .true.
      saved_model_name = 'ready_to_bounce.mod'
      !saved_model_name = 'post_bounce1.mod'

      ! output final profile
      save_model_when_terminate = .true.
      save_model_filename = 'post_bounce1.mod'
      
      ! output termination file
      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'final_profile.data'

      ! copy from ccsn but is set to zero
      set_initial_L_center = .true.
      new_L_center = 0d0

      ! copy from ccsn but is set to zero
      set_initial_v_center = .true.
      new_v_center = 0

      ! continue from previous run
      set_initial_age = .false.

      ! but reset the setno for easier reference
      set_initial_model_number = .true.
      initial_model_number = 0

      ! try to remove the low density matter if it exists at the beginning
      remove_surface_by_density = 1e-3

      ! limit on timestep
      !set_initial_dt = .true.
      !years_for_initial_dt = 1d-14

      !pgstar_flag = .true.

/ ! end of star_job namelist

&controls

      ! set up a Courant factor which does not over-limit the simulation
      dt_div_dt_explicit_limit = 5.0D3  ! 0.5 ! use 0.5 for edep phase.

      ! check for retries and backups as part of test_suite
      ! you can/should delete this for use outside of test_suite
      !   max_number_backups = 100
      !   max_number_retries = 1200
      max_model_number = 20000
	 
      ! some old setting
      initial_mass = 4.0D+01
      initial_z = 0.02d0
      Zbase = 0.02d0

      ! no termination when Fe-core collapse
      fe_core_infall_limit = 1d99 ! cm/sec

      ! original from ccsn
      velocity_logT_lower_bound = 0

      ! original from ccsn
      ! max_model_number = 5
      ! max_age = 3.01e5 ! 3.1688e-08 ! >~ 1sec (in year)

      ! equation controls
      !explicit_hydro_max_nsteps = 1
      !explicit_hydro_CFL = 0.5d0
      !use_lagrangian_riemann_solver =.true.

      use_energy_conservation_form = .true.
     
      use_ODE_var_eqn_pairing = .true.
      use_dvdt_form_of_momentum_eqn = .true.
      use_dPrad_dm_form_of_T_gradient_eqn = .true.
      use_dedt_form_of_energy_eqn = .true.
      use_momentum_outer_BC = .true.      
      use_ODE_form_of_density_eqn = .true.     ! no diffusion of rho unless this is true  

      ! No Wind only allow dynamical wind
      Dutch_scaling_factor = 0.0d0

      ! mixing controls
      !mixing_length_alpha = 1.5
      MLT_option = 'none'
      !T_mix_limit = 0
      !min_T_for_acceleration_limited_conv_velocity = 0
      !okay_to_reduce_gradT_excess = .false.

      ! atmosphere boundary conditions
      Pextra_factor = 0

      ! artificial viscosity controls
      use_artificial_viscosity = .true.
      shock_spread_linear = 0
      shock_spread_quadratic = 2d-2 

      ! mesh controls
      okay_to_remesh = .true.
      mesh_delta_coeff = 0.8
      min_dq = 1d-4
      max_center_cell_dq = 1d-4
      max_surface_cell_dq = 1d-4

      !mach1_minus_dr_factor = 1d-2
      !mach1_plus_dr_factor = 0d0

      ! output controls
      max_num_profile_models = 1000
      terminal_show_age_in_years = .false.

      photostep = 50
      profile_interval = 100 ! 50
      history_interval = 2  ! 1
      terminal_interval = 10
      write_header_frequency = 10

      ! output on screen for some insights
      num_trace_history_values = 6
      trace_history_value_name(1) = 'dt_div_dt_explicit'
      trace_history_value_name(2) = 'total_energy'
      trace_history_value_name(3) = 'log_max_abs_lgE_residual'
      trace_history_value_name(4) = 'log_max_abs_v_residual'
      trace_history_value_name(5) = 'total_mass ni56'
      trace_history_value_name(6) = 'shock_mass'

      ! timestep controls
      delta_lgL_He_limit = 0.5
      delta_lgT_cntr_limit = 0.2
      delta_lgRho_cntr_limit = 0.5

      dX_nuc_drop_limit_at_high_T = 2d-2 !1e-3

      ! relax iteration
      varcontrol_target = 1d-3
      

      ! solver controls
      newton_iterations_limit = 7
      iter_for_resid_tol2 = 4
      tol_residual_norm1 = 1d-8
      tol_max_residual1 = 1d-7

      tiny_corr_coeff_limit = 999999
      newton_itermin_until_reduce_min_corr_coeff = 999999


      ! debugging controls
      !report_hydro_solver_progress = .true. ! set true to see info about newton iterations
      !report_ierr = .true. ! if true, produce terminal output when have some internal error
      !report_all_dt_limits = .true.
      !report_why_dt_limits = .true.
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 26545
      !trace_evolve = .true.

      ! hydro debugging
      !hydro_check_everything = .true.
      !hydro_inspectB_flag = .true.
      !hydro_numerical_jacobian = .true.
      !hydro_save_numjac_plot_data = .true.
      !small_mtx_decsol = 'lapack'
      !large_mtx_decsol = 'lapack'
      !hydro_dump_call_number = 6836


      !dX_nuc_drop_limit = 1d-3


/ ! end of controls namelist




&pgstar

!pause = .true.

         Grid7_title = '15 M\d\(2281)\u CCSN - without RT'

         Grid7_win_flag = .true.
         Grid7_win_width = 11
         Grid7_win_aspect_ratio = 0.7 ! aspect_ratio = height/width

         Grid7_num_cols = 2 ! divide plotting region into this many equal width cols
         Grid7_num_rows = 5 ! divide plotting region into this many equal height rows
         Grid7_num_plots = 3 ! <= 10


         Grid7_plot_name(1) = 'Profile_Panels1'
         Grid7_plot_row(1) = 1 ! number from 1 at top
         Grid7_plot_rowspan(1) = 4 ! plot spans this number of rows
         Grid7_plot_col(1) =  1 ! number from 1 at left
         Grid7_plot_colspan(1) = 1 ! plot spans this number of columns
         Grid7_plot_pad_left(1) = 0.0 ! fraction of full window width for padding on left
         Grid7_plot_pad_right(1) = 0.10 ! fraction of full window width for padding on right
         Grid7_plot_pad_top(1) = -0.02 ! fraction of full window height for padding at top
         Grid7_plot_pad_bot(1) = -0.02 ! fraction of full window height for padding at bottom
         Grid7_txt_scale_factor(1) = 0.6 ! multiply txt_scale for subplot by this


         Grid7_plot_name(2) = 'Profile_Panels2'
         Grid7_plot_row(2) = 1 ! number from 1 at top
         Grid7_plot_rowspan(2) = 4 ! plot spans this number of rows
         Grid7_plot_col(2) =  2 ! number from 1 at left
         Grid7_plot_colspan(2) = 1 ! plot spans this number of columns
         Grid7_plot_pad_left(2) = 0.06 ! fraction of full window width for padding on left
         Grid7_plot_pad_right(2) = 0.04 ! fraction of full window width for padding on right
         Grid7_plot_pad_top(2) = -0.02 ! fraction of full window height for padding at top
         Grid7_plot_pad_bot(2) = -0.02 ! fraction of full window height for padding at bottom
         Grid7_txt_scale_factor(2) = 0.6 ! multiply txt_scale for subplot by this



         Grid7_plot_name(3) = 'Text_Summary1'
         Grid7_plot_row(3) = 5 ! number from 1 at top
         Grid7_plot_rowspan(3) = 1 ! plot spans this number of rows
         Grid7_plot_col(3) =  1 ! number from 1 at left
         Grid7_plot_colspan(3) = 2 ! plot spans this number of columns
         Grid7_plot_pad_left(3) = -0.08 ! fraction of full window width for padding on left
         Grid7_plot_pad_right(3) = -0.03 ! fraction of full window width for padding on right
         Grid7_plot_pad_top(3) = 0.09 ! fraction of full window height for padding at top
         Grid7_plot_pad_bot(3) = -0.05 ! fraction of full window height for padding at bottom
         Grid7_txt_scale_factor(3) = 0.21 ! multiply txt_scale for subplot by this



         ! file output
         Grid7_file_flag = .true.
         Grid7_file_dir = 'png'
         Grid7_file_prefix = 'ccsn_'
         Grid7_file_cnt = 25 ! output when mod(model_number,Grid7_file_cnt)==0





         Profile_Panels1_title = '' ! 'Shock Profiles'
         Profile_Panels2_title = '' ! 'Abundances + Power'


         Profile_Panels_show_Mach_1_location = .true.
         Profile_Panels_xwidth_left_of_shock = -1
	 Profile_Panels_xwidth_right_of_shock = 0.1
	 !Profile_Panels_xwidth_left_div_shock_value = 1
	 !Profile_Panels_xwidth_right_div_shock_value = 0.3


         Profile_Panels1_xaxis_name =  'mass' !  'logR_cm' !
 	 Profile_Panels1_xmin = -101d0 ! only used if /= -101d0
	 Profile_Panels1_xmax = 1.85 ! -101d0 ! 1.54 ! 2 ! 12 ! only used if /= -101d0
	 !Profile_Panels1_xmin = 3.99853 ! only used if /= -101d0
	 !Profile_Panels1_xmax = 11 ! 3.99855 ! only used if /= -101d0

         Profile_Panels1_num_panels = 4

         Profile_Panels1_yaxis_name(1) = 'logT'
         Profile_Panels1_ymin(1) = -101d0 ! -1.2 ! only used if /= -101d0
         Profile_Panels1_ymax(1) = -101d0 ! 3 ! only used if /= -101d0
         Profile_Panels1_ymargin(1) = 0.1

         Profile_Panels1_other_yaxis_name(1) = 'vel_km_per_s' 
         Profile_Panels1_other_ymin(1) = -0.1 ! -1e4 ! only used if /= -101d0
         Profile_Panels1_other_ymax(1) = -101d0 ! 2e4 ! only used if /= -101d0
         Profile_Panels1_other_ymargin(1) = 0.1

         Profile_Panels1_yaxis_name(2) = 'logRho'
         Profile_Panels1_ymin(2) = -101d0 ! -1.2 ! only used if /= -101d0
         Profile_Panels1_ymax(2) = -101d0 ! 3 ! only used if /= -101d0
         Profile_Panels1_ymargin(2) = 0.1

         Profile_Panels1_other_yaxis_name(2) = 'logP'
         Profile_Panels1_other_ymin(2) = -101d0 ! -1e4 ! only used if /= -101d0
         Profile_Panels1_other_ymax(2) = -101d0 ! 2e4 ! only used if /= -101d0
         Profile_Panels1_other_ymargin(2) = 0.1

         Profile_Panels1_yaxis_name(3) = 'log_alpha_RTI'
         Profile_Panels1_ymin(3) = -8 ! -1.2 ! only used if /= -101d0
         Profile_Panels1_ymax(3) = -101d0 ! 3 ! only used if /= -101d0
         Profile_Panels1_ymargin(3) = 0.001

         Profile_Panels1_other_yaxis_name(3) =  'log_source_plus_alpha_RTI'  ! 'logR_cm'
         Profile_Panels1_other_ymin(3) = -8 ! -1e4 ! only used if /= -101d0
         Profile_Panels1_other_ymax(3) = -101d0 ! 2e4 ! only used if /= -101d0
         Profile_Panels1_other_ymargin(3) = 0.1

         Profile_Panels1_yaxis_name(4) = 'logR_cm'
         Profile_Panels1_ymin(4) = -101d0 ! -1.2 ! only used if /= -101d0
         Profile_Panels1_ymax(4) = -101d0 ! 3 ! only used if /= -101d0
         Profile_Panels1_ymargin(4) = 0.1

         Profile_Panels1_other_yaxis_name(4) =  'entropy'  ! 'logR_cm'
         Profile_Panels1_other_ymin(4) = -101d0 ! -1e4 ! only used if /= -101d0
         Profile_Panels1_other_ymax(4) = -101d0 ! 2e4 ! only used if /= -101d0
         Profile_Panels1_other_ymargin(4) = 0.1


         Profile_Panels2_xaxis_name =  'mass' !  'logR_cm' !
         Profile_Panels2_xaxis_reversed = .false.
	 Profile_Panels2_xmin = -101d0 ! only used if /= -101d0
	 Profile_Panels2_xmax = 1.85 ! -101d0 ! 1.54 ! 2 ! 12 ! only used if /= -101d0

         Profile_Panels2_num_panels = 4

         Profile_Panels2_yaxis_name(1) = 'Abundance'
         Profile_Panels2_ymin(1) = -3 ! only used if /= -101d0
         Profile_Panels2_ymax(1) = 0.1 ! only used if /= -101d0
         Profile_Panels2_ymargin(1) = 0.1

         Profile_Panels2_yaxis_name(2) = 'Mixing'
            Profile_Panels2_ymin(2) = 0 ! only used if /= -101d0
            Profile_Panels2_ymax(2) = 20 ! only used if /= -101d0
            Profile_Panels2_ymargin(2) = 0.1

         Profile_Panels2_yaxis_name(3) = 'v_div_r'
            Profile_Panels2_ymin(3) = -101d0 ! -1.2 ! only used if /= -101d0
            Profile_Panels2_ymax(3) = -101d0 ! 3 ! only used if /= -101d0
            Profile_Panels2_ymargin(3) = 0.1

         Profile_Panels2_other_yaxis_name(3) = 'cell_IE_div_IE_plus_KE'  ! 'logR_cm'
            Profile_Panels2_other_ymin(3) = 0d0 ! -1e4 ! only used if /= -101d0
            Profile_Panels2_other_ymax(3) = -101d0 ! 2e4 ! only used if /= -101d0
            Profile_Panels2_other_ymargin(3) = 0.1

         Profile_Panels2_yaxis_name(4) = 'zone'
            Profile_Panels2_ymin(4) = -101d0 ! -1.2 ! only used if /= -101d0
            Profile_Panels2_ymax(4) = -101d0 ! 3 ! only used if /= -101d0
            Profile_Panels2_ymargin(4) = 0.1

         Profile_Panels2_other_yaxis_name(4) = 'logdq'
            Profile_Panels2_other_ymin(4) = -101d0 ! -1e4 ! only used if /= -101d0
            Profile_Panels2_other_ymax(4) = -101d0 ! 2e4 ! only used if /= -101d0
            Profile_Panels2_other_ymargin(4) = 0.1

         Power_legend_max_cnt = 10
         Power_legend_txt_scale_factor = 1.3 ! relative to other text
         Power_ymin = 6 ! only used if /= -101d0

         Text_Summary1_num_cols = 4
         Text_Summary1_num_rows = 5

         Text_Summary1_name(1,1) = 'model_number'
         Text_Summary1_name(2,1) = 'star_age_sec'
         Text_Summary1_name(3,1) = 'time_step_sec'
         Text_Summary1_name(4,1) = 'shock_mass'
         Text_Summary1_name(5,1) = 'shock_radius'
         Text_Summary1_name(6,1) = ''
         Text_Summary1_name(7,1) = ''
         Text_Summary1_name(8,1) = ''

         Text_Summary1_name(1,2) = 'log_L'
         Text_Summary1_name(2,2) = 'log_Teff'
         Text_Summary1_name(3,2) = 'log_R'
         Text_Summary1_name(4,2) = 'log_cntr_Rho'
         Text_Summary1_name(5,2) = 'log_cntr_T'
         Text_Summary1_name(6,2) = ''
         Text_Summary1_name(7,2) = ''
         Text_Summary1_name(8,2) = ''

         Text_Summary1_name(1,3) = 'log_Lnuc'
         Text_Summary1_name(2,3) = 'log_Lneu'
         Text_Summary1_name(3,3) = 'total_energy'
         Text_Summary1_name(4,3) = 'total_IE_plus_KE'
         Text_Summary1_name(5,3) = 'total_mass ni56'
         Text_Summary1_name(6,3) = ''
         Text_Summary1_name(7,3) = ''
         Text_Summary1_name(8,3) = ''

         Text_Summary1_name(1,4) = 'star_mass'
         Text_Summary1_name(2,4) = 'species'
         Text_Summary1_name(3,4) = 'num_zones'
         Text_Summary1_name(4,4) = 'num_retries'
         Text_Summary1_name(5,4) = 'num_backups'
         Text_Summary1_name(6,4) = ''
         Text_Summary1_name(7,4) = ''
         Text_Summary1_name(8,4) = ''

      Grid2_win_flag = .false.


/ ! end of pgstar namelist
