[Mesa-users] Low mass diffusion failure with remaining_source_mass
Evan Bauer
evan.bauer.astro at gmail.com
Fri Apr 7 21:06:50 UTC 2023
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> 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
> https://lists.mesastar.org/mailman/listinfo/mesa-users
>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230407/2fbc9a90/attachment.htm>
More information about the Mesa-users
mailing list