<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=Windows-1252">
<meta name="Generator" content="Microsoft Exchange Server">
<!-- converted from text --><style><!-- .EmailQuote { margin-left: 1pt; padding-left: 4pt; border-left: #800000 2px solid; } --></style>
</head>
<body>
<meta content="text/html; charset=UTF-8">
<style type="text/css" style="">
<!--
p
{margin-top:0;
margin-bottom:0}
-->
</style>
<div dir="ltr">
<div id="x_divtagdefaultwrapper" dir="ltr" style="font-size:12pt; color:#000000; font-family:Calibri,Arial,Helvetica,sans-serif">
<p>Hi Gus,</p>
<p><br>
</p>
<p>Thanks for your reply.</p>
<p><br>
</p>
<p><font face="Calibri,Arial,Helvetica,sans-serif" style="font-family:Calibri,Arial,Helvetica,sans-serif">I added a flag in my </font><font face="Calibri,Arial,Helvetica,sans-serif" style="font-family:Calibri,Arial,Helvetica,sans-serif">optfile " </font><font face="Calibri,Arial,Helvetica,sans-serif" style="font-family:Calibri,Arial,Helvetica,sans-serif">-ffpe-trap=invalid,zero,overflow</font><font face="Calibri,Arial,Helvetica,sans-serif" size="3">", will
this solve my problem too? how does this compare with your solution?</font><br>
</p>
<p><font face="Calibri,Arial,Helvetica,sans-serif" size="3"><br>
</font></p>
<p><font face="Calibri,Arial,Helvetica,sans-serif" size="3">Yangxin</font></p>
</div>
<hr tabindex="-1" style="display:inline-block; width:98%">
<div id="x_divRplyFwdMsg" dir="ltr"><font face="Calibri, sans-serif" color="#000000" style="font-size:11pt"><b>From:</b> MITgcm-support <mitgcm-support-bounces@mitgcm.org> on behalf of Gus Correa <gus@ldeo.columbia.edu><br>
<b>Sent:</b> Thursday, November 16, 2017 12:14:38 PM<br>
<b>To:</b> MITgcm Support<br>
<b>Subject:</b> Re: [MITgcm-support] NaN - stop running?</font>
<div> </div>
</div>
</div>
<font size="2"><span style="font-size:10pt;">
<div class="PlainText">Hi Yangxin<br>
<br>
Pass the --ieee and --devel flags to genmake2,<br>
when you build the model (do a 'make CLEAN' beforehand).<br>
Something like this (assuming you use MPI):<br>
<br>
/path/to/MITgcm/tools/genmake2 \<br>
--rootdir=/path/to/MITgcm \<br>
--optfile= /path/to/MITgcm/tools/build_options/your_build_options_file \<br>
--mods=/path/to/your/experiment/code \<br>
--mpi \<br>
--ieee \<br>
--devel \<br>
<br>
<br>
If your build_options file was properly written, it will<br>
turn on your compiler's floating point exceptions (as Jody suggested),<br>
and the model will crash at the first NaN.<br>
<br>
Besides, it is also useful to turn on "traceback"<br>
(or "backtrace", check the compiler man page) compiler flag,<br>
which will tell you in which preprocessed ".f" file<br>
(in your build directory)<br>
and on which line of code line the problem happened.<br>
BTW, keep the ".f" files, don't do a 'make clean' while you're<br>
chasing the error.<br>
<br>
If your build_options file doesn't have these features,<br>
you can find some example files that do have that<br>
in the MITgcm original build_options directory.<br>
Look there for "macros" such as IEEE, and DEVEL.<br>
Copy the corresponding blocks to your build_options file,<br>
and adjust to your compiler's specific flags.<br>
<br>
I hope this helps,<br>
Gus Correa<br>
<br>
On 11/15/2017 01:33 PM, Jody Klymak wrote:<br>
> You want to look up floating-point exception handling for your compiler.<br>
> <br>
> Cheers, Jody<br>
> <br>
> <br>
> <br>
>> On 15 Nov 2017, at 10:18 AM, Yangxin He <y67he@uwaterloo.ca <br>
>> <<a href="mailto:y67he@uwaterloo.ca">mailto:y67he@uwaterloo.ca</a>>> wrote:<br>
>><br>
>> Hey,<br>
>><br>
>> I am having a (probably) small problem. Sometimes my output data gives <br>
>> all NaN for various reasons. But it still kept running until the end. <br>
>> Is there any way that I can make the model stop running once it starts <br>
>> to generate NaNs?<br>
>><br>
>> Thanks<br>
>><br>
>> Yangxin<br>
>> _______________________________________________<br>
>> MITgcm-support mailing list<br>
>> MITgcm-support@mitgcm.org <<a href="mailto:MITgcm-support@mitgcm.org">mailto:MITgcm-support@mitgcm.org</a>><br>
>> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
> <br>
> --<br>
> Jody Klymak<br>
> <a href="http://web.uvic.ca/~jklymak/">http://web.uvic.ca/~jklymak/</a><br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> <br>
> _______________________________________________<br>
> MITgcm-support mailing list<br>
> MITgcm-support@mitgcm.org<br>
> <a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
> <br>
<br>
_______________________________________________<br>
MITgcm-support mailing list<br>
MITgcm-support@mitgcm.org<br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</div>
</span></font>
</body>
</html>