<div dir="ltr"><div>Dear community,</div><div><br></div><div>I'm facing the problem of Nans generation not stopping model execution. I've tried adding the -fpe0 option (on ifortran) which does the job; however it signals a <br></div><div><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>forrtl: error (73): floating divide by zero<br></div></blockquote><div><br></div><div> from kpp_routines.f right at zero timestep (after monitor printing to STDOUT).<br></div><div><br></div><div>The incriminated lines are in FIND_BETA subroutine:<br></div><div><br></div><div>               DBLOC(I,J,K-1) = gravity * (RHOK(I,J) - RHOKM1(I,J))/<br>     &           (RHOK(I,J) + rhoConst )<br>               DBSFC(I,J,K)   = gravity * (RHOK(I,J) - RHO1K (I,J)) /<br>     &           (RHOK(I,J) + rhoConst )</div><div><br></div><div>adding a small constant to the denominators solved the signalling issue. <br></div><div>But is there a better way to solve it? Are those variables RHOK rhoConst initialized correctly? I have rhoConst set in data file.</div><div><br></div><div>Following <a href="http://mailman.mitgcm.org/pipermail/mitgcm-support/2006-April/003931.html">http://mailman.mitgcm.org/pipermail/mitgcm-support/2006-April/003931.html</a></div><div>I also tried the -Ktrap=fp option, that didn't stop execution on NaN, and -fpeX X>0 with same result.<br></div><div><br></div><div><br></div><div>P.S. Unrelated to the present mail, but may be of interest to other users since I didn't find a solution on the forum: on some cluster with jobs queue handling, jobs are not killed on errors like MOM_IMPLICIT_R and similar; the model continues "not running" but consuming CPU-hours. The solution to this is to substitute STOP statements with ERROR STOP in error calls.<br></div><div><br></div><div><div>Cheers,</div><div><br></div><div>Enrico</div></div></div>