<html 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=gb2312">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<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: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;
        margin-bottom:.0001pt;
        text-align:justify;
        text-justify:inter-ideograph;
        font-size:10.5pt;
        font-family:DengXian;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:blue;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:#954F72;
        text-decoration:underline;}
p.MsoListParagraph, li.MsoListParagraph, div.MsoListParagraph
        {mso-style-priority:34;
        margin:0cm;
        margin-bottom:.0001pt;
        text-align:justify;
        text-justify:inter-ideograph;
        text-indent:21.0pt;
        font-size:10.5pt;
        font-family:DengXian;}
.MsoChpDefault
        {mso-style-type:export-only;}
/* Page Definitions */
@page WordSection1
        {size:612.0pt 792.0pt;
        margin:72.0pt 90.0pt 72.0pt 90.0pt;}
div.WordSection1
        {page:WordSection1;}
/* List Definitions */
@list l0
        {mso-list-id:87888683;
        mso-list-type:hybrid;
        mso-list-template-ids:358015150 -1 67698713 67698715 67698703 67698713 67698715 67698703 67698713 67698715;}
@list l0:level1
        {mso-level-tab-stop:none;
        mso-level-number-position:left;
        margin-left:18.0pt;
        text-indent:-18.0pt;}
@list l0:level2
        {mso-level-number-format:alpha-lower;
        mso-level-text:"%2\)";
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        margin-left:42.0pt;
        text-indent:-21.0pt;}
@list l0:level3
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        margin-left:63.0pt;
        text-indent:-21.0pt;}
@list l0:level4
        {mso-level-tab-stop:none;
        mso-level-number-position:left;
        margin-left:84.0pt;
        text-indent:-21.0pt;}
@list l0:level5
        {mso-level-number-format:alpha-lower;
        mso-level-text:"%5\)";
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        margin-left:105.0pt;
        text-indent:-21.0pt;}
@list l0:level6
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        margin-left:126.0pt;
        text-indent:-21.0pt;}
@list l0:level7
        {mso-level-tab-stop:none;
        mso-level-number-position:left;
        margin-left:147.0pt;
        text-indent:-21.0pt;}
@list l0:level8
        {mso-level-number-format:alpha-lower;
        mso-level-text:"%8\)";
        mso-level-tab-stop:none;
        mso-level-number-position:left;
        margin-left:168.0pt;
        text-indent:-21.0pt;}
@list l0:level9
        {mso-level-number-format:roman-lower;
        mso-level-tab-stop:none;
        mso-level-number-position:right;
        margin-left:189.0pt;
        text-indent:-21.0pt;}
ol
        {margin-bottom:0cm;}
ul
        {margin-bottom:0cm;}
--></style>
</head>
<body lang="ZH-CN" link="blue" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US">Thanks for your reply,</span><span lang="EN-US" style="font-size:10.0pt"> Jean-Michel.</span></p>
<p class="MsoListParagraph" style="margin-left:18.0pt;text-indent:-18.0pt;mso-list:l0 level1 lfo1">
<![if !supportLists]><span lang="EN-US"><span style="mso-list:Ignore">1.<span style="font:7.0pt "Times New Roman"">    
</span></span></span><![endif]><span lang="EN-US">Sorry, I didn’t understand what you mean by “Delta t”. I think that u** and v** is a correction of u* and v* by surface elevation, which still has a dimention of velocity (m/s) , but what is needed on the right
 side of Eq.(2.68) is a dimention of acceleration (m/s^2). Did you mean \partial_t?</span></p>
<p class="MsoListParagraph" style="margin-left:18.0pt;text-indent:-18.0pt;mso-list:l0 level1 lfo1">
<![if !supportLists]><span lang="EN-US"><span style="mso-list:Ignore">2.<span style="font:7.0pt "Times New Roman"">    
</span></span></span><![endif]><span lang="EN-US">I think to MITgcm, it is more safe to calulate by continuity, because nonhydrostatic pressure is solved by a method of iteration, errors can always exist. However, continuity is a strong constraint to stability,
 compared to forward integration. Just personal idea.</span></p>
