
&star_job

         mesa_data_dir = '../../../data'

         show_log_description_at_start = .false. 
         show_net_reactions_info = .true.

         !create_pre_main_sequence_model = .true.

         change_net = .true.      
         !new_net_name = 'net_A32.jina'
         new_net_name = 'approx21.net'

         load_saved_model = .true.
         saved_model_name = '15M_si_burn.mod'

         save_model_when_terminate = .true.
         save_model_filename = 'si_burn_all_done.mod'

         report_retries = .false.
         report_backups = .false.

         log_columns_file = 'log_columns.list' 
         profile_columns_file = 'profile_columns.list' 

         set_rate_c12ag = 'Buchmann'

         num_special_rate_factors = 1
         reaction_for_special_factor(:) = 'r_c12_ag_o16'
         special_rate_factor(1) = 1.1
         
         kappa_file_prefix = 'gs98'


         !pgstar_flag = .true.
         

/ ! end of star_job namelist


&controls

         D_norm_kappa = 1
         D_norm_alpha = 2
         D_norm_atol = 1d-3
         D_norm_rtol = 1d-3
      
         use_lnE_for_eps_grav = .true.
      
         initial_mass = 15
         initial_z = 2d-2
         
         max_model_number = 50000

         rate_factor_for_fe56ec_to_cr56 = 1d-4
         
         mesh_delta_coeff = 2         

         RGB_wind_scheme = 'Dutch'
         AGB_wind_scheme = 'Dutch'
         RGB_to_AGB_wind_switch = 1d-4
         Dutch_wind_eta = 0.8

         log_center_density_limit = 12 
         log_center_temp_limit = 10.3d0 ! max T about 2d10
         
         xa_scale = 1d-5
         newton_itermin = 2

         mass_change_full_on_dt = 1d7 ! (seconds)
         mass_change_full_off_dt = 1d6 ! (seconds)

         mixing_length_alpha = 1.6  ! 1.0 for M > 30
         screening_mode = 'extended'

         use_Ledoux_criterion = .true.
         alpha_semiconvection = 0.1

         photostep = 50
         profile_interval = 50
         log_cnt = 100
         terminal_cnt = 5
         write_header_frequency = 10
         
         !xa_central_upper_limit_species(1) = 'cr56'
         !xa_central_upper_limit(1) = 0.80
         
         max_iter_for_resid_tol1 = 3
         tol_residual_norm1 = 1d-5
         tol_max_residual1 = 1d-2
         
         max_iter_for_resid_tol2 = 12
         tol_residual_norm2 = 1d99
         tol_max_residual2 = 1d99
         
         min_timestep_limit = 1d-12 ! (seconds)

         !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_why_dt_limits = .true.
         !report_all_dt_limits = .true.
         
         !show_mesh_changes = .true.
         !mesh_dump_call_number = 1598
         
         !okay_to_remesh = .false.
         
         !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 = 195

         use_artificial_viscosity = .false.
         
         delta_lgL_He_limit = -1



         !dX_nuc_max_A_limit = 52
         !dX_limit_min_X = 1d-3
         !dX_limit = 2d-3         
         !dX_nuc_limit = 2d-3
         !dX_nuc_factor = 5d2
         !dX_nuc_min_yrs_for_dt = 5d-10


         dX_nuc_drop_max_A_limit = 52
         dX_nuc_drop_min_X_limit = 1d-4
         dX_nuc_drop_limit = 5d-3
         dX_nuc_drop_hard_limit = 1d99
         
         
         
         
         
         delta_lgTeff_limit = 0.005  
         delta_lgL_limit = 0.05
         delta_lgRho_cntr_limit = 0.02  




         mesh_dlog_pp_dlogP_extra = 0.15
         mesh_dlog_cno_dlogP_extra = 0.15
      
         mesh_dlog_3alf_dlogP_extra = 0.15
         mesh_dlog_burn_c_dlogP_extra = 0.15
         
         mesh_dlog_burn_n_dlogP_extra = 0.15
         mesh_dlog_burn_o_dlogP_extra = 0.15
      
         mesh_dlog_burn_ne_dlogP_extra = 0.15
         
         mesh_dlog_burn_na_dlogP_extra = 0.15
         mesh_dlog_burn_mg_dlogP_extra = 0.15
         
         mesh_dlog_cc_dlogP_extra = 0.15
         mesh_dlog_co_dlogP_extra = 0.15
         
         mesh_dlog_oo_dlogP_extra = 0.15
         
         mesh_dlog_burn_si_dlogP_extra = 0.15
         mesh_dlog_burn_s_dlogP_extra = 0.15
         mesh_dlog_burn_ar_dlogP_extra = 0.15
         mesh_dlog_burn_ca_dlogP_extra = 0.15
         mesh_dlog_burn_ti_dlogP_extra = 0.15
         mesh_dlog_burn_cr_dlogP_extra = 0.15
         mesh_dlog_burn_fe_dlogP_extra = 0.15
         
         mesh_dlog_pnhe4_dlogP_extra = 0.15
         mesh_dlog_other_dlogP_extra = 0.15
         mesh_dlog_photo_dlogP_extra = 1
         
         
         xa_function_species(1) = 'ne20'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'si28'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 's32'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'fe54'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'fe56'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2

 
         !overshoot_f_above_nonburn = 0.01
         !overshoot_f_below_nonburn = 0.01
         !overshoot_f_above_burn_h = 0.01
         !overshoot_f_below_burn_h = 0.01


         min_overshoot_q = 1d-3 

         T_mix_limit = 1d4


