&star_job

      create_pre_main_sequence_model = .false.  

      ! change to agb isotope network
      change_net = .true.   
      new_net_name = 'agb.net'
      
      ! set initial values
      set_initial_age = .true.
      initial_age = 0 ! in years
      set_initial_model_number = .true.
      initial_model_number = 0
      
      ! save snapshot at run termination 
      !save_model_when_terminate = .true.
      !save_model_filename = '3.0M_stop_at_TAMS.mod'

      ! begin with a pre-built model
      load_saved_model = .true.
      load_model_filename = '3.0M_stop_at_TAMS.mod'

      show_net_reactions_info = .true.
      list_net_reactions = .true.

      ! sparse history / profile columns
      history_columns_file = './history_columns.list'
      profile_columns_file = './profile_columns.list'

      ! display on-screen plots
      pgstar_flag = .true.

      ! pause just before MESA exits
      pause_before_terminate = .true.

/ ! end of star_job namelist



&eos
  ! equation-of-state options
  ! see eos/defaults/eos.defaults

/ ! end of eos namelist


&kap
  ! opacity options
  ! see kap/defaults/kap.defaults
  kap_file_prefix = 'gs98'
  use_Type2_opacities = .true.
  Zbase = 0.02

/ ! end of kap namelist




&controls

       initial_mass = 3.0d0     ! in Msun units
       initial_z = 0.020    ! Solar Metallicity 

       ! select output frequency
       photo_interval = 1000
       write_header_frequency = 2
       terminal_interval = 5


       xa_central_lower_limit_species(1) = 'he4'
       xa_central_lower_limit(1) = 1d-6


! options for energy conservation (see MESA V, Section 3)
      energy_eqn_option = 'dedt'
      use_gold_tolerances = .false.

! winds

      cool_wind_RGB_scheme   = 'Reimers'
      cool_wind_AGB_scheme   = 'Blocker'
      RGB_to_AGB_wind_switch = 1d-4
      Reimers_scaling_factor = 0.1   
      Blocker_scaling_factor = 0.5   


! resolution

       mesh_delta_coeff = 2.0

! resolution


      !varcontrol_target = 1d-4

! MS Resolution 

      delta_lg_XH_cntr_limit = 0.01d0
      delta_lg_XH_cntr_max   = 0.0d0
      delta_lg_XH_cntr_min   = -6.0d0
      delta_lg_XH_cntr_hard_limit = 0.02d0

      delta_lg_XHe_cntr_limit = 0.02d0
      delta_lg_XHe_cntr_max   = 0.0d0
      delta_lg_XHe_cntr_min   = -6.0d0
      delta_lg_XHe_cntr_hard_limit = 0.04d0

      delta_lg_XC_cntr_limit = 0.01d0
      delta_lg_XC_cntr_max   = 0.0d0
      delta_lg_XC_cntr_min   = -6.0d0
      delta_lg_XC_cntr_hard_limit = 0.02d0

      ! this is mainly to resolve properly when the star goes off the main sequence
      delta_HR_limit = 0.01d0
      delta_lgL_He_limit = 0.1



         use_Ledoux_criterion = .true.

! semiconvection

      alpha_semiconvection = 1d-2
      semiconvection_option = 'Langer_85 mixing; gradT = gradr'


! thermohaline

      thermohaline_coeff = 1d0
      thermohaline_option = 'Kippenhahn'

/ ! end of controls namelist

