<html><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;">In case you are interested: The problem is related to the (admittedly complicated) computation of the GGL90mixingLength. The simpler versions mxlMaxFlag=0,1 work. That should be possible to fix!<div><br></div><div>Martin<br><div><br><blockquote type="cite"><div>On 30. Jun 2023, at 08:44, Martin Losch <Martin.Losch@awi.de> 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;">OK, thanks. So you suspect it could have to do with the domain decomposition for KPP?<div><br></div><div>If you have any insights/pointers as to where I should start with GGL90, I would greatly appreciate it (but I am not sure when I will get to this).</div><div><br></div><div>Martin<br><div><br><blockquote type="cite"><div>On 29. Jun 2023, at 19:41, Patrick Heimbach <heimbach@mit.edu> wrote:</div><br class="Apple-interchange-newline"><div><div class="WordSection1" style="page: WordSection1; caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><div style="margin: 0in; font-size: 11pt; font-family: Calibri, sans-serif;"><span style="font-size: 12pt;">KPP is a bit of a mystery.<o:p></o:p></span></div><div style="margin: 0in; font-size: 11pt; font-family: Calibri, sans-serif;"><span style="font-size: 12pt;">We were able to use it (Gamma hat code off in adjoint) in the old vector code, but never in the “UV” (domain-decomposed) code. I did not figure out, why. I suspect a bug in the adjoint (but cannot prove it and am thus probably wrong).<o:p></o:p></span></div><div style="margin: 0in; font-size: 11pt; font-family: Calibri, sans-serif;"><span style="font-size: 12pt;">p.<o:p></o:p></span></div><div style="margin: 0in; font-size: 11pt; font-family: Calibri, sans-serif;"><span style="font-size: 12pt;"><o:p> </o:p></span></div><div id="mail-editor-reference-message-container"><div><div style="border-width: 1pt medium medium; border-style: solid none none; border-color: rgb(181, 196, 223) currentcolor currentcolor; border-image: none; padding: 3pt 0in 0in;"><p class="MsoNormal" style="margin: 0in 0in 12pt 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><b><span style="font-size: 12pt;">From:<span class="Apple-converted-space"> </span></span></b><span style="font-size: 12pt;">MITgcm-devel <mitgcm-devel-bounces@mitgcm.org> on behalf of Martin Losch <Martin.Losch@awi.de><br><b>Date:<span class="Apple-converted-space"> </span></b>Thursday, June 29, 2023 at 12:12 PM<br><b>To:<span class="Apple-converted-space"> </span></b>MITgcm-devel <MITgcm-devel@mitgcm.org><br><b>Subject:<span class="Apple-converted-space"> </span></b>Re: [MITgcm-devel] [EXTERNAL] ggl90 adjoint?<o:p></o:p></span></p></div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">Thanks Patrick and Ou,<o:p></o:p></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">that’s very helpful (so that I don’t have to debug my configuration hopelessly). Maybe I can find the time to stabilise this.<o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><o:p> </o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">How about KPP. Does that work inADmode? If so, how did you make it work (probably @Patrick)? Maybe that can be used and transfered to ggl90.<o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><o:p> </o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">Martin<o:p></o:p></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><o:p> </o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">PS. I was under the impression that ggl90 is preferred over kpp because it is more stable in ADmode, but obviously that was a misunderstanding ...<o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><br><br><o:p></o:p></div><blockquote style="margin-top: 5pt; margin-bottom: 5pt;"><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">On 29. Jun 2023, at 18:13, Wang, Ou (US 329B) <ou.wang@jpl.nasa.gov> wrote:<o:p></o:p></div></div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><o:p> </o:p></div><div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">Hi Martin,<o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"> <o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">I had also tried to do an adjoint run without turning off GGL90 in the adjoint mode a while ago, but was not successful. Similar to what you are experiencing, that run quickly crashed in the adjoint mode. I didn’t attempt to debug it at the time.<o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"> <o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">Best,<o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">Ou<o:p></o:p></div></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"> <o:p></o:p></div></div><div style="border-width: 1pt medium medium; border-style: solid none none; padding: 3pt 0in 0in; border-color: currentcolor; border-image: none;"><p class="MsoNormal" style="margin: 0in 0in 12pt 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><b><span style="font-size: 12pt;">From:<span class="apple-converted-space"> </span></span></b><span style="font-size: 12pt;">MITgcm-devel <mitgcm-devel-bounces@mitgcm.org> on behalf of Martin Losch <Martin.Losch@awi.de><br><b>Date:<span class="apple-converted-space"> </span></b>Thursday, June 29, 2023 at 8:27 AM<br><b>To:<span class="apple-converted-space"> </span></b>MITgcm-devel <MITgcm-devel@mitgcm.org><br><b>Subject:<span class="apple-converted-space"> </span></b>[EXTERNAL] [MITgcm-devel] ggl90 adjoint?</span><o:p></o:p></p></div><div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;">Hi there,<br><br>this is a question for the ECCO related crowd. Is there a successful run with pkg/ggl90 NOT turned off in ADmode? i.e. that can get along without<br> useGGL90inAdMode = .FALSE.,<br>in data.autodiff?<br><br>It keeps blowing up on me very quickly. Any insights are hightly appreciated? Did anyone try to debug this?<br><br>Martin<br>_______________________________________________<br>MITgcm-devel mailing list<br><a href="mailto:MITgcm-devel@mitgcm.org" style="color: blue; text-decoration: underline;">MITgcm-devel@mitgcm.org</a><br><a href="https://urldefense.us/v3/__http:/mailman.mitgcm.org/mailman/listinfo/mitgcm-devel__;!!PvBDto6Hs4WbVuu7!OieMzp4qr_WHIv_jv0HM9CTxq9-He6b8Qy8S2sTUa3c8Z4VUvloerrbIqgqR-rvqoN82HKMeAARjDRR9973Dqye4Hw$" style="color: blue; text-decoration: underline;">https://urldefense.us/v3/__http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel__;!!PvBDto6Hs4WbVuu7!OieMzp4qr_WHIv_jv0HM9CTxq9-He6b8Qy8S2sTUa3c8Z4VUvloerrbIqgqR-rvqoN82HKMeAARjDRR9973Dqye4Hw$</a><o:p></o:p></div></div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><span style="font-size: 9pt; font-family: Helvetica;">_______________________________________________<br>MITgcm-devel mailing list<br></span><a href="mailto:MITgcm-devel@mitgcm.org" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica;">MITgcm-devel@mitgcm.org</span></a><span style="font-size: 9pt; font-family: Helvetica;"><br></span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica;">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel</span></a><o:p></o:p></div></div></blockquote></div><div style="margin: 0in 0in 0in 0.5in; font-size: 11pt; font-family: Calibri, sans-serif;"><o:p> </o:p></div></div></div></div></div><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">MITgcm-devel mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><a href="mailto:MITgcm-devel@mitgcm.org" style="color: blue; text-decoration: underline; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;">MITgcm-devel@mitgcm.org</a><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel" style="color: blue; text-decoration: underline; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel</a></div></blockquote></div><br></div></div>_______________________________________________<br>MITgcm-devel mailing list<br>MITgcm-devel@mitgcm.org<br>http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel<br></div></blockquote></div><br></div></body></html>