
&star_job ! inlist_2

      show_log_description_at_start = .false.
      
      !create_pre_main_sequence_model = .true.

      save_model_when_terminate = .true.
      save_model_filename = 'astero_test.mod'
      
      load_saved_model = .true.
      saved_model_name = 'pre_ms.mod'

      change_initial_net = .true. ! switch nuclear reaction network
      new_net_name = 'o18_and_ne22.net'
	   
	   initial_zfracs = 6 ! AGSS09_zfracs = 6
	   !initial_zfracs = 5 ! AGS04_zfracs = 5
      set_uniform_initial_composition = .true.
      ! initial Z = 0.019
      initial_h1 = 0.741
      ! initial Y = 0.24
      initial_he3 = 3e-5
      initial_he4 = 0.23997
      initial_h2 = 0

      kappa_file_prefix = 'a09'

      !relax_initial_mass = .true.
      new_mass = 1.26866
      lg_max_abs_mdot = -100 

      !set_tau_factor = .true.
      !set_to_this_tau_factor = 15

      astero_just_call_my_extras_check_model = .true.

      !pgstar_flag = .true.
      

/ ! end of star_job namelist


&controls

      initial_mass = 1.21
      initial_z = 0.0207305

      mixing_length_alpha = 2.0d0

      log_L_lower_limit = 0.4

      ! atmosphere << run pre-ms model 100 steps before change this
      which_atm_option = 'photosphere_tables'     
      !which_atm_option = 'grey_and_kap'
      !which_atm_option = 'Eddington_grey'


      ! controls for output
      photostep = 50
      profile_interval = 50
      history_interval = 1
      terminal_cnt = 10
      write_header_frequency = 10
      
      !max_model_number = 3
      
      !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 = .false.
      report_all_dt_limits = .false.
      write_varcontrol_internals = -1
      
      report_hydro_dt_info = .false.
      
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 79

      trace_evolve = .false.
      
      ! hydro debugging
      !hydro_numerical_jacobian = .true.
      !hydro_check_everything = .true.
      !hydro_show_rcond = .true.
      !hydro_inspectB_flag = .true.
      !hydro_save_numjac_plot_data = .true.
      !hydro_dump_call_number = 1


      
      
      
      
      ! element diffusion
      !show_diffusion_info = .true. ! terminal output for diffusion
      !diffusion_dump_call_number = 228

      diffusion_dt_limit = 7d11 ! no element diffusion if dt < this limit (in seconds)
      diffusion_atol = 1d-5
      diffusion_rtol = 1d-4
      diffusion_T_full_on = 1d3
      diffusion_T_full_off = 1d3
      diffusion_calculates_ionization = .true.
      diffusion_num_classes = 3 ! number of classes of species for diffusion calculations
      diffusion_class_representative(1) = 'h1'
      diffusion_class_representative(2) = 'he4'
      diffusion_class_representative(3) = 'fe56'
      ! in ascending order.  species goes into 1st class with A_max >= species A
      diffusion_class_A_max(1) = 2
      diffusion_class_A_max(2) = 4
      diffusion_class_A_max(3) = 9999

      
      





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


      calculate_Brunt_N2 = .true. 
      brunt_N2_coefficient = 1
      

      ! mesh adjustment
      mesh_delta_coeff = 0.4



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

/ ! end of pgstar namelist
