Double meshpoints causing error?

[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

[rhtownsend]

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

[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

[rhtownsend]

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

[Mathias]

Hi Rich

Thank you, the inlist and mode file are attached.

Best regards
Mathias

[rhtownsend]

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

[Mathias]

Dear Rich

The inlists are attached, I'm using MESA version 8845.

Best regards
Mathias

[rhtownsend]

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

[Mathias]

Hi Rich

This is the MESA profile corresponding to the problematic model file.

Best regards
Mathias

[rhtownsend]

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