
! inlist_wd_ignite_1
! remove small remaining helium envelope


&star_job
      show_log_description_at_start = .false.

      read_extra_star_job_inlist1 = .true.
      extra_star_job_inlist1_name = '../inlist_test_suite'
      
      load_saved_model = .true.
      saved_model_name = '../../../data/star_data/white_dwarf_models/1.316_from_8.5_z2m2.mod'

      save_model_when_terminate = .true.
      save_model_filename = 'wd_1.mod'

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

      change_v_flag = .true.
      new_v_flag = .true.

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

      relax_tau_factor = .true. ! gradually change tau_factor, reconverging at each step.
      relax_to_this_tau_factor = 300 ! tau_factor of 1 puts outer cell at photosphere.
      dlogtau_factor = 0.2 ! 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 = 10
         max_number_retries = 80

         max_model_number = 1000

      initial_mass =  1.3
      initial_z = 0.02d0
      
      use_Type2_opacities = .true.
      Zbase = 0.02d0
      
      use_Ledoux_criterion = .true.
      alpha_semiconvection = 4d-2
      thermohaline_coeff = 2  
      
      xa_surface_lower_limit_species(1) = 'he4'
      xa_surface_lower_limit(1) = 1d-8
      
      varcontrol_target = 1d-2
      
      delta_lgTeff_limit = -1  

      logQ_limit = 1d6
      gamma_center_limit = 1d6

      photostep = 10
      profile_interval = 10
      history_interval = 1
      terminal_cnt = 1
      write_header_frequency = 10

         
      ! DISABLE THESE
      !sum_xa_tolerance = 1d99
      !min_xa_hard_limit = -1d99 
               
      xa_scale = 1d-5
         
      tol_residual_norm1 = 1d-5
      tol_max_residual1 = 1d-2
      iter_for_resid_tol2 = 3
      tol_residual_norm2 = 1d99
      tol_max_residual2 = 1d99
      iter_for_resid_tol3 = 12

      min_timestep_limit = 1d-4 ! seconds

         dX_nuc_drop_max_A_limit = 52
         dX_nuc_drop_min_X_limit = 1d-4
         dX_nuc_drop_limit = 5d-3
         dX_nuc_drop_hard_limit = 1d99

         mesh_dlog_pp_dlogP_extra = 0.15
         mesh_dlog_cno_dlogP_extra = 0.15
      
         mesh_dlog_3alf_dlogP_extra = 0.15
         mesh_dlog_burn_c_dlogP_extra = 0.15
         
         mesh_dlog_burn_n_dlogP_extra = 0.15
         mesh_dlog_burn_o_dlogP_extra = 0.15
      
         mesh_dlog_burn_ne_dlogP_extra = 0.15
         
         mesh_dlog_burn_na_dlogP_extra = 0.15
         mesh_dlog_burn_mg_dlogP_extra = 0.15
         
         mesh_dlog_cc_dlogP_extra = 0.15
         mesh_dlog_co_dlogP_extra = 0.15
         
         mesh_dlog_oo_dlogP_extra = 0.15
         
         mesh_dlog_burn_si_dlogP_extra = 0.15
         mesh_dlog_burn_s_dlogP_extra = 0.15
         mesh_dlog_burn_ar_dlogP_extra = 0.15
         mesh_dlog_burn_ca_dlogP_extra = 0.15
         mesh_dlog_burn_ti_dlogP_extra = 0.15
         mesh_dlog_burn_cr_dlogP_extra = 0.15
         mesh_dlog_burn_fe_dlogP_extra = 0.15
         
         mesh_dlog_pnhe4_dlogP_extra = 0.15
         mesh_dlog_other_dlogP_extra = 0.15
         mesh_dlog_photo_dlogP_extra = 1
         
         xa_function_species(1) = 'ne20'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'si28'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 's32'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'fe54'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2
         
         xa_function_species(1) = 'fe56'
   		xa_function_weight(1) = 20
   		xa_function_param(1) = 1d-2

 
         overshoot_f_above_nonburn_core = 0.01
         overshoot_f_above_nonburn_shell = 0.01
         overshoot_f_below_nonburn_shell = 0.01
         overshoot_f_above_burn_h_core = 0.01
         overshoot_f_above_burn_h_shell = 0.01
         overshoot_f_below_burn_h_shell = 0.01
 
         overshoot_f0_above_nonburn_core = 0.005
         overshoot_f0_above_nonburn_shell = 0.005
         overshoot_f0_below_nonburn_shell = 0.005
         overshoot_f0_above_burn_h_core = 0.005
         overshoot_f0_above_burn_h_shell = 0.005
         overshoot_f0_below_burn_h_shell = 0.005

         min_overshoot_q = 1d-3 

      T_mix_limit = 1d4
      
      which_atm_option = 'grey_and_kap'

      accrete_same_as_surface = .false. 
         accrete_given_mass_fractions = .true. 
         num_accretion_species = 1
         accretion_species_id(1) = 'c12'
         accretion_species_xa(1) = 1
         accretion_species_id(2) = 'o16'
         accretion_species_xa(2) = 0
      
      mass_change = -1d-7 ! 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
