[Mesa-users] Low mass diffusion failure with remaining_source_mass

[Saida Milena Diaz Castillo]

Dear,

Thank you for your email. I’m out of the office for the holidays from 3rd of April until 12th of April. I’ll reply to your email as soon as I can upon my return.

Cheers,

On 7. Apr 2023, at 19:07, 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-users at lists.mesastar.org
> https://lists.mesastar.org/mailman/listinfo/mesa-users
> 
-------------- next part --------------
A non-text attachment was scrubbed...
Name: not available
Type: text/enriched
Size: 1972 bytes
Desc: not available
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20230408/1606fe85/attachment.bin>