
! inlist_wd3


&star_job
      show_log_description_at_start = .false.
      
      load_saved_model = .true.
      saved_model_name = 'wd3.mod'

      change_net = .true.      
      new_net_name = 'co_burn.net'
      
      set_rate_3a = 'FL87' ! Fushiki and Lamb, Apj, 317, 368-388, 1987
      set_rate_1212 = 'G05' ! Gasques, et al. Phys Review C, 72, 025806, 2005

      change_v_flag = .true.
      new_v_flag = .false.

      set_Z_all_HELM = .true. ! if false, just use the eos default
      Z_all_HELM = 1.000001d0

      set_tau_factor = .true. ! the test model was saved with this tau factor
      set_to_this_tau_factor = 5e5

      relax_tau_factor = .true. ! gradually change tau_factor, reconverging at each step.
      relax_to_this_tau_factor = 300 
      dlogtau_factor = 0.1 ! change log10(tau_factor) by at most this amount per step

      !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 = 30
         max_number_retries = 90
         
      num_trace_history_values = 4
      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'
         
      okay_to_reduce_gradT_excess = .true.
         
            ! these are for calculation of efficiency boosted gradT
            gradT_excess_f1 = 1d-4
            gradT_excess_f2 = 1d-2
            gradT_excess_age_fraction = 0.9d0
            
            ! turn on full all the time
            gradT_excess_lambda1 = 1.0
            gradT_excess_beta1 = 0.5
            gradT_excess_lambda2 = 1.0
            gradT_excess_beta2 = 0.5
            gradT_excess_dlambda = 1
            gradT_excess_dbeta = 1


      gamma_center_limit = 1000 ! FOR THIS TEST
      
      power_he_burn_upper_limit = 1d11
      
      max_model_number = 2000

      min_center_cell_dq = 1d-7
      max_center_cell_dq = 1d-6
      max_surface_cell_dq = 1d-8
      
      photostep = 100
      profile_interval = 100
      history_interval = 10
      terminal_cnt = 10
      write_header_frequency = 10
      
      min_timestep_limit = 1d-12
      
      accrete_same_as_surface = .false. 
      accrete_given_mass_fractions = .true. 
      num_accretion_species = 1
      accretion_species_id(1) = 'he4'
      accretion_species_xa(1) = 1

      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

      delta_lgL_H_limit = 0.05 ! 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_lgL_limit = 0.05
      delta_lgL_hard_limit = 0.15

      delta_lgT_limit = 0.5
      delta_lgT_hard_limit = 1
      
      delta_lgRho_limit = 1
      delta_lgRho_hard_limit = -1

      delta_lgTeff_limit = 0.1  
      delta_lgTeff_hard_limit = 0.3

      relax_hard_limits_after_retry = .false.


		mesh_delta_coeff = 0.4 !1.2

      varcontrol_target = 1d-2
      
      super_eddington_wind_eta = 1
      super_eddington_wind_Ledd_factor = 5
      
      net_logTcut_lo = 0
      net_logTcut_lim = 0
   	
      mass_change = 1d-9 ! rate of accretion (Msun/year)

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