[Mesa-users] Low mass diffusion failure with remaining_source_mass

[Roulston, Benjamin]

Hi Evan,

Sorry for the slow response, your help is greatly appreciated! Here is the starting model for this run.



Cheers,
- Ben

On Apr 7, 2023, at 2:06 PM, Evan Bauer <evan.bauer.astro at gmail.com> wrote:

Hi Ben,

I'd be very interested in trying to help debug this. I don't think I've ever run across this particular failure mode of diffusion before, and I thought I had collected most of them by this point! Can you attach the starting model? Right now trying to run your inlists to reproduce the failure gives me

failed to open MODELS/1Gyr/1Gyr_010M_Z_000018_w001.mod

Cheers,
Evan

On Apr 7, 2023, at 1:07 PM, Roulston, Benjamin via Mesa-users <mesa-users at lists.mesastar.org<mailto:mesa-users at lists.mesastar.org>> wrote:

Hello,

I’m having an issue trying to evolve some lower mass stellar models from 1Gyr to TAMS (or 14Gyr). They seem to proceed normally, but then abruptly fail with the following message

remaining_source_mass           1         136         136                       NaN    4.2369868473723005D+25    4.2369868473723005D+25

I think this has something to do with diffusion, but I’m not are exactly what issue is happening to cause the run to fail. I looked in star/private/diffusion_procs.f90 around line #749, and this is seeming to say that the remaining_source_mass value is NaN. I’m still looking deeper to see if I can figure out where this is coming from, but if anyone who knows more has any help that would be great!

For reference, I’m doing a grid in mass from 0.1Msun to 0.9 Msun. Using the same exact inlist except for the mass, all of the models run EXCEPT the 0.1Msun, the lowest mass case.




I am using MESA 22.10.1 and also using the SDK (mesasdk-aarch64-macos-22.10.1). MESA builds and runs fine. I include the output from MESA as an attachment. As well as the inlists used.

Cheers,
- Ben
<MESA_output.txt><inlist><inlist_1Gyr_TAMS><inlist_pgstar><history_columns.list><profile_columns.list><run_star_extras.f90>_______________________________________________
mesa-users at lists.mesastar.org<mailto:mesa-users at lists.mesastar.org>
https://lists.mesastar.org/mailman/listinfo/mesa-users



-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230415/8b455f94/attachment.htm>
-------------- next part --------------
A non-text attachment was scrubbed...
Name: 1Gyr_010M_Z_000018_w001.mod
Type: audio/mod
Size: 648902 bytes
Desc: 1Gyr_010M_Z_000018_w001.mod
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230415/8b455f94/attachment.bin>