&star_job  ! inlist_astero   

      mesa_dir = '' ! empty string defaults to environment variable MESA_DIR

      show_log_description_at_start = .false.
      
      load_saved_model = .true.
      saved_model_name = 'astero_test.mod'

      change_net = .true. ! switch nuclear reaction network
      new_net_name = 'o18_and_ne22.net'

	   set_rate_n14pg = 'Imbriani'
	   set_rate_c12ag = 'Kunz'
      kappa_file_prefix = 'a09'
      
      change_lnPgas_flag = .true.
      new_lnPgas_flag = .true.

      !pgstar_flag = .true.


      set_to_this_tau_factor = 1d-4
      set_tau_factor = .true.

      

/ ! end of star_job namelist


&controls

      calculate_Brunt_N2 = .true. 
      num_cells_for_smooth_brunt_B = 0

      which_atm_option = 'simple_photosphere'    
      use_Henyey_MLT = .true.
      
      
      !max_model_number = 1295
       

      ! timesteps
      max_timestep = 3.15d15  ! in seconds.  max_timestep <= 0 means no upper limit.            
      varcontrol_target = 1d-4

      ! mesh adjustment
      mesh_delta_coeff = 0.5

      ! controls for output
      photostep = 10
      profile_interval = 10
      history_interval = 5
      terminal_interval = 5
      write_header_frequency = 10 


