
! inlist_wd2


&star_job
      show_log_description_at_start = .false.
      
      load_saved_model = .true.
      saved_model_name = 'wd2.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 = 'cno_extras.net'

      change_v_flag = .true.
      new_v_flag = .false.
      
      set_tau_factor = .true.
      set_to_this_tau_factor = 300

      !pgstar_flag = .true.

      !auto_turn_on_S_surf_when_adding_mass = .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 = 100
      max_number_retries = 800
      
      max_model_number = 1000

      use_Ledoux_criterion = .false.
      
      use_Type2_opacities = .true.
      Zbase = 0.02d0
         
      logQ_limit = 100
      gamma_center_limit = 1000 ! FOR THIS TEST
            
      star_mass_max_limit = 1.00005d0
         
      num_trace_history_values = 5
      trace_history_value_name(1) = 'k_below_just_added'
      trace_history_value_name(2) = 'logxq_below_just_added'
      trace_history_value_name(3) = 'k_const_mass'
      trace_history_value_name(4) = 'logxq_const_mass'
      trace_history_value_name(5) = 'split_mixing_choice'

      max_center_cell_dq = 1d-6
      max_surface_cell_dq = 1d-8

      photostep = 10
      profile_interval = 10
      history_interval = 1
      terminal_cnt = 1
      write_header_frequency = 10
      
      min_timestep_limit = 1d-12
      
      accrete_same_as_surface = .false. 
      accretion_h1 = 0.749D+00
      accretion_h2 = 0
      accretion_he3 = 0.292910D-04
      accretion_he4 = 0.237D+00
      accretion_zfracs = 4 ! Lodders 03

      xa_function_species(1) = 'h1'  ! name of nuclide as defined in chem_def
		xa_function_weight(1) = 20
		xa_function_param(1) = 1d-6
      xa_function_species(2) = 'he4'  ! name of nuclide as defined in chem_def
		xa_function_weight(2) = 20
		xa_function_param(2) = 1d-2

      dH_div_H_limit = 0.8d0   ! 10.7 101 3 2137

! 6.2 30 0 867

!      delta_lgL_H_limit = 0.03 ! limit for magnitude of change in lgL_H
      delta_lgL_H_limit = 0.1 ! limit for magnitude of change in lgL_H
      delta_lgL_H_hard_limit = 0.5
      lgL_H_burn_min = 1.5 ! ignore changes in lgL_H if value is less than this

      delta_lgL_He_limit = 0.025 ! limit for magnitude of change in lgL_He
      delta_lgL_He_hard_limit = 1
      lgL_He_burn_min = 2.5 ! ignore changes in lgL_He if value is less than this

      delta_lgT_limit = 0.5  ! limit for magnitude of max change in log10 temperature at any point
      delta_lgT_hard_limit = 1
      
      delta_lgRho_limit = 1 ! limit for magnitude of max change in log10 density at any point
      delta_lgRho_hard_limit = -1

      relax_hard_limits_after_retry = .false.


      delta_lgTeff_limit = 0.05  
         ! limit for magnitude of max change in log10 temperature at photosphere   
      delta_lgTeff_hard_limit = 0.1
   
      ! limit for magnitude of change in log10(L/Lsun)
      delta_lgL_limit = 0.05
      delta_lgL_hard_limit = 0.15

      varcontrol_target = 1d-2
      
      super_eddington_wind_eta = 1
      super_eddington_wind_Ledd_factor = 5
      
      net_logTcut_lo = 0
      net_logTcut_lim = 0
      
      which_atm_option = 'grey_and_kap'
      Pextra_factor = -1
      
      !use_artificial_viscosity = .true.
      !l1_coef = 100 ! increase to as much as 1d3 to suppress pulsations

      mass_change = 2.5d-7 ! 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.

      !trace_newton_bcyclic_solve_input = .true. ! input is "B" j k iter B(j,k)
      !trace_newton_bcyclic_solve_output = .true. ! output is "X" j k iter X(j,k)
      
      trace_newton_bcyclic_steplo = 1 ! 1st model number to trace
      trace_newton_bcyclic_stephi = 1 ! last model number to trace
      
      trace_newton_bcyclic_iterlo = 2 ! 1st newton iter to trace
      trace_newton_bcyclic_iterhi = 2 ! last newton iter to trace
      
      trace_newton_bcyclic_nzlo = 1 ! 1st cell to trace
      trace_newton_bcyclic_nzhi = 10000 ! last cell to trace; if < 0, then use nz as nzhi
      
      trace_newton_bcyclic_jlo = 1 ! 1st var to trace
      trace_newton_bcyclic_jhi = 100 ! last var to trace; if < 0, then use nvar as jhi

      ! 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
