Page 1 of 1

mass non-monotonicity

Posted: Sat Sep 17, 2016 4:02 pm
by ehsan
Dear Rich,

I am computing a large grid of models, where a very small fraction, yet sizeable in number, of models fail because of a mass non-monotonicity; thus, the following error:

Code: Select all

ASSERT 'ALL(m(2:) >= m(:n-1))' failed at line 514 <gyre_evol_model:evol_model_t_full_>:
Because I cannot go back to MESA to recompute these specific snapshots along the evolution (where the multi-D grid is very large), I was wondering if there is a possibility through inlist to handle such cases?

Kind regards,
Ehsan.

Re: mass non-monotonicity

Posted: Sat Sep 17, 2016 4:06 pm
by rhtownsend
Hi Ehsan --

In such cases, you might just want to modify the code so that the mass monotonicity check is commented out. I don't think this will affect execution of the code, as interpolations use radius as the abscissa, not mass.

cheers,

Rich

Re: mass non-monotonicity

Posted: Sat Sep 17, 2016 4:44 pm
by ehsan
Thanks Rich,

It totally works ;-)

Regards,
Ehsan.

Re: mass non-monotonicity

Posted: Tue Mar 20, 2018 2:05 am
by raymond
Got it! Thanks.