/ ! end of controls namelist




&pgstar
         
      ! top level controls

         !pause = .true. 
            ! if true, the code waits for user to enter a RETURN on the command line
      
      ! main window
      
         MAIN_win_flag = .true.
         
         show_TRho_Profile_with_main = .true.
      
         show_HR_TRho_with_main = .false.
         MAIN_win_width = 8.5
         MAIN_win_aspect_ratio = 1.2 ! aspect_ratio = height/width
         
         !show_HR_TRho_with_main = .false. ! if false, use aspect ratio > 1
         !MAIN_win_width = 7.5
         !MAIN_win_aspect_ratio = 1.2 ! aspect_ratio = height/width
      
         xaxis_by = 'by_mass'
            
         show_main_win_text_info = .false.
         
         ! xaxis limits -- to override system default selections
         !main_xmin = -10 !0.4
         main_xmax = 3
         
         logxq_cutoff = -12.1 ! min value when using logxm for xaxis
         !main_logL_max = 7
         !main_logL_min = 0
         
         ! control for plot showing abundances
         log_mass_frac_ymax = 0.7 ! making this > 0 helps readability
         log_mass_frac_ymin = -4.0
         num_abundance_line_labels = 5
      
         log_abund_vary = 0.1 ! don't show species if it would vary by < this in plot
      
         ! file output
         !MAIN_file_flag = .true.
         MAIN_file_dir = 'png'
         MAIN_file_prefix = 'main'
         MAIN_file_cnt = 4 ! output when mod(model_number,main_file_cnt)==0
               
               
      ! TRho Profile window -- current model in T-Rho plane
      
         TRho_Profile_win_flag = .false.

         TRho_Profile_win_width = 9.5
         TRho_Profile_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
            
         show_TRho_Profile_text_info = .true.
         show_TRho_Profile_legend = .true.
         show_HR_TRho_with_TRho_Profile = .false.
         show_TRho_Profile_burn_labels = .true.
         TRho_Profile_show_logQ_limit = .true.
         !show_TRho_Profile_cross_hair = .true.
         show_TRho_Profile_mass_locs = .false.

         !show_TRho_Profile_annotation1 = .true.
         !show_TRho_Profile_annotation2 = .true.
         !show_TRho_Profile_annotation3 = .true.
      
         ! axis limits
         TRho_Profile_xmin = -10
         TRho_Profile_xmax = 10
         TRho_Profile_ymin = 3.0
         TRho_Profile_ymax = 10
         
         ! file output
         !TRho_Profile_file_flag = .true.
         TRho_Profile_file_dir = 'png'
         TRho_Profile_file_prefix = 'trho_profile'
         TRho_Profile_file_cnt = 2 ! output when mod(model_number,TRho_Profile_file_cnt)==0
      
      
      ! Profile window
      
         Profile_win_flag = .false.

         Profile_win_width = 9.5
         Profile_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
            
         show_Profile_legend = .false.
         Profile_legend_coord = 0.65
         Profile_legend_fjust = 0.0
         Profile_legend_disp1 = -7.5
         Profile_legend_del_disp = -1.5
         Profile_show_mass_locs = .true.
         Profile_show_decorated_line = .true.
         
         show_HR_TRho_with_Profile = .false.
         show_Profile_text_info = .false.
         Profile_text_info_xfac = 0.78 ! controls x location
         Profile_text_info_dxfac = 0.02 ! controls x spacing to value from text
         Profile_text_info_yfac = 0.95 ! controls y location of 1st line
         Profile_text_info_dyfac = -0.04 ! controls line spacing
         show_Profile_cross_hair = .false.

         show_Profile_annotation1 = .false.
         show_Profile_annotation2 = .false.
         show_Profile_annotation3 = .false.
         
         ! axis choices
            ! to get a list of axis identifiers, edit your star_job controls
               ! to set show_profile_column_numbers = .true.
            ! you can use any identifier for x or y
         Profile_xaxis_name = 'logP'
         !Profile_xaxis_reversed = .true.
         !Profile_xmin = 0 ! only used if > -100
         !Profile_xmax = 1.2 !-101 ! only used if > -100
         
         Profile_yaxis_name = 'entropy' !'grada'
         Profile_yaxis_reversed = .false.
         Profile_ymin = -111 !5.4 ! only used if > -100
         Profile_ymax = -111 !8.6 ! only used if > -100    
         Profile_dymin = 0.025    
         
         Profile_other_yaxis_name = 'logQ'
         Profile_other_yaxis_reversed = .false.
         Profile_other_ymin = -101 ! only used if > -100
         Profile_other_ymax = -101 ! only used if > -100        
         
         ! file output
         !Profile_file_flag = .true.
         Profile_file_dir = 'png'
         Profile_file_prefix = 'profile'
         Profile_file_cnt = 5 ! output when mod(model_number,Profile_file_cnt)==0
      
      
      ! Convection window -- history of convection and more.
      
         !CONV_win_flag = .true.
      
         CONV_win_width = 9
         CONV_win_aspect_ratio = 0.62 ! aspect_ratio = height/width
      
         ! axis limits
         CONV_xmax = -1 ! step number.  negative means use default.
         CONV_xmin = -1 ! step number.  negative means use default.
         CONV_max_width = 0 ! only used if > 0.  causes xmin to move with xmax.
         CONV_mmax = 1.00007  !-1 ! (Msun units) negative means use start initial mass
         CONV_mmin = 1  !0 ! (Msun units)
         CONV_lgLmax = -101 ! only used if > -100; (L in Lsun units)
         CONV_lgLmin = -2 ! only used if > -100; (L in Lsun units)

         CONV_show_burn = .true.
         CONV_show_mixing = .true.
         CONV_show_log_radius = .true.
         CONV_show_luminosities = .true.
         CONV_show_mass_boundaries = .true.

         show_CONV_annotation1 = .false.
         show_CONV_annotation2 = .false.
         show_CONV_annotation3 = .false.
         
         ! file output
         CONV_file_flag = .false.
         CONV_file_dir = 'pgstar_out'
         CONV_file_prefix = 'conv'
         CONV_file_cnt = 5 ! output when mod(model_number,CONV_file_cnt)==0
         CONV_file_width = -1 ! negative means use same value as for window
         CONV_file_aspect_ratio = -1 ! negative means use same value as for window
      
      
         


/ ! end of pgstar namelist
