<div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr">Martin,</div><div dir="ltr"><br></div><div dir="ltr">    Thanks much for the reply.  The .png picture is a vertical cross section.  The initial temperature field is cold blubs at the top.  Suddenly high-temperature disturbances occur around the lateral boundary in every simulation with different grids of .bin data with which I have worked.  Thus, I believe something wrong with the lateral b.c.  If the setup for Orlanski b.c. is right, why would the temperature patterns be the same in two different b.c.?</div><div dir="ltr"><br></div><div dir="ltr">    I turned off <span class="gmail-gr_ gmail-gr_55 gmail-gr-alert gmail-gr_gramm gmail-gr_inline_cards gmail-gr_disable_anim_appear gmail-Style gmail-replaceWithoutSep" id="gmail-55" style="display:inline;border-bottom:2px solid transparent;background-repeat:no-repeat">centrifugal</span> force by setting amplCF=0.0 in data.  The cold blubs are only testing simulations.</div><div dir="ltr"><br></div><div>    I am new to MITgcm and would like to know about the relationship between radius and grid information.  Do you think the given <span class="gmail-gr_ gmail-gr_58 gmail-gr-alert gmail-gr_spell gmail-gr_inline_cards gmail-gr_disable_anim_appear gmail-ContextualSpelling gmail-ins-del gmail-multiReplace" id="gmail-58" style="display:inline;border-bottom:2px solid transparent;background-repeat:no-repeat">rSphere</span> or <span style="color:rgb(0,0,0)">radius_fromHorizGrid in the data file would change the grid information in CS, LLC, or the curvilinear coordinates when </span><font color="#000000">OLD_GRID_IO is set in CPP_OPTIONS.h</font><span style="color:inherit;font-size:inherit">?  Certainly, I will try the C</span><span class="gmail-gr_ gmail-gr_51 gmail-gr-alert gmail-gr_spell gmail-gr_inline_cards gmail-gr_disable_anim_appear gmail-ContextualSpelling gmail-ins-del gmail-multiReplace" id="gmail-51" style="display:inline;border-bottom:2px solid transparent;background-repeat:no-repeat">artesian</span> coordinates you suggested.</div><div><br></div><div>    Many thanks,</div><div><br></div><div>                Yi-Chih</div><div dir="ltr"><div><br></div><div>Date: Wed, 19 Dec 2018 13:57:24 +0100<br>From: Martin Losch <<a href="mailto:Martin.Losch@awi.de" target="_blank">Martin.Losch@awi.de</a>><br>To: MITgcm Support <<a href="mailto:mitgcm-support@mitgcm.org" target="_blank">mitgcm-support@mitgcm.org</a>><br>Subject: Re: [MITgcm-support] is Orlanski radiation condition working<br>        in the rotating rank of MITgcm?<br>Message-ID: <<a href="mailto:2ECB81F3-A681-4324-9108-D4065A7E778F@awi.de" target="_blank">2ECB81F3-A681-4324-9108-D4065A7E778F@awi.de</a>><br>Content-Type: text/plain; charset="utf-8"<br><br>Yi-Chih,<br><br>I still don?t really know what the problem is. in your png file, which I assume is a vertical section along a boundary, there are some smooth variations of the otherwise homogeneous temperature (scaled density) field. What do you expect to see? Constant temperature?<br><br>I would try without the ?centrifugal force? in apply_forcing.F, which I assume is your addition (BTW, the masking in the ?centrifugal force? is not correct. I should be maskW for gU_arr and maskS vor gV_arr, you have maskW for both), and see if you get what you expect.<br><br>More general comments:<br><br>The obcs-settings appear to be OK although I still don?t see why you want to use open boundary conditions for your simluation at all.<br><br>You seem to use z-coordinates (default). In your data you have very high viscosities (1e6), but very low diffusivities (5e-6, 1e-7). Is that on purpose? I would expect that this would put a strong contraint on your timestep deltaT. Also you use rSphere = 60000.E3, rotationPeriod=38520.0, so a planet rotating about twice as fast as the earth with a radius 10 times larger. Is the larger radius reflected in the grid input fields?<br><br>For your experiment, I would just use a cartesian grid centered at the North Pole of the sphere (constant f0), or if the variation of coriolis is important, you can read the (precomputed) coriolis parameter via files fCoriC.bin, fCorCs.bin (for the 2*om*cos(lat)),fCoriG.bin, for that you need to set selectCoriMap = 3 in data. It would make things much easier for you to control.<br><br>Martin</div></div><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div dir="ltr"><div class="gmail_quote"><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex">
</blockquote></div></div></div></div></div>
</blockquote></div></div></div></div></div></div></div></div></div></div></div></div></div></div>
</blockquote></div></div></div></div></div>