! inlist_ns


&star_job
      show_log_description_at_start = .false.
   
      load_saved_model = .true.
      saved_model_name = 'ns_1.4M_c.mod'
      !saved_model_name = 'test.mod'

      save_model_when_terminate = .true.
      !save_model_number = 51
      save_model_filename = 'final.mod'
      !save_model_filename = 'test.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'

      ! shut off PC since get bad convergence for high gamma
      set_eos_PC_parameters = .true.
      mass_fraction_limit_for_PC = 1d-2 ! skip species if abundance < this
      logRho1_PC_limit = 20 ! 3.7d0 ! okay for pure PC for logRho > this
      logRho2_PC_limit = 20 ! 2.8d0 ! don't use PC for logRho < this (>= 2.8 or so because of PPT)
      log_Gamma_all_HELM = 1.60206d0 ! HELM for log_Gamma <= this
      log_Gamma_all_PC = 1.90309d0 ! PC for log_Gamma >= this

      set_to_this_tau_factor = 300 ! tau_factor of 1 puts outer cell at photosphere.
      set_tau_factor = .true. ! change tau_factor without reconverging.

      kappa_file_prefix = 'gs98'

      set_initial_age = .true.
      initial_age = 0 ! in years

      set_initial_model_number = .true.
      initial_model_number = 0

      !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 = 20
         max_number_retries = 80
         max_model_number = 2000

      initial_mass = 1
      initial_z = 0.02d0


      max_abar_for_burning = 55
      
      use_Type2_opacities = .true.
      Zbase = 0.02d0
      
      photostep = 10
      profile_interval = 10
      history_interval = 1
      terminal_cnt = 1
      write_header_frequency = 10
      
      min_timestep_limit = 1d-9
      gamma_center_limit = 1d99

      power_nuc_burn_upper_limit = 6d4
      
      merge_if_dlnR_too_small = .true.
      mesh_min_dlnR = 1d-9

      
      which_atm_option = 'grey_and_kap'
      !mesh_delta_coeff = 0.5
      
      varcontrol_target = 1d-3 
      
      RGB_wind_scheme = ''
      AGB_wind_scheme = ''
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.5d0  
      Blocker_wind_eta = 0.1d0 
      
      use_GR_factors = .true.
      
      !super_eddington_wind_eta = 1
         
      !mass_change_full_on_dt = 1d-4 ! (seconds)
      !mass_change_full_off_dt = 1d-6 ! (seconds)

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

!      net_pp_rate_factor = 0
!      net_cno_rate_factor = 0
!      max_abar_for_burning = 55
      
      mass_change = 1.7d-9  ! rate of accretion (Msun/year)

      ! 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
         
      ! top level controls

         !pause = .true. 
            ! if true, the code waits for user to enter a RETURN on the command line
      


/ ! end of pgstar namelist
