! inlist_c13_pocket
! start at 53150 and end at 53350

&star_job

      show_log_description_at_start = .false.
      
      load_saved_model = .true.
      saved_model_name =  'AGBpreDUP2M.mod'

      save_model_when_terminate = .true.
      save_model_filename = 'final.mod'
      
      write_profile_when_terminate = .true.
      filename_for_profile_when_terminate = 'final_profile.data'

      set_rates_preference = .true.
      new_rates_preference = 1          
      set_rate_c12ag = 'Kunz'


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






      !redo_mix_info_at_each_iteration = .true.








      initial_mass = 2.00
      initial_z = 0.010

      mixing_length_alpha = 1.73

! some logistics
      max_model_number = 5000 ! negative means no maximum

      photostep = 10
      profile_interval = 10
      history_interval = 1
      terminal_cnt = 1
      write_header_frequency = 10
      
      max_num_profile_models = 1000 ! maximum number of saved profiles

! mesh
      max_allowed_nz = 5000 ! maximum number of grid points allowed

! resolve the C13 pockets 
      xa_function_species(1) = 'h1'  ! name of nuclide as defined in chem_def
        xa_function_weight(1) = 5
        xa_function_param(1) = 3d-4
      xa_function_species(2) = 'he4'  ! name of nuclide as defined in chem_def
        xa_function_weight(2) = 5
        xa_function_param(2) = 1d-2
      xa_function_species(3) = 'c13'  ! name of nuclide as defined in chem_def
        xa_function_weight(3) = 4
        xa_function_param(3) = 1d-3
      xa_function_species(4) = 'n14'  ! name of nuclide as defined in chem_def
        xa_function_weight(4) = 3
        xa_function_param(4) = 1d-2

! timesteps
! general
      varcontrol_target = 1d-4
        ! this is the target value for relative variation in the structure from one model to the next.
        ! the default timestep adjustment is to increase or reduce the timestep depending on whether
        ! the actual variation was smaller or greater than this value.

      delta_lgRho_cntr_limit = 0.05  
      delta_lgRho_cntr_hard_limit = 0.1
   
      delta_lgT_cntr_limit = 0.01  
      delta_lgT_cntr_hard_limit = 0.02

      delta_lgTeff_limit = 0.01  
      delta_lgTeff_hard_limit = 0.02
   
      delta_lgL_limit = 0.1
      delta_lgL_hard_limit = 0.2

! TP
      delta_lgL_He_limit = 0.01 ! for TP resolution
      !delta_lgL_He_hard_limit = 0.05
      lgL_He_burn_min = 2.0 ! ignore changes in lgL_He if value 
      		            ! is less than this

! DUP
      DUP_varcontrol_factor = 1.0
        ! multiply varcontrol_target by this when TP_state == 2, i.e., during dredge up during thermal pulse
      
! time step controls for interpulse

      dH_limit_min_H = 1.e-2
      dH_limit = 0.1

      dH_div_H_limit_min_H = 1d-5
      dH_div_H_limit = 0.5

! Overshooting
      mass_for_overshoot_full_on = 1.8 ! Msun units
      mass_for_overshoot_full_off = 1.1 ! Msun units

      overshoot_f_above_nonburn_core = 0.014
      overshoot_f_above_nonburn_shell = 0.014
      overshoot_f_below_nonburn_shell = 0.014
      overshoot_f_above_burn_h_core = 0.014
      overshoot_f_above_burn_h_shell = 0.014
      overshoot_f_below_burn_h_shell = 0.014
      overshoot_f_above_burn_he_core = 0.014
      overshoot_f_above_burn_he_shell = 0.014
      overshoot_f_below_burn_he_shell = 0.014
      overshoot_f_above_burn_z_core = 0.014
      overshoot_f_above_burn_z_shell = 0.014
      overshoot_f_below_burn_z_shell = 0.014

      overshoot_f0_above_nonburn_core = 0.004
      overshoot_f0_above_nonburn_shell = 0.004
      overshoot_f0_below_nonburn_shell = 0.004
      overshoot_f0_above_burn_h_core = 0.004
      overshoot_f0_above_burn_h_shell = 0.004
      overshoot_f0_below_burn_h_shell = 0.004
      overshoot_f0_above_burn_he_core = 0.004
      overshoot_f0_above_burn_he_shell = 0.004
      overshoot_f0_below_burn_he_shell = 0.004
      overshoot_f0_above_burn_z_core = 0.004
      overshoot_f0_above_burn_z_shell = 0.004
      overshoot_f0_below_burn_z_shell = 0.004

! C13 pocket      
      overshoot_below_noburn_shell_factor = 9
! He-shell flash convection zone
      ovr_below_burn_he_shell_factor = 0.72

! atmosphere option
      which_atm_option = 'photosphere_tables'
      which_atm_off_table_option = 'grey_and_kap'

! mass loss       
      RGB_wind_scheme = 'Reimers'
      AGB_wind_scheme = 'Blocker'
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_wind_eta = 0.5d0  
      Blocker_wind_eta = 0.05d0  

! opacities
      use_Type2_opacities = .true.
      Zbase = 0.01d0

/ ! end of controls namelist

