<div>Hi,</div><div><br></div><div>I notice that there are discontinuities in your surface plot of the residual of the two terms in PDF file.  It is a little surprising to me, but your PDF plots are really blurred I cannot read the numerics.</div><div><br></div><div>I do not know much about normal mode.  But the classical Kelvin wave in the equatorial beta plane is derived from a simplified dynamics e.g., ignoring the nonlinear advection terms.  So I guess you cannot expect a perfect Kelvin wave in the MITgcm as it is based on the primitive equations (there should be exceptions that some solutions perfectly fit the nonlinear terms, rendering an exact cancellation of the nonlinear terms).  You may verify your normal mode perfectly only in the simplified model that your solution is based on.</div><div><sign signid="0"><div style="color:#909090;font-family:Arial Narrow;font-size:12px"><br><br><br>------------------</div><div style="font-size:14px;font-family:Verdana;color:#000;"><div> <font size="2">Best regards <br> <br>Yu-Kun Qian (钱钰坤) <br>Center for Monsoon and Environment Research</font> </div> <font size="2">Department of Atmospheric Sciences<br>School of Environmental Science and Engineering<font size="2"> <br> </font>Sun Yat-sen University <br>No. 135 Xingang West Road, Haizhu District <br>Guangzhou, 510275, P.R. China <br>Tel; 020-84115227 <br>Email: <a href="mailto:qianyk@mail2.sysu.edu.cn">qianyk@mail3.sysu.edu.cn</a> </font> 
<div> <font size="2"> </font> </div></div></sign></div><div> </div><div><includetail><div> </div><div> </div><div style="font:Verdana normal 14px;color:#000;"><div style="FONT-SIZE: 12px;FONT-FAMILY: Arial Narrow;padding:2px 0 2px 0;">------------------ 原始邮件 ------------------</div><div style="FONT-SIZE: 12px;background:#efefef;padding:8px;"><div id="menu_sender"><b>发件人:</b> "Bertrand Louis Rene Delorme"<bdelorme@stanford.edu>;</div><div><b>发送时间:</b> 2017年11月17日(星期五) 上午10:58</div><div><b>收件人:</b> "mitgcm-support@mitgcm.org"<mitgcm-support@mitgcm.org>; <wbr></div><div></div><div><b>主题:</b> Re: [MITgcm-support] Non-Zero Meridional Velocity Kelvin Wave</div></div><div> </div><div style="position:relative;">





<style></style>


<div class="WordSection1">
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">A little update on my problem (note also that the plot of Vm_dPHdy_Surf in my previous message is wrong):<o:p></o:p></span></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo3"><!--[if !supportLists]--><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"><span style="mso-list:Ignore">-<span style="  ; ">         
</span></span></span><!--[endif]--><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">The full baroclinic pressure gradient force (ie: Vm_dPHdy-gravity*(ETAN(i) - ETAN(i-1) )/DYC)  and
<o:p></o:p></span></p>
<p class="MsoNormal" style="margin-left:18.0pt;text-indent:17.4pt"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">the Coriolis force (in Vm_Advec or Vm_Cori -- same) are the only forces occurring in the v-momentum
 equation <o:p></o:p></span></p>
