<html><head><meta http-equiv="content-type" content="text/html; charset=utf-8"></head><body style="overflow-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;">Hi Elian,<div><br></div><div>this is odd, as I had thought that the tutorial_baroclinic_gyre experiment was well tuned.</div><div><br></div><div>I had a look at the verification experiment and it looks like there are unresolved gravity waves emerging from the boundaries that fill the domain: Even after turning off momentum and tracer advection, these features remain.</div><div><br></div><div>I tried using very large horizontal viscosity (viscAhGrid = 0.1 or viscA4Grid = 0.1) and it didn’t go away until (and this is the most suprising part) I modified the grid spacing. As I said before, the acutal dx in meters in spherical grids is really the prescribed delta in degrees (here 1) times cos(lat), this means that the grid cells are not nearly quadratic in higher latitudes (here the grid spans from 14 to 76degN). I changed the delY to include the cos factor like this (python):</div><div><p style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-variant-emoji: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span style="font-variant-ligatures: no-common-ligatures;">yg=[</span><span style="font-variant-ligatures: no-common-ligatures; color: rgb(47, 180, 29);">14</span><span style="font-variant-ligatures: no-common-ligatures;">]</span></p>
<p style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-variant-emoji: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal; color: rgb(47, 180, 29);"><span style="font-variant-ligatures: no-common-ligatures; color: #2d961e"><b>for</b></span><span style="font-variant-ligatures: no-common-ligatures; color: #000000"> k </span><span style="font-variant-ligatures: no-common-ligatures; color: #c200ff"><b>in</b></span><span style="font-variant-ligatures: no-common-ligatures; color: #000000"> </span><span style="font-variant-ligatures: no-common-ligatures; color: #2d961e">range</span><span style="font-variant-ligatures: no-common-ligatures; color: #000000">(</span><span style="font-variant-ligatures: no-common-ligatures">62</span><span style="font-variant-ligatures: no-common-ligatures; color: #000000">):</span></p>
<p style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-variant-emoji: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><span style="font-variant-ligatures: no-common-ligatures; color: #2fb41d"> </span><span style="font-variant-ligatures: no-common-ligatures"> yg.append(yg[-</span><span style="font-variant-ligatures: no-common-ligatures; color: #2fb41d">1</span><span style="font-variant-ligatures: no-common-ligatures">]+cos(yg[-</span><span style="font-variant-ligatures: no-common-ligatures; color: #2fb41d">1</span><span style="font-variant-ligatures: no-common-ligatures">]*pi/</span><span style="font-variant-ligatures: no-common-ligatures; color: #2fb41d">180</span><span style="font-variant-ligatures: no-common-ligatures">))</span></p><p style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-variant-emoji: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><br></p><p style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-variant-emoji: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;">delY = np.diff(np.asarray(yg))</p><p style="margin: 0px; font-style: normal; font-variant-caps: normal; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-variant-emoji: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"><br></p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> delY=0.97029573, 0.96605988, 0.96156728, 0.95682412, 0.95183686,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.94661222, 0.94115716, 0.93547886, 0.92958469, 0.92348219,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.91717905, 0.9106831 , 0.90400227, 0.89714458, 0.89011811,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.88293099, 0.87559138, 0.86810744, 0.86048733, 0.85273916,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.844871 , 0.83689089, 0.82880674, 0.82062641, 0.81235764,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.80400804, 0.79558512, 0.78709623, 0.77854857, 0.7699492 ,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.76130498, 0.75262264, 0.7439087 , 0.73516951, 0.72641122,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.71763981, 0.70886103, 0.70008046, 0.69130348, 0.68253526,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.67378076, 0.66504477, 0.65633185, 0.64764638, 0.63899254,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.63037432, 0.6217955 , 0.61325969, 0.6047703 , 0.59633058,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.58794358, 0.57961218, 0.5713391 , 0.56312688, 0.55497789,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.54689438, 0.53887841, 0.53093191, 0.52305665, 0.51525428,</p><p style="margin: 0px; font-stretch: normal; font-size: 11px; line-height: normal; font-family: Menlo; font-size-adjust: none; font-kerning: auto; font-variant-alternates: normal; font-variant-ligatures: normal; font-variant-numeric: normal; font-variant-east-asian: normal; font-variant-position: normal; font-feature-settings: normal; font-optical-sizing: auto; font-variation-settings: normal;"> 0.50752631, 0.49987412,</p><div><br></div><div>And reran with large viscosity (viscA4Grid = 0.1 or viscAhGrid) and the noise goes away (except at the boundary itself). Maybe that’s a way forward for you. Note, that I ran this only for 100 days.</div><div><br></div><div>At the same time, we should inspect this and maybe explain this in the tutorial description (@Jeff)?</div><div><br></div><div>Martin<br><br><blockquote type="cite"><div>On 9. Apr 2025, at 10:42, Vanderborght, E.Y.P. (Elian) <e.y.p.vanderborght@uu.nl> wrote:</div><br class="Apple-interchange-newline"><div><meta charset="UTF-8"><div class="WordSection1" style="page: WordSection1; caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span style="font-size: 11pt;"><o:p> </o:p></span></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span style="font-size: 11pt;">Dear Martin and Jody,<o:p></o:p></span></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span style="font-size: 11pt;"><o:p> </o:p></span></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span lang="EN-US" style="font-size: 11pt;">Thank you for your suggestions and help. I have simplified the problem, removing practically everything and adding all processes step by step. The moment I obtain a non-zero flow field, the irregularities appear. I realized that this may just be an error in my source code, so I reinstalled the MITgcm and tried running again. Sadly, the error persisted. As a check I ran the baroclinic ocean gyre tutorial (outlined in<a href="https://mitgcm.readthedocs.io/en/latest/examples/baroclinic_gyre/baroclinic_gyre.html" style="color: blue; text-decoration: underline;">https://mitgcm.readthedocs.io/en/latest/examples/baroclinic_gyre/baroclinic_gyre.html</a>). While I do obtain similar results as outlined by the figures 4, a closer look at the velocity field reveals that they are again irregular (see attached figure)… At this point the only thing I can think of is that something goes wrong during compilation or the post-processing. But I am curious for your suggestions.<o:p></o:p></span></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span lang="EN-US"><o:p> </o:p></span></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span lang="EN-US">Kind regards,<o:p></o:p></span></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span lang="EN-US">ELian<o:p></o:p></span></div><div><div style="border-width: 1pt medium medium; border-style: solid none none; border-color: rgb(225, 225, 225) currentcolor currentcolor; border-image: none; padding: 3pt 0cm 0cm;"><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><b><span style="font-size: 11pt; font-family: Calibri, sans-serif;">Van:</span></b><span style="font-size: 11pt; font-family: Calibri, sans-serif;"><span class="Apple-converted-space"> </span>MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" style="color: blue; text-decoration: underline;">mitgcm-support-bounces@mitgcm.org</a>><span class="Apple-converted-space"> </span><b>Namens<span class="Apple-converted-space"> </span></b>Jody Klymak<br><b>Verzonden:</b><span class="Apple-converted-space"> </span>dinsdag 8 april 2025 15:53<br><b>Aan:</b><span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;">mitgcm-support@mitgcm.org</a><br><b>Onderwerp:</b><span class="Apple-converted-space"> </span>Re: [MITgcm-support] Irregular Velocity Fields Near Basin Boundaries in MITgcm<o:p></o:p></span></div></div></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div><table class="MsoNormalTable" border="0" cellspacing="0" cellpadding="0" align="left" width="100%" style="width: 1007px;"><tbody><tr><td width="0" style="width: 0.3pt; background: rgb(166, 166, 166); padding: 5.25pt 1.5pt;"></td><td width="100%" style="aspect-ratio: revert !important; background: revert !important; block-size: revert !important; border: revert !important; bottom: revert !important; color: revert !important; color-scheme: revert !important; content-visibility: revert !important; cursor: revert !important; direction: revert !important; display: revert !important; font-size: revert !important; height: revert !important; hyphens: revert !important; letter-spacing: revert !important; line-height: revert !important; margin: revert !important; opacity: revert !important; order: revert !important; outline: revert !important; overflow: revert !important; padding: revert !important; position: revert !important; resize: revert !important; rotate: revert !important; scale: revert !important; tab-size: revert !important; table-layout: revert !important; text-align: revert !important; text-indent: revert !important; text-orientation: revert !important; text-overflow: revert !important; text-shadow: revert !important; text-transform: revert !important; text-wrap: revert !important; top: revert !important; transition: revert !important; vertical-align: revert !important; visibility: revert !important; white-space: revert !important; width: revert !important; word-break: revert !important; word-spacing: revert !important; writing-mode: revert !important; zoom: revert !important;"><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span style="font-size: 9pt; font-family: "Segoe UI", sans-serif; color: rgb(33, 33, 33);">U ontvangt niet vaak e-mail van<span class="Apple-converted-space"> </span><a href="mailto:jklymak@uvic.ca" style="color: blue; text-decoration: underline;">jklymak@uvic.ca</a>.<span class="Apple-converted-space"> </span><a href="https://aka.ms/LearnAboutSenderIdentification" style="color: blue; text-decoration: underline;">Ontdek waarom dit belangrijk is</a><o:p></o:p></span></div></div></td><td width="75" style="aspect-ratio: revert !important; background: revert !important; block-size: revert !important; border: revert !important; bottom: revert !important; color: revert !important; color-scheme: revert !important; content-visibility: revert !important; cursor: revert !important; direction: revert !important; display: revert !important; font-size: revert !important; height: revert !important; hyphens: revert !important; letter-spacing: revert !important; line-height: revert !important; margin: revert !important; opacity: revert !important; order: revert !important; outline: revert !important; overflow: revert !important; padding: revert !important; position: revert !important; resize: revert !important; rotate: revert !important; scale: revert !important; tab-size: revert !important; table-layout: revert !important; text-align: revert !important; text-indent: revert !important; text-orientation: revert !important; text-overflow: revert !important; text-shadow: revert !important; text-transform: revert !important; text-wrap: revert !important; top: revert !important; transition: revert !important; vertical-align: revert !important; visibility: revert !important; white-space: revert !important; width: revert !important; word-break: revert !important; word-spacing: revert !important; writing-mode: revert !important; zoom: revert !important;"></td></tr></tbody></table><div><div style="border: 1pt solid rgb(156, 101, 0); padding: 2pt;"><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif; line-height: 12pt; background: rgb(255, 235, 156);"><span style="font-size: 10pt; font-family: Calibri, sans-serif; color: rgb(156, 101, 0);">CAUTION:</span><span style="font-size: 10pt; font-family: Calibri, sans-serif;"><span class="Apple-converted-space"> </span>This email originated from outside of Utrecht University. Do not click links or open attachments unless you recognize the sender and know the content is safe.<o:p></o:p></span></div></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div><div><div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;">Good advice from Martin<o:p></o:p></div></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;">Also check that your forcing fields are the correct order. I once had wind forcing transposed so it varied in x direction instead of y and that definitely gave some strange results. <o:p></o:p></div></div></div><div id="ms-outlook-mobile-body-separator-line"><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div></div><div id="ms-outlook-mobile-signature"><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;">Get<span class="Apple-converted-space"> </span><a href="https://aka.ms/o0ukef" style="color: blue; text-decoration: underline;">Outlook for iOS</a><o:p></o:p></div></div><div class="MsoNormal" align="center" style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif; text-align: center;"><hr size="2" width="98%" align="center"></div><div id="divRplyFwdMsg"><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><b><span style="font-size: 11pt; font-family: Calibri, sans-serif;">From:</span></b><span style="font-size: 11pt; font-family: Calibri, sans-serif;"><span class="Apple-converted-space"> </span>MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" style="color: blue; text-decoration: underline;">mitgcm-support-bounces@mitgcm.org</a>> on behalf of Martin Losch <<a href="mailto:Martin.Losch@awi.de" style="color: blue; text-decoration: underline;">Martin.Losch@awi.de</a>><br><b>Sent:</b><span class="Apple-converted-space"> </span>Tuesday, April 8, 2025 2:57:53 AM<br><b>To:</b><span class="Apple-converted-space"> </span>MITgcm Support <<a href="mailto:mitgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;">mitgcm-support@mitgcm.org</a>><br><b>Subject:</b><span class="Apple-converted-space"> </span>Re: [MITgcm-support] Irregular Velocity Fields Near Basin Boundaries in MITgcm</span><o:p></o:p></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"> <o:p></o:p></div></div></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;">Hi again,<span class="Apple-converted-space"> </span><o:p></o:p></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;">my approach in these kind of problems: Try to remove everything, forcing, intital conditions, etc. and simplify the problem as much as possible (remove convection, advection of tracer, of momentum, etc.) and see if the problem persists (shouldn’t) and then add them back one by one.<o:p></o:p></div></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;">You problem is simple enough that it should not have these issues, so I can only guess that it has to do with forcing or initial fields (unless you modified the code anywhere, what does "git diff” say? what’s in your code directory?)<o:p></o:p></div></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;">Martin<o:p></o:p></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><br><br><o:p></o:p></div><blockquote style="margin-top: 5pt; margin-bottom: 5pt;"><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;">On 7. Apr 2025, at 19:14, Vanderborght, E.Y.P. (Elian) <<a href="mailto:e.y.p.vanderborght@uu.nl" style="color: blue; text-decoration: underline;">e.y.p.vanderborght@uu.nl</a>> wrote:<o:p></o:p></div></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div><div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">Dear Martin,<br><br>Thank you very much for your response and suggestions. I’ve implemented all the proposed changes. Unfortunately, the output still shows irregularities near the northern boundary, even after significantly increasing viscAhGrid. Switching to synchronous time stepping also did not resolve the issue. The striping persists across all three components of the velocity field.<br><br>Any further ideas are appreciated!<br><br>Best regards,<br>Elian<br>-----Oorspronkelijk bericht-----<br>Van: MITgcm-support <</span><a href="mailto:mitgcm-support-bounces@mitgcm.org" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">mitgcm-support-bounces@mitgcm.org</span></a><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">> Namens Martin Losch<br>Verzonden: vrijdag 4 april 2025 8:12<br>Aan: MITgcm Support <</span><a href="mailto:mitgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">mitgcm-support@mitgcm.org</span></a><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">><br>Onderwerp: Re: [MITgcm-support] Irregular Velocity Fields Near Basin Boundaries in MITgcm<br><br>[U ontvangt vaak geen e-mail van<span class="xapple-converted-space"> </span></span><a href="mailto:martin.losch@awi.de" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">martin.losch@awi.de</span></a><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">. Informatie over waarom dit belangrijk is op<span class="xapple-converted-space"> </span></span><a href="https://aka.ms/LearnAboutSenderIdentification%5d" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">https://aka.ms/LearnAboutSenderIdentification]</span></a><span style="font-size: 9pt; font-family: Helvetica, sans-serif;"><br><br>CAUTION: This email originated from outside of Utrecht University. Do not click links or open attachments unless you recognize the sender and know the content is safe.<br><br><br>Hi Elian,<br><br>just a hunch, but your domain stretches from 3degN to 67degN, where the zonal (dx) spacing will have decreased from 1*cos(3deg) ~100km to 1*cos(64deg)~50km, the dy spacing will always be ~111km, so you’ll have non-square grid cells. At 1deg resolution it shouldn’t matter that much, but this could be a problem for the cfl-number, especially the one for viscosity. I’d try to use viscAhGrid instead of viscAh (see mom_common/mom_calc_visc.F or the corresponding part in the documentation), instead of viscAh=2e4, you could use viscAhGrid=0.0065, if I got it right. The vertical velocity is just a consequence of the noise in the horizontal velocities (as it is computed from the divergence of the horizontal velocities)<br><br>I would turn off the CD-scheme (this may increase the noise, but a smoothing method may mask problems that you want to see for debugging). In general, I usually use viscA4Grid ~ 0.01 instead.<br><br>Did you try synchronous time stepping (i.e. only use deltaT, but not deltaTtracer/Clock, and what’s missing is deltaTFreeSurf)?<br><br>Is the striping only in zonal direction?<br><br>unrelated: globalFiles does not always work, is not necessary here, useSingleCpuIO is enough.<br><br>Martin<br><br style="font-variant-caps: normal; text-align: start; word-spacing: 0px;"><br></span><o:p></o:p></div><blockquote style="margin-top: 5pt; margin-bottom: 5pt;"><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">On 3. Apr 2025, at 21:18, Vanderborght, E.Y.P. (Elian) <<a href="mailto:e.y.p.vanderborght@uu.nl" style="color: blue; text-decoration: underline;">e.y.p.vanderborght@uu.nl</a>> wrote:<br><br>Hi Jody,<br>Thanks for the reply. The tracer field look fine, I also checked the convection adjustment index and there are no irregularities found here.<br>Best regards,<br>Elian<br>Van: MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" style="color: blue; text-decoration: underline;">mitgcm-support-bounces@mitgcm.org</a>> Namens Jody<br>Klymak<br>Verzonden: donderdag 3 april 2025 20:41<br>Aan:<span class="Apple-converted-space"> </span><a href="mailto:mitgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;">mitgcm-support@mitgcm.org</a><br>Onderwerp: Re: [MITgcm-support] Irregular Velocity Fields Near Basin<br>Boundaries in MITgcm U ontvangt niet vaak e-mail van<span class="Apple-converted-space"> </span><a href="mailto:jklymak@uvic.ca" style="color: blue; text-decoration: underline;">jklymak@uvic.ca</a>.<br>Ontdek waarom dit belangrijk is<br>CAUTION: This email originated from outside of Utrecht University. Do not click links or open attachments unless you recognize the sender and know the content is safe.<br>Hi Elian,<br>I would look at your tracer fields. Is that maybe just (poorly resolved) convection?<br>Cheers, Jody<br>From: MITgcm-support <<a href="mailto:mitgcm-support-bounces@mitgcm.org" style="color: blue; text-decoration: underline;">mitgcm-support-bounces@mitgcm.org</a>> on behalf of<br>"Vanderborght, E.Y.P. (Elian)" <<a href="mailto:e.y.p.vanderborght@uu.nl" style="color: blue; text-decoration: underline;">e.y.p.vanderborght@uu.nl</a>><br>Reply-To: "<a href="mailto:mitgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;">mitgcm-support@mitgcm.org</a>" <<a href="mailto:MITgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;">MITgcm-support@mitgcm.org</a>><br>Date: Thursday, April 3, 2025 at 8:25 AM<br>To: MITgcm Support <<a href="mailto:MITgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;">MITgcm-support@mitgcm.org</a>><br>Subject: [MITgcm-support] Irregular Velocity Fields Near Basin<br>Boundaries in MITgcm You don't often get email from<br><a href="mailto:e.y.p.vanderborght@uu.nl" style="color: blue; text-decoration: underline;">e.y.p.vanderborght@uu.nl</a>. Learn why this is important Dear MITgcm<br>community, I am running MITgcm in a rectangular basin with a single-hemisphere setup. However, I am encountering irregularities in both vertical and horizontal velocity fields near the basin boundaries (see attached figure).<br>Initially, I suspected insufficient horizontal viscosity. While increasing it to unrealistically high values reduced the issue, it did not resolve it entirely. I have also tested the following:<br><br> • Decreasing time-step,<br> • Using different advection schemes,<br> • Increasing horizontal diffusion,<br> • Increasing vertical viscosity.<br>None of these approaches have fully addressed the problem. I have also reviewed studies using similar setups, but they do not report encountering this issue. Does anyone have insights into what might be causing this or suggestions for a fix? I have attached my data file and SIZE.h for reference.<br>Any help would be greatly appreciated.<br>Best regards,<br>Elian<br>_______________________________________________<br>MITgcm-support mailing list<br><a href="mailto:MITgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;">MITgcm-support@mitgcm.org</a><br><a href="http://mailm/" style="color: blue; text-decoration: underline;">http://mailm/</a><br><a href="http://an.mitgcm.org/" style="color: blue; text-decoration: underline;">an.mitgcm.org</a>%2Fmailman%2Flistinfo%2Fmitgcm-support&data=05%7C02%7Ce.y<br>.p.vanderborght%<a href="http://40uu.nl/" style="color: blue; text-decoration: underline;">40uu.nl</a>%7Ce0d45c7644744079e0be08dd733faf50%7Cd72758a0a<br>4464e0fa0aa4bf95a4a10e7%7C0%7C0%7C638793439556237863%7CUnknown%7CTWFpb<br>GZsb3d8eyJFbXB0eU1hcGkiOnRydWUsIlYiOiIwLjAuMDAwMCIsIlAiOiJXaW4zMiIsIkF<br>OIjoiTWFpbCIsIldUIjoyfQ%3D%3D%7C0%7C%7C%7C&sdata=AL4JLgPcbLRGF54MU5Z6m<br>zdGh9OjhJQ0MkmzJkocVd4%3D&reserved=0<o:p></o:p></span></div></blockquote><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;"><br><br>_______________________________________________<br>MITgcm-support mailing list<br></span><a href="mailto:MITgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">MITgcm-support@mitgcm.org</span></a><span style="font-size: 9pt; font-family: Helvetica, sans-serif;"><br></span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</span></a><span style="font-size: 9pt; font-family: Helvetica, sans-serif;"><br>_______________________________________________<br>MITgcm-support mailing list<br></span><a href="mailto:MITgcm-support@mitgcm.org" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">MITgcm-support@mitgcm.org</span></a><span style="font-size: 9pt; font-family: Helvetica, sans-serif;"><br></span><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" style="color: blue; text-decoration: underline;"><span style="font-size: 9pt; font-family: Helvetica, sans-serif;">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</span></a><o:p></o:p></div></div></blockquote></div><div style="margin: 0cm; font-size: 12pt; font-family: Aptos, sans-serif;"><o:p> </o:p></div></div></div></div></div></div><span id="cid:EF507C02-8CD3-4808-A2CA-908E6E1CB221"><UVEL.png></span><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">_______________________________________________</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><span style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none; float: none; display: inline !important;">MITgcm-support mailing list</span><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><a href="mailto:MITgcm-support@mitgcm.org" style="color: blue; text-decoration: underline; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;">MITgcm-support@mitgcm.org</a><br style="caret-color: rgb(0, 0, 0); font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; text-decoration: none;"><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" style="color: blue; text-decoration: underline; font-family: Helvetica; font-size: 12px; font-style: normal; font-variant-caps: normal; font-weight: 400; letter-spacing: normal; orphans: auto; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; widows: auto; word-spacing: 0px; -webkit-text-stroke-width: 0px;">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a></div></blockquote></div><br></div></body></html>