<html xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<!--[if !mso]><style>v\:* {behavior:url(#default#VML);}
o\:* {behavior:url(#default#VML);}
w\:* {behavior:url(#default#VML);}
.shape {behavior:url(#default#VML);}
</style><![endif]--><style><!--
/* Font Definitions */
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:DengXian;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
@font-face
{font-family:Aptos;
panose-1:2 11 0 4 2 2 2 2 2 4;}
@font-face
{font-family:"Segoe UI";
panose-1:2 11 5 2 4 2 4 2 2 3;}
@font-face
{font-family:"\@DengXian";
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"\@SimSun";
panose-1:2 1 6 0 3 1 1 1 1 1;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
font-size:12.0pt;
font-family:"Aptos",sans-serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
p.xmsonormal, li.xmsonormal, div.xmsonormal
{mso-style-name:x_msonormal;
margin:0cm;
font-size:12.0pt;
font-family:"Aptos",sans-serif;}
span.EmailStyle21
{mso-style-type:personal-reply;
font-family:"Aptos",sans-serif;
color:windowtext;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;
mso-ligatures:none;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:72.0pt 72.0pt 72.0pt 72.0pt;}
div.WordSection1
{page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
</head>
<body lang="EN-GB" link="blue" vlink="purple" style="word-wrap:break-word">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt">Hi Flynn,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">You are right, thank you! <o:p>
</o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt">I found the ‘new wet grid-pts ‘ you mentioned and extrapolated the temperature and salinity, the error is solved.<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div>
<p class="MsoNormal" style="margin-bottom:6.0pt"><b><span style="font-size:11.0pt;font-family:"Arial",sans-serif;color:#007FA3"><o:p> </o:p></span></b></p>
<p class="MsoNormal" style="margin-bottom:6.0pt"><b><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#007FA3">Kind Regards<o:p></o:p></span></b></p>
<p class="MsoNormal" style="margin-bottom:6.0pt"><b><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:#1E3765">Juncheng Zhang</span></b><span style="font-size:10.0pt;font-family:"Arial",sans-serif"><br>
<span style="color:#007FA3">PhD Student, <o:p></o:p></span></span></p>
<p class="MsoNormal" style="margin-bottom:6.0pt"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#007FA3">Earth, Atmospheric, and Planetary Physics</span><span style="font-size:10.0pt;font-family:"Arial",sans-serif"><br>
<br>
<span style="color:#1E3765">Department of Physics<br>
University of Toronto<br>
McLennan Physical Laboratories<br>
60 St Gorge St<br>
Toronto, Ontario, M5S 1A7<br>
<br>
</span><a href="mailto:juncheng.zhang@mail.utoronto.ca"><span style="color:#0563C1">juncheng.zhang@mail.utoronto.ca</span></a><br>
<br>
</span><a href="https://www.physics.utoronto.ca/"><b><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#007894;text-decoration:none"><img border="0" width="277" height="62" style="width:2.8854in;height:.6458in" id="Picture_x0020_1" src="cid:image002.jpg@01DBDB25.D0EF05C0" alt="Pursue STEM"></span></b></a><span style="font-size:10.0pt;font-family:"Arial",sans-serif"><o:p></o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt"><o:p> </o:p></span></p>
<div id="mail-editor-reference-message-container">
<div>
<div>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="margin-bottom:12.0pt"><b><span style="color:black">From:
</span></b><span style="color:black">MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Flynn Ames <f.ames@pgr.reading.ac.uk><br>
<b>Date: </b>Tuesday, 10 June 2025 at 04:09<br>
<b>To: </b>mitgcm-support@mitgcm.org <mitgcm-support@mitgcm.org><br>
<b>Subject: </b>Re: [MITgcm-support] Interpolation issue: wet grid-pts with theta=0 identically<o:p></o:p></span></p>
</div>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" align="left" width="100%" style="width:100.0%;background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important;border:0!important;display:table!important;width:100%!important;table-layout:fixed!important;border-collapse:seperate!important;float:none!important;border-spacing:0px 0px!important">
<tbody>
<tr style="background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important">
<td style="background:#A6A6A6;padding:0cm 0cm 0cm 0cm;background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important;padding:2px!important;background-color:#A6A6A6!important;width:0px!important">
</td>
<td width="100%" style="width:100.0%;background:#EAEAEA;padding:0cm 0cm 0cm 0cm;background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important;width:100%!important;background-color:#EAEAEA!important;padding:15px!important;font-size:12px!important;font-weight:normal!important;color:#212121!important;text-align:left!important;word-wrap:break-word!important">
<div>
<p class="MsoNormal" style="mso-element:frame;mso-element-frame-hspace:2.25pt;mso-element-wrap:around;mso-element-anchor-vertical:paragraph;mso-element-anchor-horizontal:column;mso-height-rule:exactly">
<span style="font-family:"Segoe UI",sans-serif;color:black">You don't often get email from f.ames@pgr.reading.ac.uk.
<a href="https://aka.ms/LearnAboutSenderIdentification">Learn why this is important</a>
</span><span style="font-family:"Segoe UI",sans-serif"><o:p></o:p></span></p>
</div>
</td>
<td style="background:#EAEAEA;padding:0cm 0cm 0cm 0cm;background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important;width:75px!important;background-color:#EAEAEA!important;padding:5px!important;font-size:12px!important;font-weight:normal!important;color:#212121!important;text-align:left!important;word-wrap:break-word!important">
</td>
</tr>
</tbody>
</table>
<div>
<div>
<p class="MsoNormal"><span style="color:black">Hi Juncheng,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">The error might be occurring because, as you increase resolution while maintaining the same bathymetry profile, you will resolve more of the detail of the bathymetry profile.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">As a result, at the bottom boundary in the 1km simulation, you have 'new' wet grid points that were not present in the 5km simulation (i.e., vertical levels at some lat-long positions that were previously dry will
now be wet).<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">My guess is that when you interpolate to 1km resolution, you aren't filling the 'new' wet grid points with any non-zero value values, which is why the model throws the error.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">When ignoring the error, you are initialising the model with a temperature of zero degrees at 2981 grid points at the bottom boundary.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">If you are you are providing other initial conditions (e.g., velocity fields), you will also likely be initialising these with values of zero. This could potentially be causing the numerical artifacts.<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">Admittedly I'm not sure what the best way round this. Perhaps creating a 'buffer layer' at the bottom of your interpolated fields (e.g., where you take the interpolated temperature field and create an extra layer
of non-zero values at the ocean bottom, perhaps just copy-pasting the value of temperature from the grid cell above) could be something to try?<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black"><o:p> </o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">Cheers,<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><span style="color:black">Flynn<o:p></o:p></span></p>
</div>
<div class="MsoNormal" align="center" style="text-align:center">
<hr size="2" width="98%" align="center">
</div>
<div id="divRplyFwdMsg">
<p class="MsoNormal"><b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black">From:</span></b><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:black"> MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of
Juncheng Zhang <juncheng.zhang@mail.utoronto.ca><br>
<b>Sent:</b> 09 June 2025 5:49 PM<br>
<b>To:</b> MITgcm Support <mitgcm-support@mitgcm.org><br>
<b>Subject:</b> [MITgcm-support] Interpolation issue: wet grid-pts with theta=0 identically</span>
<o:p></o:p></p>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
</div>
<div>
<table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" align="left" width="100%" style="width:100.0%;background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important;border:0!important;display:table!important;width:100%!important;table-layout:fixed!important;border-collapse:seperate!important;float:none!important;border-spacing:0px 0px!important">
<tbody>
<tr style="background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important">
<td style="background:#A6A6A6;padding:0cm 0cm 0cm 0cm;background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important;padding:2px!important;background-color:#A6A6A6!important;width:0px!important">
</td>
<td width="100%" style="width:100.0%;background:#EAEAEA;padding:0cm 0cm 0cm 0cm;background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important;width:100%!important;background-color:#EAEAEA!important;padding:15px!important;font-size:12px!important;font-weight:normal!important;color:#212121!important;text-align:left!important;word-wrap:break-word!important">
<div>
<p class="MsoNormal" style="mso-element:frame;mso-element-frame-hspace:2.25pt;mso-element-wrap:around;mso-element-anchor-vertical:paragraph;mso-element-anchor-horizontal:column;mso-height-rule:exactly">
<span style="font-family:"Segoe UI",sans-serif;color:black">You don't often get email from juncheng.zhang@mail.utoronto.ca.
<a href="https://aka.ms/LearnAboutSenderIdentification">Learn why this is important</a>
</span><span style="font-family:"Segoe UI",sans-serif"><o:p></o:p></span></p>
</div>
</td>
<td style="background:#EAEAEA;padding:0cm 0cm 0cm 0cm;background:revert!important;border:revert!important;bottom:revert!important;color:revert!important;direction:revert!important;display:revert!important;font-size:revert!important;height:revert!important;letter-spacing:revert!important;line-height:revert!important;margin:revert!important;opacity:revert!important;order:revert!important;outline:revert!important;overflow:revert!important;padding:revert!important;position:revert!important;tab-size:revert!important;table-layout:revert!important;text-align:revert!important;text-indent:revert!important;text-orientation:revert!important;text-overflow:revert!important;text-transform:revert!important;top:revert!important;vertical-align:revert!important;visibility:revert!important;white-space:revert!important;width:revert!important;word-break:revert!important;word-spacing:revert!important;writing-mode:revert!important;zoom:revert!important;width:75px!important;background-color:#EAEAEA!important;padding:5px!important;font-size:12px!important;font-weight:normal!important;color:#212121!important;text-align:left!important;word-wrap:break-word!important">
</td>
</tr>
</tbody>
</table>
<div>
<div>
<p class="xmsonormal"><span style="font-size:11.0pt">Dear all, </span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt">I am using MITgcm to run a 1km high-resolution simulation based on the tutorial reentrant channel,
</span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt">I first run a 5 km simulation, then use the output of 5km run as the initial condition of my 1km run,
</span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt">However, after I did the interpolation of 5 km output, and run the model, the error said below:</span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt">(PID.TID 0000.0001) *** ERROR *** INI_THETA: found 2891 wet grid-pts with theta=0 identically.</span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt">(PID.TID 0000.0001) *** ERROR *** If this is intentional, you need to set checkIniTemp=.FALSE. in "data", namelist PARM05</span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt">I know wet grid-pts relates to the parameter hFac? But I don’t know how to fix it,
</span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt">And I tried to ignore this error by adding checkIniTemp=.FALSE., however, the simulation is not stable, and I had weird numerical artifacts appearing.</span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt"> </span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt">Thank you in advance, I appreciate any help.</span><o:p></o:p></p>
<p class="xmsonormal"><span style="font-size:11.0pt"> </span><o:p></o:p></p>
<div>
<p class="xmsonormal" style="margin-bottom:6.0pt"><b><span style="font-size:11.0pt;font-family:"Arial",sans-serif;color:#007FA3"> </span></b><o:p></o:p></p>
<p class="xmsonormal" style="margin-bottom:6.0pt"><b><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#007FA3">Kind Regards</span></b><o:p></o:p></p>
<p class="xmsonormal" style="margin-bottom:6.0pt"><b><span style="font-size:10.5pt;font-family:"Arial",sans-serif;color:#1E3765">Juncheng Zhang</span></b><span style="font-size:10.0pt;font-family:"Arial",sans-serif"><br>
<span style="color:#007FA3">PhD Student, </span></span><o:p></o:p></p>
<p class="xmsonormal" style="margin-bottom:6.0pt"><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#007FA3">Earth, Atmospheric, and Planetary Physics</span><span style="font-size:10.0pt;font-family:"Arial",sans-serif"><br>
<br>
<span style="color:#1E3765">Department of Physics<br>
University of Toronto<br>
McLennan Physical Laboratories<br>
60 St Gorge St<br>
Toronto, Ontario, M5S 1A7<br>
<br>
</span><a href="mailto:juncheng.zhang@mail.utoronto.ca"><span style="color:#0563C1">juncheng.zhang@mail.utoronto.ca</span></a><br>
<br>
</span><a href="https://www.physics.utoronto.ca/"><span style="color:windowtext;text-decoration:none"><b><span style="font-size:10.0pt;font-family:"Arial",sans-serif;color:#007894"><img border="0" width="277" height="62" style="width:2.8854in;height:.6458in" id="x_Picture_x0020_1" src="cid:image001.jpg@01DBD93B.AB0C00B0" alt="Pursue STEM"></span></b></span></a><o:p></o:p></p>
</div>
<p class="xmsonormal"> <o:p></o:p></p>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</div>
</body>
</html>