[Mesa-users] Stop condition at the end of the RGB tip

[Adrian Ayala]

Hi!

I am trying to evolve a 1 solar mass model up to the start of the HB 
(horizontal branch phase) after the RGB tip. I tried the stop condition 
related to the central He mass abundance:

! stop passed horizontal branch
     HB_limit = 0.95d0


  I found very dramatic differences between mesa12 and 10 versions. In 
the second, the simulation lasted almost twice.

Should I use another stop criterion in MESA 10 to end just before HB? 
Please, find the inlist_project file in attach.

Best wishes,

Adrián.

-------------- next part --------------
! inlist to evolve a 0.82 solar mass star, Reimers 0.4 and Y 0.248

! For the sake of future readers of this file (yourself included),
! ONLY include the controls you are actually using.  DO NOT include
! all of the other controls that simply have their default values.

&star_job

  ! begin with a pre-main sequence model
    create_pre_main_sequence_model = .true.

  ! save a model at the end of the run
    save_model_when_terminate = .false.
    save_model_filename = 'M1.22solar_Reimers04_Y0273_HP_Chi0'

  ! display on-screen plots
    pgstar_flag = .true.

  ! Set of Castelli 2003 bolometric corrections in UBV (Johnson), RI (Cousins) and JHKL (Bessell)
  ! Can be used at the same time as the lcb98cor.dat file
  ! Filter names BC_v U-B B-V V-I V-R V-J V-K V-L H-K
    color_num_files=2
    color_file_names(2)='vmag_tables_Adrian.dat'
    color_num_colors(2)=9

/ !end of star_job namelist


&controls

  ! starting specifications
    initial_mass = 1.0d0
    initial_y = 0.2676d0
    initial_z = 0.0125d0

  ! use_Type2_opacities = .true.
  ! Zbase = 0.001

  ! Use neutrino hook
    use_other_neu = .true.

  ! Activate semi convection while mixing
    use_Ledoux_criterion = .true.
    alpha_semiconvection = 1    
    semiconvection_option = 'Langer_85'
    thermohaline_coeff = 1
    thermohaline_option = 'Brown_Garaud_Stellmach_13'

  ! mass gain or loss
    mdot_omega_power = 0.43
    cool_wind_RGB_scheme = 'Reimers'
    Reimers_scaling_factor = 0.2d0
    cool_wind_AGB_scheme = 'Blocker'
    Blocker_scaling_factor = 0.2d0  
    RGB_to_AGB_wind_switch = 1d-4
  !  cool_wind_full_on_T = 8d3
  !  cool_wind_full_off_T = 1d4

  ! relax_mass = .true. ! gradually change total mass by a wind ! Comentado porque la actual version no reconoce el parametro
  ! relax_initial_mass = .true. ! Comentado porque la actual version no reconoce el parametro

    
  ! MLT
    mixing_length_alpha = 1.9658
    MLT_option = 'Henyey'

  ! atmosphere
    atm_table = 'photosphere'

  ! predictive mixing (to avoid breathing pulses and allow core expansion)
    predictive_mix(1) = .true.
    predictive_zone_type(1) = 'any'
    predictive_zone_loc(1) = 'any'
    predictive_bdy_loc(1) = 'any'
    predictive_superad_thresh(1) = 0.01
    predictive_avoid_reversal(1) = 'he4'

  ! stop passed horizontal branch 
    HB_limit = 0.95d0

    write_pulse_data_with_profile = .true.
    pulse_data_format = 'GYRE'


      
/ ! end of controls namelist