<div dir="ltr"><div>Hi Jody,<br></div><div><br></div><div>Your reply ended up in spam, I noticed only now.. <br></div><div>Eventually I managed to run it with a 5 s timestep. It was indeed a stability issue; I have switched to a cartesian grid and things are working now with a reasonable (150 s) time step.</div><div><br></div><div>Thank you and best regards,</div><div><br></div><div>Enrico<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno mar 22 ott 2019 alle ore 18:05 Jody Klymak <<a href="mailto:jklymak@uvic.ca">jklymak@uvic.ca</a>> ha scritto:<br></div><blockquote class="gmail_quote" style="margin:0px 0px 0px 0.8ex;border-left:1px solid rgb(204,204,204);padding-left:1ex"><div style="overflow-wrap: break-word;">Hi Enrico,<div><br></div><div>You can check the cfl criteria by looking at STDOUT.0000 (advu_cfl, advw_cfl etc). You don’t say anything about your vertical resolution, but if it is fine enough to use non-hydrostatic then you are quite possibly violating the CFL criteria in the vertical velocity.  </div><div><br></div><div>For testing purposes you may want to turn nonHydrostatic off (and maybe consider whether 0.5 degree runs can properly represent non-hydrostatic effects).  That may not help your cfl criteria, but the model will run faster.  </div><div><br></div><div>Cheers,   Jody<br><div><br><blockquote type="cite"><div>On 22 Oct 2019, at 06:53, Pochini, Enrico <<a href="mailto:epochini@inogs.it" target="_blank">epochini@inogs.it</a>> wrote:</div><br><div><div dir="ltr"><div>Dear MITgcm support group,</div><div><br></div><div>I'm trying to develop a model of the Ross Sea for my PhD. Until now I've generated the grid and I've started testing it with constant atmospheric forcing, boundary conditions, and initial conditions, but I'm crashing a few kilo-seconds after initialization.</div><div><br></div><div>To generate it I followed what is written in this page:<br></div><div><a href="https://github.com/knaughten/mitgcm_python/wiki/Creating-a-new-MITgcm-domain#grid" target="_blank">https://github.com/knaughten/mitgcm_python/wiki/Creating-a-new-MITgcm-domain#grid</a></div><div></div><div><br></div><div>However even minimal runs with disabled SHELFice, OBCS packages, imposing no surface forcing, with uniform initial conditions, still it develops instabilities and crashes.<br></div><div><br></div><div>I'm absolutely new to MITgcm and ocean modeling and I would need your help to understand what is causing the issues, in particular if the bathymetry has some pathologies or the (spherical) grid violates CFL somewhere .. (?) I'm using 150 seconds timestep, 0.5 degree grid.<br></div><div><br></div><div>I've attached four minimal test outputs,  with and without ice shelf, with and without smoothing of ice shelf draft and bathymetry. The number with the file is the last iteration printed. <br></div><div>Thank you.<br></div><div><br></div><div>Best regards,<div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:rgb(245,245,245);padding:5px;color:rgb(34,34,34);font-family:arial;font-style:normal;font-weight:bold;font-size:13px;border:1px solid rgb(221,221,221);line-height:1"><a href="https://drive.google.com/a/inogs.it/file/d/1SM5mtzq3WgF61EbLh4404N5cz3czIuSI/view?usp=drive_web" style="display:inline-block;max-width:366px;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0px;border:medium none" target="_blank"><img style="vertical-align: bottom; border: medium none;" src="https://ssl.gstatic.com/docs/doclist/images/icon_10_generic_list.png"> <span dir="ltr" style="color:rgb(17,85,204);text-decoration:none;vertical-align:bottom">file_support.tar.gz</span></a><img style="opacity: 0.55; float: right; display: none;"></div></div><div><br></div><div>Enrico Pochini<br></div><div><br></div><div>PhD student, University of Trieste, Italy<br></div></div>
_______________________________________________<br>MITgcm-support mailing list<br><a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br><a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br></div></blockquote></div><br><div>
<span style="border-collapse:separate;border-spacing:0px;color:rgb(0,0,0);font-family:"Lucida Sans Typewriter";font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div style="overflow-wrap: break-word;"><span style="border-collapse:separate;border-spacing:0px;color:rgb(0,0,0);font-family:"Lucida Sans Typewriter";font-size:12px;font-style:normal;font-variant:normal;font-weight:normal;letter-spacing:normal;line-height:normal;text-indent:0px;text-transform:none;white-space:normal;word-spacing:0px"><div>--</div><div>Jody Klymak    </div><div><a href="http://web.uvic.ca/~jklymak/" target="_blank">http://web.uvic.ca/~jklymak/</a></div><div><br></div><div><br></div><br></span></div></span><br>
</div>
<br></div></div>_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>