
! inlist_wd3


&star_job
      show_log_description_at_start = .false.
      
      load_saved_model = .true.
      saved_model_name = 'wd3.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 = '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 = 300 !5e5

      !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 = 190
      
      max_model_number = 2000
      
      use_Type2_opacities = .true.
      Zbase = 0.02d0
         
      num_trace_history_values = 0
      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_center_cell_dq = 1d-6
      max_surface_cell_dq = 1d-8
      
      photostep = 50
      profile_interval = 500
      history_interval = 50
      terminal_cnt = 5
      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.8 ! 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)


      ! 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
