[Mesa-users] crashes at Psurf=1 Ba

[Sho Shibata]

Dear Philip and all,

Thank you very much for your help. After applying the fixes, my code is now
running successfully. I really appreciate everyone’s kind support!

Best regards,
Sho

2025年11月24日(月) 8:50 Philip Mocz <pmocz at flatironinstitute.org>:

> Hi Sho,
>
> I'm suspecting we have a long-standing bug in the code when calculating
> residuals near Psurf = 1. I just found these lines in the code:
>
> https://github.com/MESAHub/mesa/pull/890/files
>
> We should have this fixed in time for the next official release.
>
> Thanks Sho for the excellent reporting and write-up of the problem, and Eb
> and Warrick for excellent diagnostics and suggestions that helped track
> this issue down quickly.
>
> Best,
> Philip
>
> On Sat, Nov 22, 2025 at 9:05 PM Sho Shibata via Mesa-users <
> mesa-users at lists.mesastar.org> wrote:
>
>> Dear EbF,
>>
>> Thank you for creating the GitHub thread. I really appreciate your kind
>> support.
>> Let’s continue the discussion there.
>>
>> Best regards,
>> Sho
>>
>>
>> 22.11.2025 16:41、Farag, Ebraheem <ebraheem.farag at yale.edu>のメール:
>>
>> Hello Sho,
>>
>> I can reproduce your issue, but I am unaware of a solution. As such, I've
>> made an issue on the MESA github for this thread.
>> https://github.com/MESAHub/mesa/issues/889
>> <https://urldefense.com/v3/__https://github.com/MESAHub/mesa/issues/889__;!!DSb-azq1wVFtOg!S6nlaV0ALsHRJkDoNeP_hiZM9rgbrxr7CSdsm8150-3Hs-2_rDkqxoOE1fJOjUIVX02ICyTF9QQaktetfLSz7iYRB6RO__ag$>
>>
>> Perhaps the discussion can continue there.
>>
>> -EbF
>> ------------------------------
>> *From:* Sho Shibata <s.shibata423 at gmail.com>
>> *Sent:* Saturday, November 22, 2025 9:08 AM
>> *To:* Farag, Ebraheem <ebraheem.farag at yale.edu>
>> *Cc:* Jeremy Sakstein via Mesa-users <mesa-users at lists.mesastar.org>
>> *Subject:* Re: [Mesa-users] crashes at Psurf=1 Ba
>>
>> Hi EbF,
>>
>> I am sharing the link to download my code. You can compile it by running:
>> ./mk, and run the simulation with: ./rn1.
>> change_atm.zip
>> <https://urldefense.com/v3/__https://drive.google.com/file/d/1S1QeYyWI2Psii4jyT1GNYpXVkoYCtR-m/view?usp=sharing__;!!DSb-azq1wVFtOg!S6nlaV0ALsHRJkDoNeP_hiZM9rgbrxr7CSdsm8150-3Hs-2_rDkqxoOE1fJOjUIVX02ICyTF9QQaktetfLSz7iYRB_krGwfO$>
>> drive.google.com
>> <https://urldefense.com/v3/__https://drive.google.com/file/d/1S1QeYyWI2Psii4jyT1GNYpXVkoYCtR-m/view?usp=sharing__;!!DSb-azq1wVFtOg!S6nlaV0ALsHRJkDoNeP_hiZM9rgbrxr7CSdsm8150-3Hs-2_rDkqxoOE1fJOjUIVX02ICyTF9QQaktetfLSz7iYRB_krGwfO$>
>> <drive_2020q4_32dp.png>
>> <https://urldefense.com/v3/__https://drive.google.com/file/d/1S1QeYyWI2Psii4jyT1GNYpXVkoYCtR-m/view?usp=sharing__;!!DSb-azq1wVFtOg!S6nlaV0ALsHRJkDoNeP_hiZM9rgbrxr7CSdsm8150-3Hs-2_rDkqxoOE1fJOjUIVX02ICyTF9QQaktetfLSz7iYRB_krGwfO$>
>> <https://urldefense.com/v3/__https://drive.google.com/file/d/1S1QeYyWI2Psii4jyT1GNYpXVkoYCtR-m/view?usp=sharing__;!!DSb-azq1wVFtOg!S6nlaV0ALsHRJkDoNeP_hiZM9rgbrxr7CSdsm8150-3Hs-2_rDkqxoOE1fJOjUIVX02ICyTF9QQaktetfLSz7iYRB_krGwfO$>
>>
>> The compilation should work in any standard MESA environment, but if you
>> encounter issues, please try copying run_star_extras.f90 into your working
>> directory.
>>
>> I would appreciate it if you could try running the code and let me know
>> if you notice anything or run into any problems.
>>
>> Sho Shibata
>> Department of Earth, Environmental and Planetary Sciences,
>> Rice University, Houston, TX, USA
>> Email: s.shibata423 at gmail.com
>>
>>
>>
>> 21.11.2025 23:43、Farag, Ebraheem <ebraheem.farag at yale.edu>のメール:
>>
>> Hello Sho,
>>
>> I am not aware of such an issue., but from what you've shared, it could
>> be a bug, but it could also be related to your model setup. Perhaps you can
>> share a working clean directory with us and, maybe, we can find the source.
>>
>> -EbF
>> ------------------------------
>> *From:* Sho Shibata <s.shibata423 at gmail.com>
>> *Sent:* Friday, November 21, 2025 11:30 PM
>> *To:* Farag, Ebraheem <ebraheem.farag at yale.edu>
>> *Cc:* Jeremy Sakstein via Mesa-users <mesa-users at lists.mesastar.org>
>> *Subject:* Re: [Mesa-users] crashes at Psurf=1 Ba
>>
>> Dear EbF,
>>
>> I see your point. The case with Tsurf=100 K was not good to show my
>> problem.
>>
>> I attached the log file of the first 1000 lines with Tsurf=300 K. In this
>> case, the EOS should be in the domain. I don’t set up pgstar, so I cannot
>> use "show_TRho_Profile_eos_regions = .true. “ for now. I will try it later.
>>
>> In the new log file, MESA got the EOS error at line 180, but after
>> several retries, MESA can proceed further. So this is not a big problem.
>>
>> From line 698, MESA retries again, and this is not related to the EOS
>> error. Most problematic crashes do not relate to EOS errors. So I think EOS
>> is not the main cause.
>>
>> Do you know any modules that could relate to 1 Ba in MESA?
>>
>>
>> Sho Shibata
>> Department of Earth, Environmental and Planetary Sciences,
>> Rice University, Houston, TX, USA
>> Email: s.shibata423 at gmail.com
>>
>>
>>
>> 21.11.2025 21:50、Farag, Ebraheem <ebraheem.farag at yale.edu>のメール:
>>
>> show_TRho_Profile_eos_regions = .true.
>>
>>
>> _______________________________________________
>> mesa-users at lists.mesastar.org
>>
>> https://urldefense.com/v3/__https://lists.mesastar.org/mailman/listinfo/mesa-users__;!!DSb-azq1wVFtOg!S6nlaV0ALsHRJkDoNeP_hiZM9rgbrxr7CSdsm8150-3Hs-2_rDkqxoOE1fJOjUIVX02ICyTF9QQaktetfLSz7iYRB8sC_R40$
>>
>>
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <https://lists.mesastar.org/pipermail/mesa-users/attachments/20251124/b9e46f97/attachment.htm>