<p class="MsoNormal"><span lang="EN-US">Regards,</span></p>
<p class="MsoNormal"><span lang="EN-US">Pengyang Song</span></p>
<p class="MsoNormal"><span lang="EN-US"><o:p> </o:p></span></p>
<div style="mso-element:para-border-div;border:none;border-top:solid #E1E1E1 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal" style="border:none;padding:0cm"><b>发件人<span lang="EN-US">: </span>
</b><span lang="EN-US"><a href="mailto:mitgcm-support-request@mitgcm.org">mitgcm-support-request@mitgcm.org</a><br>
</span><b>发送时间<span lang="EN-US">: </span></b><span lang="EN-US">2017</span>年<span lang="EN-US">9</span>月<span lang="EN-US">23</span>日<span lang="EN-US"> 0:00<br>
</span><b>收件人<span lang="EN-US">: </span></b><span lang="EN-US"><a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a><br>
</span><b>主题<span lang="EN-US">: </span></b><span lang="EN-US">MITgcm-support Digest, Vol 171, Issue 9</span></p>
</div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:SimSun"><o:p> </o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:10.0pt">Send MITgcm-support mailing list submissions to<br>
        mitgcm-support@mitgcm.org<br>
<br>
To subscribe or unsubscribe via the World Wide Web, visit<br>
        <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
or, via email, send a message with subject or body 'help' to<br>
        mitgcm-support-request@mitgcm.org<br>
<br>
You can reach the person managing the list at<br>
        mitgcm-support-owner@mitgcm.org<br>
<br>
When replying, please edit your Subject line so it is more specific<br>
than "Re: Contents of MITgcm-support digest..."<br>
<br>
<br>
Today's Topics:<br>
<br>
   1. Re: questions about non-hydrostiatic module (Jean-Michel Campin)<br>
<br>
<br>
----------------------------------------------------------------------<br>
<br>
Message: 1<br>
Date: Thu, 21 Sep 2017 18:50:03 -0400<br>
From: Jean-Michel Campin <jmc@mit.edu><br>
To: mitgcm-support@mitgcm.org<br>
Subject: Re: [MITgcm-support] questions about non-hydrostiatic module<br>
Message-ID: <20170921225003.GA6558@ocean.mit.edu><br>
Content-Type: text/plain; charset=iso-8859-1<br>
<br>
Hi Pengyang,<br>
<br>
Regarding question 1, I think a "Delta t" is missing there.<br>
Thanks for pointing to this.<br>
But probably better to keep u**, v** & w** (not wrong) to remain<br>
closer to what is done in the code.<br>
<br>
And regarding question 2, the fluid is incompressible, so the <br>
continuity equation has to be satisfied as well as the prognostic equation for w(n+1).<br>
In the analytical form, the system of equations guarantees both.<br>
And indeed, when/if the solver converges perfectly, the two equations<br>
give the same answer.<br>
<br>
Cheers,<br>
Jean-Michel<br>
<br>
On Thu, Sep 21, 2017 at 07:45:58AM +0000, ??? ?????? wrote:<br>
> Hi all,<br>
> <br>
> I am reading the manual of MITgcm about non-hydrostatic part. Here are <br>
> some questions:<br>
> <br>
> <a href="http://mitgcm.org/public/r2_manual/latest/online_documents/node39.html">
http://mitgcm.org/public/r2_manual/latest/online_documents/node39.html</a><br>
> <br>
> 1.? In equation 2.68, is it written by mistake at the right side? By <br>
> Marshall et al.(1997), it should be Gu, Gv, Gw instead of u**, v**, w**.<br>
> <br>
> 2.? In equation 2.71, in non-hydrostatic scheme, w should be prognostic <br>
> instead of diagnostic, should it be something like w(n+1)=w(*)+dt * ( <br>
> \partial_z p_nh ) ?<br>
> <br>
> Regards,<br>
> <br>
> Pengyang Song<br>
> <br>
> Ocean University of China<br>
> <br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> MITgcm-support@mitgcm.org<br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
<br>
------------------------------<br>
<br>
Subject: Digest Footer<br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
MITgcm-support@mitgcm.org<br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
<br>
<br>
------------------------------<br>
<br>
End of MITgcm-support Digest, Vol 171, Issue 9<br>
**********************************************<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:12.0pt;font-family:SimSun"><o:p> </o:p></span></p>
</div>
</body>
</html>