/ ! 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


      ! dPg_dnu window -- delta_Pg vs. delta_nu
      
         !dPg_dnu_win_flag = .true.
      
         dPg_dnu_win_width = 9
         dPg_dnu_win_aspect_ratio = 0.62 ! aspect_ratio = height/width
         
         show_dPg_dnu_cross_hair = .false.
      
         ! axis limits -- to override system default selections
         dPg_dnu_delta_nu_min = -101 ! only used if > -100
         dPg_dnu_delta_nu_max = -101 ! only used if > -100
         dPg_dnu_delta_Pg_min = -101 ! only used if > -100
         dPg_dnu_delta_Pg_max = -101 ! only used if > -100
      
         dPg_dnu_age_min = -1 ! only plot models with age >= this
         dPg_dnu_step_min = -1 ! only plot models with model number >= this
         dPg_dnu_step_max = 999999 ! only plot models with model number <= this

         show_dPg_dnu_annotation1 = .false.
         show_dPg_dnu_annotation2 = .false.
         show_dPg_dnu_annotation3 = .false.
         
         ! file output
         dPg_dnu_file_flag = .false.
         dPg_dnu_file_dir = 'png'
         dPg_dnu_file_prefix = 'dPg_dnu'
         dPg_dnu_file_cnt = 5 ! output when mod(model_number,dPg_dnu_file_cnt)==0
         dPg_dnu_file_width = -1 ! negative means use same value as for window
         dPg_dnu_file_aspect_ratio = -1 ! negative means use same value as for window

      
      
      ! HR window -- lg_L vs. lg_Teff
      
         HR_win_flag = .false.
      
         HR_win_width = 9
         HR_win_aspect_ratio = 0.62 ! aspect_ratio = height/width
         
         show_HR_cross_hair = .false.
      
         ! axis limits -- to override system default selections
         HR_logT_min = 3.74 ! -101 ! only used if > -100
         HR_logT_max = 3.83 ! only used if > -100
         HR_logL_min = 0.25 ! -101 ! only used if > -100
         HR_logL_max = 0.82 ! -101 ! only used if > -100
      
         HR_age_min = -1 ! only plot models with age >= this
         HR_step_min = -1 ! only plot models with model number >= this
         HR_step_max = 999999 ! only plot models with model number <= this

         show_HR_annotation1 = .false.
         show_HR_annotation2 = .false.
         show_HR_annotation3 = .false.
         
         ! file output
         HR_file_flag = .false.
         HR_file_dir = 'pgstar_out'
         HR_file_prefix = 'hr'
         HR_file_cnt = 5 ! output when mod(model_number,HR_file_cnt)==0
         HR_file_width = -1 ! negative means use same value as for window
         HR_file_aspect_ratio = -1 ! negative means use same value as for window
         
      
               
               
      ! TRho Profile window -- current model in T-Rho plane
      
         TRho_Profile_win_flag = .false.

         TRho_Profile_win_width = 7.5
         TRho_Profile_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
            
         show_TRho_Profile_text_info = .false.
         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_kap_regions = .false.
         show_TRho_Profile_eos_regions = .true.
         show_TRho_Profile_degeneracy_line = .true.
         show_TRho_Profile_Pgas_Prad_line = .true.
         show_TRho_Profile_burn_lines = .true.

         show_TRho_Profile_annotation1 = .true.
         annotation1_ci = 1 ! color index
         annotation1_ch = 1.1 ! character height
         annotation1_lw = 1 ! line width
         annotation1_cf = 1 ! character font (1=sans serif, 2=roman, 3=italic, 4=script)
         annotation1_side = 'T' ! T, B, L, or R
         annotation1_disp = 1.8 ! displacement perpendicular to axis
         annotation1_coord = 0.5 ! fractional location along axis
         annotation1_fjust = 0.5 ! 0 => left justified, 0.5 => centered, 1 => right justified
         annotation1_text = 'HD49385'
      
         ! 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
         
         
      ! Mode Propagation window (for asterseismology)
      
         !Mode_Prop_win_flag = .true.
      
         Mode_Prop_win_width = 8
         Mode_Prop_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
         
         Mode_Prop_nu_max_obs = 1000
      
         Mode_Prop_xaxis_by = 'by_r_div_R' ! select xaxis (same options as for main window)
         
         ! power xaxis limits -- to override system default selections
         Mode_Prop_xmin = -101 ! only used if > -100
         Mode_Prop_xmax = 0.3 !-101 ! only used if > -100
         
         ! power yaxis limits -- to override system default selections
         Mode_Prop_ymin = -101 ! only used if > -100
         Mode_Prop_ymax = -101 ! only used if > -100

         show_Mode_Prop_cross_hair = .false.
         
         ! file output
         Mode_Prop_file_flag = .false.
         Mode_Prop_file_dir = 'pgstar_out'
         Mode_Prop_file_prefix = 'mode_prop'
         Mode_Prop_file_cnt = 5 ! output when mod(model_number,Mode_Prop_file_cnt)==0
         Mode_Prop_file_width = -1 ! negative means use same value as for window
         Mode_Prop_file_aspect_ratio = -1 ! negative means use same value as for window      
      
      
      ! loggTe window -- logg vs. Teff
      
         loggTe_win_flag = .false.
      
         loggTe_win_width = 7
         loggTe_win_aspect_ratio = 0.62 ! aspect_ratio = height/width         
         
         show_loggTe_cross_hair = .false.


         show_loggTe_target_box = .true.
         loggTe_target_n_sigma = -3
         loggTe_target_logg = 4.00d0
         loggTe_target_logg_sigma = 0.06d0
         loggTe_target_Teff = 6095d0
         loggTe_target_Teff_sigma = 65  
      
         ! axis limits -- to override system default selections
         loggTe_logg_min = 3.8 ! only used if > -100
         loggTe_logg_max = 4.4 ! only used if > -100
         loggTe_Teff_min = 5700 ! only used if > -100
         loggTe_Teff_max = 6500 ! only used if > -100
      
         loggTe_age_min = -1 ! only plot models with age >= this
         loggTe_step_min = -1 ! only plot models with model number >= this
         loggTe_step_max = 999999 ! only plot models with model number <= this

         show_loggTe_annotation1 = .false.
         show_loggTe_annotation2 = .false.
         show_loggTe_annotation3 = .false.
         
         ! file output
         loggTe_file_flag = .false.
         loggTe_file_dir = 'pgstar_out'
         loggTe_file_prefix = 'loggTe'
         loggTe_file_cnt = 5 ! output when mod(model_number,loggTe_file_cnt)==0
         loggTe_file_width = -1 ! negative means use same value as for window
         loggTe_file_aspect_ratio = -1 ! negative means use same value as for window
      
      
      ! Profile window
      
         Profile_win_flag = .false.

         Profile_win_width = 7.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 = .false.
         Profile_show_decorated_line = .false.
         
         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 = 'logxm'
         Profile_xaxis_reversed = .true.
         Profile_xmin = -12 ! 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 = 'logL'
         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 = 8
         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 ! (Msun units) negative means use start initial mass
         CONV_mmin = 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 = .true.
         CONV_file_dir = 'png_extra'
         CONV_file_prefix = 'conv'
         CONV_file_cnt = 1 ! 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
      
      
         
         
      ! Surface History window
      
         Surf_Hist_win_flag = .false.
      
         Surf_Hist_win_width = 6.5
         Surf_Hist_win_aspect_ratio = 1.2 ! aspect_ratio = height/width
      
         ! axis limits
         Surf_Hist_xmax = -1 ! step number.  negative means use default.
         Surf_Hist_xmin = -1 ! step number.  negative means use default.
         Surf_Hist_max_width = 0 ! only used if > 0.  causes xmin to move with xmax.
         
         ! abundance limits
         Surf_Hist_mass_frac_min = 0
         Surf_Hist_mass_frac_max = 1.1 ! > 0 helps legibility

         show_Surf_Hist_annotation1 = .false.
         show_Surf_Hist_annotation2 = .false.
         show_Surf_Hist_annotation3 = .false.
         
         ! file output
         Surf_Hist_file_flag = .false.
         Surf_Hist_file_dir = 'pgstar_out'
         Surf_Hist_file_prefix = 'conv'
         Surf_Hist_file_cnt = 5 ! output when mod(model_number,Surf_Hist_file_cnt)==0
         Surf_Hist_file_width = -1 ! negative means use same value as for window
         Surf_Hist_file_aspect_ratio = -1 ! negative means use same value as for window
      
         
         
      ! Power window
      
         Power_win_flag = .false.
      
         Power_win_width = 8.5
         Power_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
      
         Power_xaxis_by = 'by_mass'
         
         ! power xaxis limits -- to override system default selections
         Power_xmin = -101 ! only used if > -100
         Power_xmax = -101 ! only used if > -100
         
         ! power yaxis limits -- to override system default selections
         Power_ymin = -101 ! only used if > -100
         Power_ymax = -101 ! only used if > -100

         show_Power_cross_hair = .false.
         
         ! file output
         Power_file_flag = .false.
         Power_file_dir = 'pgstar_out'
         Power_file_prefix = 'power'
         Power_file_cnt = 5 ! output when mod(model_number,Power_file_cnt)==0
         Power_file_width = -1 ! negative means use same value as for window
         Power_file_aspect_ratio = -1 ! negative means use same value as for window
         
         
      ! Abundance window
      
         Abundance_win_flag = .false.
      
         Abundance_win_width = 8.5
         Abundance_win_aspect_ratio = 0.618 ! aspect_ratio = height/width
      
         Abundance_xaxis_by = 'by_mass' ! select xaxis
         
         ! power xaxis limits -- to override system default selections
         Abundance_xmin = -101 ! only used if > -100
         Abundance_xmax = -101 ! only used if > -100
         
         ! power yaxis limits -- to override system default selections
         Abundance_log_mass_frac_min = 1 ! only used if < 0
         Abundance_log_mass_frac_max = 1 ! only used if < 0

         show_Abundance_cross_hair = .false.
         
         ! file output
         Abundance_file_flag = .false.
         Abundance_file_dir = 'pgstar_out'
         Abundance_file_prefix = 'abund'
         Abundance_file_cnt = 5 ! output when mod(model_number,Abundance_file_cnt)==0
         Abundance_file_width = -1 ! negative means use same value as for window
         Abundance_file_aspect_ratio = -1 ! negative means use same value as for window
         
         


/ ! end of pgstar namelist
