<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=iso-8859-1">
<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:"Cambria Math";
panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
{font-family:Calibri;
panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:#0563C1;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:#954F72;
text-decoration:underline;}
p
{mso-style-priority:99;
margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
pre
{mso-style-priority:99;
mso-style-link:"HTML Preformatted Char";
margin:0cm;
margin-bottom:.0001pt;
font-size:10.0pt;
font-family:"Courier New";}
p.msonormal0, li.msonormal0, div.msonormal0
{mso-style-name:msonormal;
margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman",serif;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:"Calibri",sans-serif;
color:#1F497D;}
span.HTMLPreformattedChar
{mso-style-name:"HTML Preformatted Char";
mso-style-priority:99;
mso-style-link:"HTML Preformatted";
font-family:"Courier New";}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@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-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Dear all,<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">I got the warning and found these lines about KPP and the use of the non-hydrostatic capabilities of MITgcm (see below).<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">I’m new with MITgcm. We are trying to run MITgcm in a small lake (2.5 km in length and 0.5 km wide) and we would like to resolve the convective cells developing
during night-time cooling in the lake. The grid size that it’s needed to resolve or at least partially resolve these cells (dx = dy = 5m), taking into account non-hydrostatic effects, is still coarse (compared to the grid size that would be needed for a LES),
so, a turbulence closure is still needed to model the mixing. <o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">So please, could you explain why the non-hydrostatic capabilities cannot be used with a turbulence closure scheme? Or is this only specifically related to the
KPP formulation?<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Thanks in advance<o:p></o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D"><o:p> </o:p></span></p>
<p class="MsoNormal"><span style="font-size:11.0pt;font-family:"Calibri",sans-serif;color:#1F497D">Cintia<o:p></o:p></span></p>
<p class="MsoNormal"><o:p> </o:p></p>
<div id="divtagdefaultwrapper">
<div>
<div>
<div class="MsoNormal" align="center" style="text-align:center">
<hr size="5" width="100%" align="center">
</div>
<pre>Hi Sherry,<o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre>thanks to version control, you can easily checkout old versions of the code and try to find, when (in “code-time”) the instabilities occur. 2009 ended with checkpoint62 (get it with "cvs co -r checkpoint62 MITgcm”, but you can also use dates, see cvs documentation, see MITgcm/doc/tag-index for record of what has been done over time). With “divide and conquer” (e.g. if 2009 code works try 2012 code, if that works 2014, if not 2010.5, etc), it doesn’t take that long, but it’s tedious anyway. Once you’ve found the point, where the model starts to have problems, you can use <a href="http://mitgcm.org/viewvc/MITgcm/MITgcm/">http://mitgcm.org/viewvc/MITgcm/MITgcm/</a> to look at differences.<o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre>KPP does not make too much sense together with a non-hydrostratic simulation, where you (usually) want to resolve the mixing processes, rather than parameterise them. I would take the warnings seriously.<o:p></o:p></pre>
<pre><o:p> </o:p></pre>
<pre>Martin<o:p></o:p></pre>
<pre>><i> On 18 Aug 2016, at 07:59, Sherry <<a href="http://mitgcm.org/mailman/listinfo/mitgcm-support">schou at hawaii.edu</a>> wrote:<o:p></o:p></i></pre>
<pre>><i> <o:p></o:p></i></pre>
<pre>><i> Hi all,<o:p></o:p></i></pre>
<pre>><i> <o:p></o:p></i></pre>
<pre>><i> I have a 2-D nonhydrostatic configuration of an internal wave beam forced through specifying u, v, w, and temperature at the western boundary. At the moment I am trying to reproduce an experiment that was done in 2009 (Grisouard and Staquet, 2010), but for some reason I am getting instabilities after about 10 wave cycles, even though I am using the same input parameters. Could there be some changes in the model since 2009 that would affect default values? Any advice about how to track these down?<o:p></o:p></i></pre>
<pre>><i> <o:p></o:p></i></pre>
<pre>><i> To control the instabilities I am trying higher diffusion values, as well as turning on KPP. However, it seems there is some internal inconsistency with the configuration being nonhydrostatic, and I get the following warnings (from config_check.F):<o:p></o:p></i></pre>
<pre>><i> <o:p></o:p></i></pre>
<pre>><i> WARNING CONFIG_CHECK: Implicit viscosity applies to provisional u,vVel<o:p></o:p></i></pre>
<pre>><i> WARNING => not consistent withfinal vertical shear (after appling 3-D solver solution<o:p></o:p></i></pre>
<pre>><i> WARNING CONFIG_CHECK: Implicit viscosity not implemented in CALC_GW<o:p></o:p></i></pre>
<pre>><i> WARNING CONFIG_CHECK: Explicit viscosity might become unstable if too large<o:p></o:p></i></pre>
<pre>><i> <o:p></o:p></i></pre>
<pre>><i> Should I be concerned about these warnings? Is there any way around them since nonhydrostatic requires 3-D solver and KPP requires implicit viscosity?<o:p></o:p></i></pre>
<pre>><i> <o:p></o:p></i></pre>
<pre>><i> Thanks for reading!<o:p></o:p></i></pre>
<pre>><i> <o:p></o:p></i></pre>
<pre>><i> Sherry Chou<o:p></o:p></i></pre>
<pre>><i> University of Hawaii<o:p></o:p></i></pre>
<pre>><i> _______________________________________________<o:p></o:p></i></pre>
<pre>><i> MITgcm-support mailing list<o:p></o:p></i></pre>
<pre>><i> <a href="http://mitgcm.org/mailman/listinfo/mitgcm-support">MITgcm-support at mitgcm.org</a><o:p></o:p></i></pre>
<pre>><i> <a href="http://mitgcm.org/mailman/listinfo/mitgcm-support">http://mitgcm.org/mailman/listinfo/mitgcm-support</a><o:p></o:p></i></pre>
<pre><o:p> </o:p></pre>
<pre><o:p> </o:p></pre>
<p class="MsoNormal"><span style="font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
<p class="MsoNormal"><span style="font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
<p><span style="font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
<p><span style="font-family:"Calibri",sans-serif;color:black"><o:p> </o:p></span></p>
</div>
</div>
</body>
</html>