! inlist_wd_3
! accrete until mg capture heating kicks in


&star_job
      show_log_description_at_start = .false.

      mesa_dir = '../../..'
      
      load_saved_model = .true.
      saved_model_name = 'wd_2.mod'

      save_model_when_terminate = .true.
      save_model_filename = 'wd_3.mod'

      set_tau_factor = .true.
      set_to_this_tau_factor = 300

      ! pgstar_flag = .true.

      set_initial_dt = .true.
      years_for_initial_dt = 10 ! in years

!      relax_mass_change = .true.
!      relax_mass_change_min_steps = 25
!      relax_mass_change_max_yrs_dt = 1 ! years
!      relax_mass_change_init_mdot = 1d-8 ! Msun/year
!      relax_mass_change_final_mdot = 1d-6 ! Msun/year
      
/ ! 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
         
         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'

      initial_mass =  1.3
      initial_z = 0.02d0

      varcontrol_target = 1e-3
      
      mdot_eps_grav_limit = -1
      delta_lgTeff_limit = -1  

      dX_div_X_limit_min_X = 0.01
      dX_div_X_limit = 0.05

      which_atm_option = 'grey_and_kap'

      log_center_temp_limit = 8.5 ! Tc about 3e8
      log_center_density_limit = 9.7

      mass_change = 1d-7 ! 1d-6 ! rate of accretion (Msun/year)

      accrete_same_as_surface = .false. 

      accrete_given_mass_fractions = .true. 
      num_accretion_species = 2
      accretion_species_id(1) = 'c12'
      accretion_species_xa(1) = 0.5
      accretion_species_id(2) = 'o16'
      accretion_species_xa(2) = 0.5

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


/ ! end of controls namelist



&pgstar

/ ! end of pgstar namelist
