Page 1 of 2
Double meshpoints causing error?
Posted: Wed Mar 01, 2017 10:05 am
by Mathias
Dear Rich
Using GYRE5.0 I get the following error message:
ASSERT 'ALL(x(2:) > x(:SIZE(x)-1))' failed at line 20 <gyre_r_interp:r_interp_t_>:
Non-monotonic abscissa
I guess this is because of the double meshpoints in the pulse info files from MESA.
I'm using MESA8845, with "add_double_points_to_pulse_info = .true".
What do you recommend to solve this issue?
Best regards
Mathias
Re: Double meshpoints causing error?
Posted: Wed Mar 01, 2017 12:48 pm
by rhtownsend
Mathias wrote:Dear Rich
Using GYRE5.0 I get the following error message:
ASSERT 'ALL(x(2:) > x(:SIZE(x)-1))' failed at line 20 <gyre_r_interp:r_interp_t_>:
Non-monotonic abscissa
I guess this is because of the double meshpoints in the pulse info files from MESA.
I'm using MESA8845, with "add_double_points_to_pulse_info = .true".
What do you recommend to solve this issue?
Best regards
Mathias
Hi Mathias --
This shouldn't happen with GYRE 5.0 -- one of the big changes since 4.x is the ability to handle double mesh points. Can you post an inlist and a model file?
cheers,
Rich
Re: Double meshpoints causing error?
Posted: Thu Mar 02, 2017 11:01 am
by Mathias
Dear Rich
The model file (4.6MB) is apparently too large to add as an attachement, even if I split the file and compress it (resulting in a file of about 0.6MB), I still get the error message "File too large".
Best regards
Mathias
Re: Double meshpoints causing error?
Posted: Thu Mar 02, 2017 4:26 pm
by rhtownsend
Mathias wrote:Dear Rich
The model file (4.6MB) is apparently too large to add as an attachement, even if I split the file and compress it (resulting in a file of about 0.6MB), I still get the error message "File too large".
Best regards
Mathias
Hi Mathias --
I've increased the file upload limit. Can you try submitting again (in compressed form)?
cheers,
Rich
Re: Double meshpoints causing error?
Posted: Fri Mar 03, 2017 8:17 am
by Mathias
Hi Rich
Thank you, the inlist and mode file are attached.
Best regards
Mathias
Re: Double meshpoints causing error?
Posted: Fri Mar 03, 2017 12:12 pm
by rhtownsend
Mathias wrote:Dear Rich
The model file (4.6MB) is apparently too large to add as an attachement, even if I split the file and compress it (resulting in a file of about 0.6MB), I still get the error message "File too large".
Best regards
Mathias
Hi Mathias --
I've tracked down the problem -- in your model, point number 9223 (counting from 1) has a smaller radius than the previous point.
I'm guessing that this point is where the atmosphere joins on to the interior. Can I ask what MESA version you're using, and could you send an inlist for the MESA model? This is a bug in MESA, and I'd like to get it fixed.
cheers,
Rich
Re: Double meshpoints causing error?
Posted: Mon Mar 06, 2017 5:05 am
by Mathias
Dear Rich
The inlists are attached, I'm using MESA version 8845.
Best regards
Mathias
Re: Double meshpoints causing error?
Posted: Mon Mar 06, 2017 9:32 am
by rhtownsend
Mathias wrote:Dear Rich
The inlists are attached, I'm using MESA version 8845.
Best regards
Mathias
Hi Mathias --
One more thing -- can you post a MESA profile for the problem model?
cheers,
RIch
Re: Double meshpoints causing error?
Posted: Tue Mar 07, 2017 12:32 pm
by Mathias
Hi Rich
This is the MESA profile corresponding to the problematic model file.
Best regards
Mathias
Re: Double meshpoints causing error?
Posted: Tue Mar 07, 2017 1:33 pm
by rhtownsend
Mathias wrote:Hi Rich
This is the MESA profile corresponding to the problematic model file.
Best regards
Mathias
Hi Mathias --
Comparing the profile and the GYRE model, it seems the GYRE model has an extra 50 points in the surface layers, which presumably represent the stellar atmosphere. However, looking through your inlists I can't see any flags which tell MESA to add an atmosphere onto the GYRE model.
Is it possible that the inlists you sent me aren't exactly the same ones that you use to calculate the profile & GYRE model?
cheers,
Rich