! inlist_ns_he


&star_job
      
      load_saved_model = .true.
      saved_model_name = 'ns_1.4M.mod'

      save_model_when_terminate = .true.
      save_model_filename = 'final.mod'
   
      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'final_profile.data'

      change_initial_net = .true.      
      new_net_name = 'approx21_plus_co56.net'

      kappa_file_prefix = 'gs98'
      
      !pgstar_flag = .true.



/ ! end of star_job namelist

&controls

      ! check for retries and backups as part of test_suite
      ! you can/should delete this for use outside of test_suite
         max_number_backups = 150
         max_number_retries = 500

         max_model_number = 4000


      max_abar_for_burning = 55
      
      initial_mass = 1.4
      initial_z = 0.02d0
      
      use_Type2_opacities = .true.
      Zbase = 0.02d0
      
      photostep = 50
      profile_interval = 100
      history_interval = 100
      terminal_cnt = 10
      write_header_frequency = 10
      
		max_age = 0.05
      min_timestep_limit = 1d-9
      gamma_center_limit = 1d99
      
      varcontrol_target = 1d-3
      
      mesh_min_dlnR = 1d-9
      merge_if_dlnR_too_small = .true.
      
      which_atm_option = 'simple_photosphere'  
      Pextra_factor = 2

      accrete_same_as_surface = .false. 
      accrete_given_mass_fractions = .true. 
      num_accretion_species = 1
      accretion_species_xa(1) = 1
      accretion_species_id(1) = 'he4'

      mass_change = 3d-9  ! rate of accretion (Msun/year)
		okay_to_reduce_gradT_excess = .true.
		super_eddington_wind_eta = 1

      mixing_length_alpha = 1.5
      MLT_option = 'Henyey'

      use_Ledoux_criterion = .false.
      
      alpha_semiconvection = 0.1
      thermohaline_coeff = 2
      thermohaline_option = 'Kippenhahn'

      mesh_dlog_pp_dlogP_extra = 1
      mesh_dlog_cno_dlogP_extra = 1
   
      mesh_dlog_3alf_dlogP_extra = 1
      mesh_dlog_burn_c_dlogP_extra = 1
      mesh_dlog_burn_n_dlogP_extra = 1
      mesh_dlog_burn_o_dlogP_extra = 1
   
      mesh_dlog_burn_ne_dlogP_extra = 1
      mesh_dlog_burn_na_dlogP_extra = 1
      mesh_dlog_burn_mg_dlogP_extra = 1
      
      mesh_dlog_cc_dlogP_extra = 1
      mesh_dlog_co_dlogP_extra = 1
      mesh_dlog_oo_dlogP_extra = 1
      
      mesh_dlog_burn_si_dlogP_extra = 1
      mesh_dlog_burn_s_dlogP_extra = 1
      mesh_dlog_burn_ar_dlogP_extra = 1
      mesh_dlog_burn_ca_dlogP_extra = 1
      mesh_dlog_burn_ti_dlogP_extra = 1
      mesh_dlog_burn_cr_dlogP_extra = 1
      mesh_dlog_burn_fe_dlogP_extra = 1
      
      mesh_dlog_pnhe4_dlogP_extra = 1
      mesh_dlog_other_dlogP_extra = 1
      mesh_dlog_photo_dlogP_extra = 1





      ! FOR DEBUGGING

      !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
      
      !max_years_for_timestep = 3.67628942044319d-05

      !report_why_dt_limits = .true.
      !report_all_dt_limits = .true.
      
      !show_mesh_changes = .true.
      !mesh_dump_call_number = 5189
      !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


/ ! end of controls namelist









&pgstar
         


/ ! end of pgstar namelist
