&star_job
  ! see star/defaults/star_job.defaults

  ! begin with a pre-main sequence model
    create_pre_main_sequence_model = .false.
    load_saved_model = .true.
    load_model_filename = 'stop_at_RG.mod'

    relax_initial_mass_to_remove_H_env = .true.
    extra_mass_retained_by_remove_H_env = 0.4
    lg_max_abs_mdot = 0

  ! save a model at the end of the run
    save_model_when_terminate = .true.
    save_model_filename = 'stripped_at_TPAGB.mod'

  ! display on-screen plots
    pgstar_flag = .false.

/ ! end of star_job namelist


&eos
  ! eos options
  ! see eos/defaults/eos.defaults

/ ! end of eos namelist


&kap
  ! kap options
  ! see kap/defaults/kap.defaults
/ ! end of kap namelist


&controls
  ! see star/defaults/controls.defaults

  ! when to stop
    stop_at_phase_TP_AGB = .true.
    stop_at_phase_WDCS = .true.
  ! wind
    cool_wind_rgb_scheme = 'Reimers'
    cool_wind_agb_scheme = 'Blocker'
    Reimers_scaling_factor = 0.477
    Blocker_scaling_factor = 0.1
  ! atmosphere

  ! rotation

  ! element diffusion

  ! mlt
    okay_to_reduce_gradT_excess = .true.

  ! mixing
  !  overshoot_scheme(1) = 'exponential'
  !  overshoot_zone_type(1) = 'any'
  !  overshoot_zone_loc(1) = 'any'
  !  overshoot_bdy_loc(1) = 'any'
  !  overshoot_f(1) = 0.014
  !  overshoot_f0(1) = 0.004

  ! timesteps
    min_timestep_limit = 1d-20
  !  varcontrol_target = 1d-4

  ! mesh
  !  delta_lgTeff_limit = 0.01d0
  !  delta_lgTeff_hard_limit = 0.05d0
  !  relax_hard_limits_after_retry = .false.

  !  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
  !  mesh_delta_coeff = 0.5


  ! timesteps

  ! mesh

  ! solver
     ! options for energy conservation (see MESA V, Section 3)
     energy_eqn_option = 'eps_grav'
     convergence_ignore_equL_residuals = .true.

  ! output

/ ! end of controls namelist
