[MITgcm-support] shelfice_thermodynamics.F crash

[Martin Losch]

I don’t know. If both solutions give the same results, I’d go with SHELFICE_USE_HEATBALANCE_FOR_FRESHWATERFLUX defined, as this is in compact form what is implemented in pkg/steep_icecavity. I have no experience with this version. The other option is a hack that just avoids the division by zero, this is only physically correct if sLoc is really zero at the same time.

M.

> On 22. Jan 2026, at 21:05, Dimitris Menemenlis <dimitris.menemenlis at sjsu.edu> wrote:
> 
> Done.  I left a review.  Both options work.
> Which one do you recommend that Kayhan and I use?
> with or without SHELFICE_USE_HEATBALANCE_FOR_FRESHWATERFLUX defined?
> 
> Thanks a million for super-quick turnaround on this request.
> 
> D.
> 
> 
>> On Jan 22, 2026, at 4:58 AM, Martin Losch <Martin.Losch at awi.de> wrote:
>> 
>> Please have a look at PR #968
>> https://github.com/MITgcm/MITgcm/pull/968
>> 
>> @Dimitry, it would be great, if you could try this code in your little test, and also review it in github.
>> 
>> M
> 
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20260123/648a96df/attachment.html>