[Mesa-users] rel_run_E_err Warning

[Jennifer Van Saders]

Dear Frank,

When should we be concerned that the rel_run_E_err is thrown? Is there a
community best practice for solar models as to what this threshold should
be? And, if we're not happy with where it is now, is the best path
forward increasing the mess points and decreasing timesteps?

Thanks!

Jen

On Sat, May 2, 2020 at 7:33 AM Francis Timmes via Mesa-users <
mesa-users at lists.mesastar.org> wrote:

> 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://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!KGKeukY!ipaYn4IY1MG7yu0_mxlNicrfkO3gmZSILAJFUAwWSseP034Uq7x91fd7YDgqhGttrCMgcq1FUDvpYw$
> >>
> >
>
> _______________________________________________
> mesa-users at lists.mesastar.org
>
> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!KGKeukY!ipaYn4IY1MG7yu0_mxlNicrfkO3gmZSILAJFUAwWSseP034Uq7x91fd7YDgqhGttrCMgcq1FUDvpYw$
>
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20200504/ccaa9ec7/attachment.htm>