[Mesa-users] rel_run_E_err Warning

[Francis Timmes]

hi erica,

of course one should be aware of the potential issues when turning off warnings.

for the second item, in your inlist_astero_search_controls, there is

       ! terminate if simplex gets "small" enough
       simplex_x_atol = 1d-10 ! tolerance for absolute differences
       simplex_x_rtol = 1d-10 ! tolerance for relative differences

these might be too small relative to when you seek to terminate.
perhaps experiment with larger values, at least until one of the astero experts comment.

fxt






> On May 2, 2020, at 8:34 AM, Francis Timmes <fxt44 at mac.com> wrote:
> 
> hi erica,
> 
> for the first item, one can consider using
> 
> ! silence energy warnings
>     warn_when_large_rel_run_E_err = 1.0d99
>     warn_when_stop_checking_residuals = .false.
> 
> fxt
> 
> 
> 
> 
>> On May 1, 2020, at 12:45 PM, Erica Sawczynec <ericasaw at hawaii.edu> wrote:
>> 
>> Hi all!
>> 
>> I am working with the astero module using MESA 12115 and MESASDK 20190830. I am trying to model a star with individual pulsation frequencies, but I am running into a few problems:
>> 1) Throughout each iteration (from pre-main sequence to within the error box) I noticed an energy conservation warning: "WARNING: rel_run_E_err". The warning continues even when fitting the seiesmic frequencies.   When looking through some of the documentation for MESA I noticed a few things that might help this error: a. using gold tolerances and b. using the dedt form of the energy eq. After turning both of these parameters on, the warning still appears. Further, to check it was not a problem with my particular inlist, I also ran the MESA test suit for example-astero and found the same warning all the way up the main sequence. 
>> 
>> 2) I am also having difficultly getting the astero module to stop modeling once the parameters have converged. I am unsure what parameters to change that will end the entire run. I have tried: changing the chi2 search limits, changing the number of simplex iterations, changing the sigmas_coeff values for all the parameters, and changing the chi2 relative increase limit. Right now I just plot all the output values from the simplex results file and wait until I see the values have converged before I end the run (see the image convergence.png).
>> 
>> Here is some information about my os:
>> NAME="CentOS Linux"
>> VERSION="7 (Core)"
>> ID="centos"
>> ID_LIKE="rhel fedora"
>> VERSION_ID="7"
>> PRETTY_NAME="CentOS Linux 7 (Core)"
>> 
>> Below I have attached a text file that shows the terminal outputs for my star on the main sequence to right before the astero module begins a new model with a different combination of starting parameters.
>> I have also attached my inlist_astero and my inlist_astero_search_controls so you can see all the parameters that are set.
>> 
>> Thanks so much!
>> Erica
>> 
>> <covergence.png><inlist_astero.txt><inlist_astero_search_controls.txt><terminal_output.txt>_______________________________________________
>> mesa-users at lists.mesastar.org
>> https://lists.mesastar.org/mailman/listinfo/mesa-users
>> 
>