<p class="MsoNormal" style="margin-left:18.0pt;text-indent:17.4pt"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">(AB_gV is non-zero but much smaller).<o:p></o:p></span></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l0 level1 lfo3"><!--[if !supportLists]--><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US"><span style="mso-list:Ignore">-<span style="  ; ">         
</span></span></span><!--[endif]--><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">I plotted those at the first time step (see enclosed). We can see that they do seem to have opposite spatial structures.
<o:p></o:p></span></p>
<p class="MsoListParagraph"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Yet, looking at the difference of their amplitude ( ABS(f*u) - ABS(dPHI/dy) ), we can see that they do not balance
<o:p></o:p></span></p>
<p class="MsoListParagraph"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">each other perfectly. The induced acceleration is consistent with the plot of V.<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>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">So, I managed to find where the nonzero values in V come from, but still struggle to know why.
<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">If anyone has an idea, please let me know.
<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>
<p class="MsoNormal"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Thanks,<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>
<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 lang="EN-US" style="font-family:Calibri;color:black">De :
</span></b><span lang="EN-US" style="font-family:Calibri;color:black">MITgcm-support <mitgcm-sup</span><span style="font-family:Calibri;color:black">port-bounces@mitgcm.org> au nom de Bertrand DELORME <bdelorme@stanford.edu><br>
<b>Répondre à : </b>"mitgcm-support@mitgcm.org" <mitgcm-support@mitgcm.org><br>
<b>Date : </b>jeudi 16 novembre 2017 à 12:20 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"><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Thank you again, Yu-Kun, for your precious help.</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 just checked the model output after the first time-step at t1 (where t0=initial conditions)., and right at t1 v is nonzero.
</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">In the v-momentum equation, the 3 diagnostics terms that I found to be nonzero at t1 are:</span><o:p></o:p></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l1 level1 lfo2"><!--[if !supportLists]--><span style="font-family:Calibri"><span style="mso-list:Ignore">-<span style="  ; ">         
</span></span></span><!--[endif]--><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">AB_gV (~10^-12)</span><o:p></o:p></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l1 level1 lfo2"><!--[if !supportLists]--><span style="font-family:Calibri"><span style="mso-list:Ignore">-<span style="  ; ">         
</span></span></span><!--[endif]--><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Vm_Advec — which is exactly equal to Vm_Cori — (~10^-9)</span><o:p></o:p></p>
<p class="MsoListParagraph" style="text-indent:-18.0pt;mso-list:l1 level1 lfo2"><!--[if !supportLists]--><span style="font-family:Calibri"><span style="mso-list:Ignore">-<span style="  ; ">         
</span></span></span><!--[endif]--><span lang="EN-US" style="font-size:11.0pt;font-family:Calibri;mso-fareast-language:EN-US">Vm_dPHdy (~10^-9)</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 also computed separately the surface pressure horizontal gradient (which is ~10^-8) using the relation:</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">Vm_dPHdy_Surf = - gravity * ( Eta(i) - Eta(i-1) ) / dyC</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 other terms (Vm_Diss and Vm_Ext) are equal to zero.</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">In the Kelvin wave case, I would expect Vm_Cori and Vm_dPHdy+ Vm_dPHdy_Surf to balance each other perfectly.</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">Yet, as you can see in the plots enclosed to this email, their spatial structures are very different.</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 shape of the Coriolis term looks consistent with the prescribed U velocity at t0.</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">And from the plot of V at t1, it seems that it is the quantity accelerating a flow in the meridional direction.</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">Surprisingly, it seems that even if Vm_Cori and Vm_dPHdy are the same magnitude, Vm_dPHdy has no influence on V.</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 second problematic thing is the fact that the baroclinic pressure gradient seems to have a different vertical
</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">wavelength than the temperature anomaly that I prescribe at t=0 (see also enclosed).</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">Yet, this gradient should depend only on that quantity right (I do not prescribe any S 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"> </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">Please let me know what your thoughts are at the view of this plot.</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 regards,</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"> </span><o:p></o:p></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 à 5:52 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</span><o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<p class="MsoNormal">Hi, <o:p></o:p></p>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">You may want to find the time v becomes nonzero.  You can output probably the 2nd step to see if v becomes nonzero, and then using diagnostics to do the full v momentum budget to find the reason.<o:p></o:p></p>
</div>
<div>
<p class="MsoNormal"> <o:p></o:p></p>
</div>
<div>
<p class="MsoNormal">A possible reason is that you use analytic normal mode as IC, if specified to c-grid mitgcm, may generate truncation error, so that at the 2nd step v becomes nonzero.<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>
<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="  ; ">在</span> 2017<span style="  ; ">年</span>11<span style="  ; ">月</span>16<span style="  ; ">日，</span>09:27<span style="  ; ">，</span>Bertrand Louis
 Rene Delorme <<a href="mailto:bdelorme@stanford.edu">bdelorme@stanford.edu</a>> <span style="  ; ">
写道：</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">Thank you very much, Yu-Kun, for your suggestion. Unfortunately, it is not the origin of my problem.
<br>
My 1^st grid number is the boundary, and my 76^th 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?).</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">Let me know if you have any other suggestions, and thank you again.</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>
<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 <<a href="mailto:mitgcm-support-bounces@mitgcm.org">mitgcm-support-bounces@mitgcm.org</a>> au nom de SYSU <<a href="mailto:qianyk@mail3.sysu.edu.cn">qianyk@mail3.sysu.edu.cn</a>><br>
<b>Répondre à : </b>"<a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a>" <<a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a>><br>
<b>Date : </b>mercredi 15 novembre 2017 à 4:33 PM<br>
<b>À : </b>"<a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a>" <<a href="mailto:mitgcm-support@mitgcm.org">mitgcm-support@mitgcm.org</a>><br>
<b>Objet : </b>Re: [MITgcm-support] Non-Zero Meridional Velocity Kelvin Wave</span><o:p></o:p></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>
<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="  ; ">在</span> 2017<span style="  ; ">年</span>11<span style="  ; ">月</span>16<span style="  ; ">日，</span>03:05<span style="  ; ">，</span>Bertrand Louis
 Rene Delorme <<a href="mailto:bdelorme@stanford.edu">bdelorme@stanford.edu</a>> <span style="  ; ">
写道：</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>
</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>


</div></div><!--<![endif]--></includetail></div>