<div dir="ltr"><div>Some more info on this ..</div><div><br></div><div>I've written snapshots a few timesteps before the crash and discovered that sea ice velocities explode at point i,j= (235,355) (see netcdf attached) but instabilities in seaice velocities start to develop since iter 320486.</div><div><br></div><div>Enrico<br></div><div><br></div><div></div><div class="gmail_chip gmail_drive_chip" style="width:396px;height:18px;max-height:18px;background-color:#f5f5f5;padding:5px;color:#222;font-family:arial;font-style:normal;font-weight:bold;font-size:13px;border:1px solid #ddd;line-height:1"><a href="https://drive.google.com/file/d/1R7k3xtN_zPEdcrjcUfafiQ3WjNJ2lEwx/view?usp=drive_web" target="_blank" style="display:inline-block;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0px;border:medium none;width:100%"><img style="vertical-align: bottom; border: none;" src="https://ssl.gstatic.com/docs/doclist/images/icon_10_generic_list.png"> <span dir="ltr" style="color:#15c;text-decoration:none;vertical-align:bottom">snapshots_last.nc</span></a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno gio 7 gen 2021 alle ore 15:50 Pochini, Enrico <<a href="mailto:epochini@inogs.it" target="_blank">epochini@inogs.it</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 dir="ltr"><div dir="ltr">The attachment didn't get through <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/file/d/1brDxfDiu35GlcTrkiqRYmpPT2LL0urqx/view?usp=drive_web" style="display:inline-block;overflow:hidden;text-overflow:ellipsis;white-space:nowrap;text-decoration:none;padding:1px 0px;border:medium none;width:100%" 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">STDOUT.log</span></a></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno gio 7 gen 2021 alle ore 15:47 Pochini, Enrico <<a href="mailto:epochini@inogs.it" target="_blank">epochini@inogs.it</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 dir="ltr"><div>This is the log of the longer run.. 1900 is a fake initial date, the real time is 21'000 years ago <br></div><div><br></div><div>Enrico<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">Il giorno gio 7 gen 2021 alle ore 15:33 Pochini, Enrico <<a href="mailto:epochini@inogs.it" target="_blank">epochini@inogs.it</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 dir="ltr"><div>Dear community,</div><div><br></div><div>I am trying to run a regional model with SHELFICE and SEAICE packages. <br></div><div>However I'm still having problems with  crashes while employing the EVP solver. <br></div><div><br></div><div>The error I get is the general <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>ERROR STOP MOM_IMPLICIT_R: error when solving 3-Diag problem.</div></blockquote><div>and I only get it with the EVP; with the LSR solver the model runs fine for more years but the ice I get seems to be less realistic. <br></div><div><br></div><div>I was able to run the model for just two years with EVP without crashes with the namelist files attached.<br></div><div><br></div><div>To restart the model I have tried playing with OBCS_SEAICE options, sponge layer thickness and relaxation times and EVP alpha/beta parameters (from 10 to 600). I tried also disabling OB for seaice. I have 5km resolution, 300 seconds deltaT, vertical levels with min 10m thickness.<br></div><div><br></div><div></div><div>Can you point me to some tests to fix this problem? <br></div><div>How can I print debugging information in the STDOUT? In the SEAICE_OPTIONS.h I set #define SEAICE_DEBUG and have debugLevel>0 in "data" but it doesn't seem to work with EVP. <br></div><div></div><div><br></div><div>Kind regards,</div><div><br></div><div>Enrico<br></div><div><br></div></div>
</blockquote></div>
</blockquote></div></div>
</blockquote></div>