<div dir="ltr"><div dir="ltr"><font face="arial, helvetica, sans-serif" color="#000000">Hi, </font><div><font face="arial, helvetica, sans-serif" color="#000000">I'm running a lake simulation with the nonlinear free surface and the r* coordinate, an eventual aim is to include seiche behaviour in the model. The model is breaking consistently as I mentioned in an email before, and this is the error we have been getting:</font></div><div><font face="arial, helvetica, sans-serif" color="#000000"><br></font></div><div><div><p style="margin-top:0px;margin-bottom:0px"><font face="arial, helvetica, sans-serif" color="#000000"> fail at i,j=  14  11 ; rStarFacC,H,eta =  0.199897  1.939551E-01 -1.551841E-01</font></p><p style="margin-top:0px;margin-bottom:0px"><font face="arial, helvetica, sans-serif" color="#000000">WARNING: r*FacC < hFacInf at       1 pts : bi,bj,Thid,Iter=   1   1   1     63844</font></p><p style="margin-top:0px;margin-bottom:0px"><font face="arial, helvetica, sans-serif" color="#000000">STOP in CALC_R_STAR : too SMALL rStarFac[C,W,S] !</font></p><p style="margin-top:0px;margin-bottom:0px"><font face="arial, helvetica, sans-serif" color="#000000"><br></font></p><p style="margin-top:0px;margin-bottom:0px"><font face="arial, helvetica, sans-serif" color="#000000">My understanding is that this is due to near-drainage of the cell. What can we do to solve this problem? Should we change the bathymetry to have a minimum depth, or change the grid spacing? Is it possible with the MITgcm to simulate large-amplitude seiches or is that pointless with the current model? </font></p><p style="margin-top:0px;margin-bottom:0px"><font face="arial, helvetica, sans-serif" color="#000000"><br></font></p><p style="margin-top:0px;margin-bottom:0px"><span style="color:rgb(0,0,0);font-family:arial,helvetica,sans-serif">Regards, </span><br></p></div><div dir="ltr" class="gmail_signature">Senja</div></div></div></div>