! inlist_agb2M: this is the Falk approved inlist for a 2Msun Z=0.01
!  stellar evolution track to the end of the AGB (25JUN2010). The
!  track from this inlist should get to the end in about 65000 cycles
!  and under 24hrs runtime (OpenMP_THREADS=2) to finish. The track has
!  been checked against the latest (last?) 2Msun, Z=0.01 stellar
!  evolution track E0146D4 from Falk's EVOL stellar evolution
!  code. That track is very similar to the ones published by Herwig &
!  Austin (2004ApJ...613L..73). In particular this track has been
!  checked for:
!
!   * MS lifetime (dt_MS<3%)
!   * occurence (dM_H1<3%) and evolution of TP
!   * peak-He-shell flash luminosity
!   * intershell abundances, dependence on overshooting
!   * occurence and evolution of the 3rd DUP
!   * convergence of 3rd DUP
!   * C star formation
!   * formation and evolution of C13 pocket
!   * evolution of core mass
!   * mass loss
!   * interpulse duration (<6%)
!
! With respect to all these criteria the MESA track is agreeing with
!  the EVOL code results. However, MESA tracks are run in 1/20th of
!  the time needed for the same run done with EVOL. MESA runs do not
!  need any intervention while EVOL runs need several
!  'interventions'. Finally, MESA results are much smoother
!  demonstrating the better MESA performance and overall convergence.
!
! I therfore declare MESA to be verified for low-mass AGB stellar
!  evolution.

&star_job
      
      mesa_data_dir = '/home/umberto/mesa_July_2012/data' ! for eos and opacity tables
      create_pre_main_sequence_model = .true.

      save_model_number = 100
      save_model_filename = '100.mod'

      !load_saved_model = .true.
      !saved_model_name =  'AGBpreDUP2M.mod'
      !saved_model_name =  '2MpreMS'

      change_net = .true.      
      new_net_name = 'agb.net'

      set_uniform_xa_from_file = .true.
      file_for_uniform_xa = 'xa_2.0E-02GN93.mesa_agb'

      set_rates_preference = .true. ! for use by net + rates modules
      new_rates_preference = 1 
         ! 1 = NACRE rates -- this is the default
         ! 2 = jina reaclib rates -- to match jina where possible
         
      set_rate_c12ag = 'Kunz' ! empty string means ignore this control
         ! one of 'NACRE', 'Buchmann', 'Kunz', or 'CF88'
         ! note: original CF88 rate is actually multiplied by 1.7 as in Timmes' rates
      set_rate_n14pg = 'Imbriani' ! empty string means ignore this control
         ! one of 'NACRE', 'Imbriani', or 'CF88'
      set_rate_3a = 'Fynbo' ! empty string means ignore this control
         ! one of 'NACRE', 'Fynbo', or 'CF88'
	 
      use_se_output = .true.
      change_v_flag = .true.
      new_v_flag = .true.
      !relax_tau_factor = .true.	     
      !relax_to_this_tau_factor = 10.


      !log_columns_file = 'log_columns.list' 

/ ! end of star_job namelist

&controls

         ! artificial viscosity -- only applies when using velocity variables
         use_artificial_viscosity = .true.
         l1_coef = 1.e3 ! increase to as much as 1d3 to suppress pulsations

      initial_mass = 2.00
      initial_z = 0.020
      mixing_length_alpha = 1.73

      operator_coupling_choice = 0
          ! 0 -- fully coupled

! some logistics
      max_model_number = 320000 ! negative means no maximum
      !max_model_number = 95000 ! negative means no maximum
      !max_model_number = 10000 ! negative means no maximum

      photostep = 500
      profile_interval = 100
      max_num_profile_models = 100 ! maximum number of saved profiles
      !star_mass_min_limit = 1.17
      star_mass_min_limit = 0.69
      log_cnt = 1
      terminal_cnt = 10
      write_header_frequency = 10

! mesh
      max_allowed_nz = 8000 ! 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) = 5
        xa_function_param(3) = 1d-4
      xa_function_species(4) = 'n14'  ! name of nuclide as defined in chem_def
        xa_function_weight(4) = 4
        xa_function_param(4) = 1d-6

! 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

      !radiation_turbulence_coeff = 1 ! 1 is reasonable value for this coefficient

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

      overshoot_f_above_nonburn = 0.014
      overshoot_f_below_nonburn = 0.014
      overshoot_f_above_burn_h = 0.014
      overshoot_f_below_burn_h = 0.014
      overshoot_f_above_burn_he = 0.024
      overshoot_f_below_burn_he = 0.024
      overshoot_f_above_burn_z = 0.014
      overshoot_f_below_burn_z = 0.014

      !overshoot_D2_above_nonburn = 100000000000
      overshoot_D2_below_nonburn = 100000000000
      !overshoot_D2_above_burn_h = 100000000000
      overshoot_D2_below_burn_h = 100000000000
      !overshoot_D2_above_burn_he = 100000
      overshoot_D2_below_burn_he = 100000
      !overshoot_D2_above_burn_z = -0.0001
      overshoot_D2_below_burn_z = -0.0001
         
      !overshoot_f2_above_nonburn = 0.24
      overshoot_f2_below_nonburn = 0.24
      !overshoot_f2_above_burn_h = 0.24
      overshoot_f2_below_burn_h = 0.24
      !overshoot_f2_above_burn_he = 0.14
      overshoot_f2_below_burn_he = 0.14
      !overshoot_f2_above_burn_z = -0.00001
      overshoot_f2_below_burn_z = -0.00001   
! C13 pocket      
      overshoot_below_noburn_factor = 1
! He-shell flash convection zone
      !ovr_below_burn_he_factor = 0.5714

! atmosphere option
      which_atm_option = 'simple_photosphere'
       ! 'simple_photosphere'	
          ! don't integrate, just estimate for tau=2/3

! 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.04d0  
      !Blocker_wind_eta = 0.01d0  

! opacities
      use_CO_enhanced_opacities = .true.
      base_Z = 0.02
      base_fC = 0.17123643187d0 ! base value of XC/Z
      base_fO = 0.4812175d0 ! base value of XO/Z

      max_logT_for_special_kap = 5.7
      special_kap_use_CO_enhanced = .false.
      special_kap_Z_base = 0.0
      !opacity_max = 110.
      opacity_factor = 0.95

/ ! end of controls namelist