Hi All,
I've posted this same question to the MESA listserv, but perhaps this is a better place:
I'm stopping and starting MESA on very short time intervals (on the order of half a year, in some cases) and generating GYRE profiles at each of those times.
Proceeding chronologically, for the first 5 or so such .gyre models out of a set of, say, 20, the stellar radius at the appropriate value of t in MESA's history file, say R(t), is in near exact agreement with the radius passed to the .gyre file.
With increasing evolutionary time, however, an offset between R(t) in the history file and the radius given to .gyre (say R_gyre) emerges. The discrepancy between R(t) and R_gyre worsens the more .gyre models I try to produce from a single run of ./star.
All of the .gyre files in a set are generated from a single evolutionary run---though that run may proceed discontinuously---using the same inlist and controls. I've copied one of my inlists below for example.
My instinct is that this is due to some quirk or error which iteratively compounds a loss of precision or numerical misalignment each time I stop and restart the models.
Has anyone encountered this before or know a way around it?
Thank you
----------------------------------------------------------------------------------------------------------------------
#example inlist
----------------------------------------------------------------------------------------------------------------------
&star_job
mesa_dir = '/avatar/mjoyce/mesa-r10398/'
history_columns_file='/avatar/mjoyce/mesa-r10398/star/work/work_columns.list'
load_saved_model = .true.
saved_model_name = 'M1.7_A2.0_Z0.011_1.64749870014Gyr.inlist.mod '
create_pre_main_sequence_model = .false.
save_model_when_terminate = .true.
save_model_filename= 'M1.7_A2.0_Z0.011_1.64749874806Gyr.inlist.mod'
write_profile_when_terminate = .true.
filename_for_profile_when_terminate = 'profile_M1.7_A2.0_Z0.011_1.64749874806Gyr.inlist.data'
profile_columns_file = 'profile_columns.list'
save_pulse_data_when_terminate = .true.
!save_pulse_data_filename = 'M1.7_A2.0_Z0.011_1.64749874806Gyr.inlist.gong'
save_pulse_data_filename = 'M1.7_A2.0_Z0.011_1.64749874806Gyr.inlist.gyre'
pgstar_flag = .false.
save_pgstar_files_when_terminate = .false.
change_net = .true. ! switch nuclear reaction network
new_net_name = 'cno_extras_o18_to_mg26.net'
/ !end of star_job namelist
&controls
star_history_name = 'history_M1.7_A2.0_Z0.011_1.64749874806Gyr.inlist.data'
!--------- new pulsation output parameters 12/21/18 ------------
add_atmosphere_to_pulse_data = .true. !! maybe an issue with pulses not going low enough?
pulse_data_format = 'GYRE'
!max_number_backups = 0
!max_number_retries = 5 !! models keep terminating with this on
!max_model_number = 500
initial_mass =1.7d0
initial_y = 0.27d0 !solar
initial_z = 0.011 !solar
mixing_length_alpha = 2.0
stop_at_TP=.false.
!envelope_mass_limit = 0.005d0 !! NEW 1/4/19, to auto-stop after TP #fxt ! changed from 0.01d0
!! turning off bc not compatible with GYRE !! which_atm_option = 'photosphere_tables' !! NEW 1/4/19
terminal_show_age_in_years = .true.
max_age = 1.64749874806d9 ! age of Sun in Gyr
cool_wind_AGB_scheme = 'Blocker' !! recommended; also given in Fadeyev paper
Blocker_scaling_factor = 0.5
RGB_to_AGB_wind_switch = 1d-4
mesh_dlog_pp_dlogP_extra = 0.25
mesh_dlog_cno_dlogP_extra = 0.25
mesh_dlog_3alf_dlogP_extra = 0.225
mesh_dlog_burn_c_dlogP_extra = 0.225
mesh_dlog_burn_n_dlogP_extra = 0.225
mesh_dlog_burn_o_dlogP_extra = 0.225
mesh_logX_species(1) = 'h1'
mesh_logX_min_for_extra(1) = -6
mesh_dlogX_dlogP_extra(1) = 0.25
mesh_logX_species(1) = 'he4'
mesh_logX_min_for_extra(1) = -6
mesh_dlogX_dlogP_extra(1) = 0.25
xtra_coef_os_above_nonburn = 0.5
xtra_coef_os_below_nonburn = 0.25
xtra_coef_os_above_burn_h = 0.25
xtra_coef_os_below_burn_h = 0.15
xtra_coef_os_above_burn_he = 0.25
xtra_coef_os_below_burn_he = 0.25
xtra_coef_os_above_burn_z = 0.25
xtra_coef_os_below_burn_z = 0.25
xtra_dist_os_above_nonburn = 0.1
xtra_dist_os_below_nonburn = 0.2
xtra_dist_os_above_burn_h = 0.2
xtra_dist_os_below_burn_h = 0.2
xtra_dist_os_above_burn_he = 0.2
xtra_dist_os_below_burn_he = 0.2
xtra_dist_os_above_burn_z = 0.2
xtra_dist_os_below_burn_z = 0.2
history_interval = 1
/ ! end of controls namelist