<html aria-label="message body"><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">I’ll take care of it.<div>M<br id="lineBreakAtBeginningOfMessage"><div><br><blockquote type="cite"><div>On 21. Jan 2026, at 19:49, Dimitris Menemenlis <dmenemenlis@gmail.com> wrote:</div><br class="Apple-interchange-newline"><div>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8"><div style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Martin, for time being I added your suggested modification to our local repo:<div><a href="https://github.com/MITgcm-contrib/llc_hires/blob/master/trillium/llc_4320/code/shelfice_thermodynamics.F">https://github.com/MITgcm-contrib/llc_hires/blob/master/trillium/llc_4320/code/shelfice_thermodynamics.F</a></div><div><br></div><div>But yes, I think this probably needs to make its way back to the main branch.</div><div>Will you take care of the PR? Or shall I try to figure out how to submit one?</div><div><br></div><div>D.</div><div><br id="lineBreakAtBeginningOfMessage"><div><br><blockquote type="cite"><div>On Jan 21, 2026, at 1:04 AM, Martin Losch <Martin.Losch@awi.de> wrote:</div><br class="Apple-interchange-newline"><div><div style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Do you think that this should be a PR?<div><br></div><div>Martin<br id="lineBreakAtBeginningOfMessage"><div><br><blockquote type="cite"><div>On 20. Jan 2026, at 18:38, Dimitris Menemenlis <dmenemenlis@gmail.com> wrote:</div><br class="Apple-interchange-newline"><div>
<div style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Martin, I just tested your modification below in my baby set-up and it works fine.<div><br></div><div>You just need to add a “-d 0” in "IF ( saltFreeze .EQ. 0. _d 0 ) THEN"<div><br></div><div>Thanks for quick response.</div><div><br></div><div>D.</div><div><br></div><div>P.S.: I still owe you a response for the “Supercooled waters” post. I don’t know how you and Oliver and Jean-Michel are so efficient!</div><div><br></div><div>P.P.S.: Happy New Year and belated Merry Christmas to you and family.</div><div><br></div><div><br><blockquote type="cite"><div>On Jan 20, 2026, at 5:37 AM, Martin Losch <Martin.Losch@awi.de> wrote:</div><br class="Apple-interchange-newline"><div><div style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Dimitris and Kayhan,<div><br></div><div>great catch. I suggest a different fix: Even though it’s unlikely, we may want to allow saltFreeze=0. So If sLoc = 0, then we have freshWaterFlux = rUnit2mass * shiTransCoeffS(i,j,bi,bj) * ( 1. _d 0 - 0. / 0. ) and could drop the part in the parentheses, i.e.</div><div><br></div><div><div> IF ( saltFreeze .EQ. 0. _d ) THEN</div><div> freshWaterFlux = rUnit2mass</div><div> & * shiTransCoeffS(i,j,bi,bj)</div><div> ELSE</div><div> freshWaterFlux = rUnit2mass</div><div> & * shiTransCoeffS(i,j,bi,bj)</div><div> & * ( 1. _d 0 - sLoc(I,J) / saltFreeze )</div><div> ENDIF</div></div><div><br></div><div>What do you think? Does that work for you?</div><div><br></div><div>Martin</div><div><div><div><br><blockquote type="cite"><div>On 20. Jan 2026, at 00:46, Dimitris Menemenlis <dimitris.menemenlis@sjsu.edu> wrote:</div><br class="Apple-interchange-newline"><div><div style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">shelfice_thermodynamics.F can crash when the line<br> saltFreeze = (- bqe - SQRT(discrim))*recip_aqe<br>returns 0, so that<br> freshWaterFlux = rUnit2mass<br> & * shiTransCoeffS(i,j,bi,bj)<br> & * ( 1. _d 0 - sLoc(I,J) / saltFreeze )<br>returns NaN.<br><br>This appears to be a consequence of "sLoc(I,J) = 0”<br>One possible fix is to add:<br> sLoc(I,J) = MAX(salt(I,J,K,bi,bj), 0.1 _d 0)<br><br>>>>>>>><br><br><div>Kayhan Momeni, from his ongoing llc4320 simulation has extracted all the model state variables</div><div>at the time step before (3205150) and during (3205151) the first NaN:</div><div><a href="https://nasa-ext.box.com/s/5ah3cffma49vkvfwgjsown7ytgs8i5i1">https://nasa-ext.box.com/s/5ah3cffma49vkvfwgjsown7ytgs8i5i1</a></div><div><br></div><div>I have created a small (17x17x173) regional set-up that (approximately) recreates this crash:</div><div><a href="https://nasa-ext.box.com/s/varcd6iu3ilg44bzly9hnzee1wnigfod">https://nasa-ext.box.com/s/varcd6iu3ilg44bzly9hnzee1wnigfod</a></div><div><a href="https://nasa-ext.box.com/s/bpq9p9yksksdhgr2knofj9fsqj0ntv5x">https://nasa-ext.box.com/s/bpq9p9yksksdhgr2knofj9fsqj0ntv5x</a></div><div><a href="https://nasa-ext.box.com/s/atgqgp43st6o4u4976gfr5not465zbdr">https://nasa-ext.box.com/s/atgqgp43st6o4u4976gfr5not465zbdr</a></div><div><br></div><div>I have verified that the small set-up is able to reproduce the SHICE_fwFlux NaN, the white grid cell at center of:</div><div><a href="https://nasa-ext.box.com/s/trhs4kd8vpt9tx4t8riq934olrks3h30">https://nasa-ext.box.com/s/trhs4kd8vpt9tx4t8riq934olrks3h30</a></div><div><br></div><div>All of the above are under: <a href="https://nasa-ext.box.com/s/q4sx5w4bnwui4v1afec9k2kb41ejml3v">https://nasa-ext.box.com/s/q4sx5w4bnwui4v1afec9k2kb41ejml3v</a></div><div><br></div></div>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></blockquote></div><br></div></div></div>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></blockquote></div><br></div></div>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></blockquote></div><br></div></div>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></blockquote></div><br></div></div>_______________________________________________<br>MITgcm-support mailing list<br>MITgcm-support@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br></div></blockquote></div><br></div></body></html>