<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=utf-8">
<meta name="Titre" content="">
<meta name="Mots clés" content="">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<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;}
@font-face
{font-family:SimSun;
panose-1:2 1 6 0 3 1 1 1 1 1;}
@font-face
{font-family:"MS Mincho";
panose-1:2 2 6 9 4 2 5 8 3 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
{margin:0cm;
margin-bottom:.0001pt;
font-size:12.0pt;
font-family:"Times New Roman";}
a:link, span.MsoHyperlink
{mso-style-priority:99;
color:blue;
text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
{mso-style-priority:99;
color:purple;
text-decoration:underline;}
span.apple-style-span
{mso-style-name:apple-style-span;}
span.EmailStyle18
{mso-style-type:personal;
font-family:Calibri;
color:windowtext;}
span.EmailStyle19
{mso-style-type:personal-reply;
font-family:Calibri;
color:windowtext;}
span.msoIns
{mso-style-type:export-only;
mso-style-name:"";
text-decoration:underline;
color:teal;}
.MsoChpDefault
{mso-style-type:export-only;
font-size:10.0pt;}
@page WordSection1
{size:612.0pt 792.0pt;
margin:70.85pt 70.85pt 70.85pt 70.85pt;}
div.WordSection1
{page:WordSection1;}
--></style>
</head>
<body bgcolor="white" lang="FR" link="blue" vlink="purple">
<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Thank you very much, Yu-Kun, for your suggestion. Unfortunately, it is not the origin of my problem.
<br>
My 1<sup>st</sup> grid number is the boundary, and my 76<sup>th</sup> grid number is centered exactly at the Equator.
<br>
I also checked my IC variables and those are perfectly symmetric about the equator (assuming that I am defining my initial U and V velocity at cell center, right?).<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Let me know if you have any other suggestions, and thank you again.<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Best,<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Bertrand<o:p></o:p></span></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"><o:p> </o:p></span></p>
<div style="border:none;border-top:solid #B5C4DF 1.0pt;padding:3.0pt 0cm 0cm 0cm">
<p class="MsoNormal"><b><span style="font-family:Calibri;color:black">De : </span>
</b><span style="font-family:Calibri;color:black">MITgcm-support <mitgcm-support-bounces@mitgcm.org> au nom de SYSU <qianyk@mail3.sysu.edu.cn><br>
<b>Répondre à : </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
<b>Date : </b>mercredi 15 novembre 2017 à 4:33 PM<br>
<b>À : </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
<b>Objet : </b>Re: [MITgcm-support] Non-Zero Meridional Velocity Kelvin Wave<o:p></o:p></span></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<p class="MsoNormal">Hi Bertrand, <o:p></o:p></p>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Not sure but I guess the symmetry of your grid about the equator is important.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
<div>
<p class="MsoNormal">Your Y-grid number is 150, so make sure the 75th grid is located exactly at equator, not the 75.5th (center btw 75th and 76th). Make sure the IC setup by normal mode is also symmetric about equator. That is the 1st and 149th, 2nd and
148th, ..., are exactly the same, as the 150th is the boundary.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><o:p> </o:p></p>
<div id="AppleMailSignature">
<p class="MsoNormal">------------ <o:p></o:p></p>
<div>
<p class="MsoNormal">Best regards <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"><br>
<br>
<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Yu-Kun Qian (<span style="font-family:SimSun">钱钰坤</span>)<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Center for Monsoon and Environment Research <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Department of Atmospheric Sciences<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Sun Yat-sun University<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">Guangzhou, China<o:p></o:p></p>
</div>
</div>
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt"><br>
<span style="font-family:"MS Mincho"">在</span> 2017<span style="font-family:"MS Mincho"">年</span>11<span style="font-family:"MS Mincho"">月</span>16<span style="font-family:"MS Mincho"">日，</span>03:05<span style="font-family:"MS Mincho"">，</span>Bertrand Louis
Rene Delorme <<a href="mailto:bdelorme@stanford.edu">bdelorme@stanford.edu</a>> <span style="font-family:"MS Mincho"">
写道：</span><o:p></o:p></p>
</div>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Dears,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">I am trying to set-up a very simple experiment to generate baroclinic Kelvin waves over an equatorial beta-plane.
</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">My approach so far has been to prescribe initial conditions for U, V, T and ETA based on the normal-mode solutions of the wave.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">I managed to get a nice modal structure (both in the vertical and horizontal plane), which propagates without attenuation throughout the simulation.
</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Here is my problem: in my initial conditions, I set V=0, but right after the first iteration V get non-zero values.
</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">As I am generating a Kelvin wave, I would expect V to stay equal to zero throughout the simulation, if the flow is well-balanced.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">In my case, V ends up with an amplitude approximately one order of magnitude lower than U, which is still pretty high.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">From the patterns I get in the outputs, it seems that the model is using dV/dy balance dU/dx in the continuity equation.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">As I prescribed an initial free-surface, I was expecting W to play this role alone.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">My issue can be link to an error in the normal-mode solutions that I use, which I checked carefully without finding any mistakes.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Another possibility would be that it comes from the set-up of my experiment.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">I am using a domain centered at the equator with 60x150x200 grid points, with a resolution high enough to resolve the wave.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">The only package that I am using is the MNC package, and my data file is pretty simple (see below).</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">If any of you has an idea of what I should check/test to understand where those non-zero values come from, please let me know.</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Best,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Bertrand</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&PARM01</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">sBeta=0.0,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">selectCoriMap=1,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">f0=0,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">beta=2.29E-11,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">readBinaryPrec=64,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">writeBinaryPrec=64,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&PARM02</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&PARM03</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">startTime=0.,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">endTime=1600000,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">deltaT=200,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">dumpFreq=16000,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">monitorFreq=1600000,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">monitorSelect=1,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&PARM04</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">usingCartesianGrid=.TRUE.,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">ygOrigin=-1050000,</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">delXfile='delXvar',</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">delYfile='delYvar',</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">delRfile='delZvar',</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&PARM05</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">bathyFile='topog.flat',</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">hydrogThetaFile='Tini.bin',</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">uVelInitFile='Uini.bin',</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">vVelInitFile='Vini.bin',</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">pSurfInitFile='ETAini.bin',</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">&</span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"> </span><o:p></o:p></p>
<p class="MsoNormal"><span lang="EN-US"> </span><o:p></o:p></p>
</div>
</blockquote>
<blockquote style="margin-top:5.0pt;margin-bottom:5.0pt">
<div>
<p class="MsoNormal" style="margin-bottom:12.0pt">_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><o:p></o:p></p>
</div>
</blockquote>
</div>
</div>
</body>
</html>