<html><head><meta http-equiv="Content-Type" content="text/html; charset=us-ascii"></head><body style="word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class="">Hi Kunal.<div class=""><br class=""></div><div class="">Check STDOUT.0000 as that is more relevant.  Does the cfl criteria blow up?  If so, you will simply have to reduce the time step. Its possible you have also set the model up incorrectly and it is convecting initially, which will easily violate the cfl criteria.  It is also hard to tell fro what you are giving us when the model blows up.  time step 1?  time step 100? Finally, its often useful to plot the fields to see where the instability is happening.  That may require you to save quite a bit of data, but its hard to debug in the absence of information. </div><div class=""><br class=""></div><div class="">Best of luck!  Jody<br class=""><div><br class=""><blockquote type="cite" class=""><div class="">On 13 Sep 2020, at 07:22, kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com" class="">kunal.madkaiker02@gmail.com</a>> wrote:</div><br class="Apple-interchange-newline"><div class=""><div dir="ltr" class=""><div class="">Hi Gus,</div><div class=""><br class=""></div><div class="">As per your suggestion, I made the respective changes and tried to run the executable again. Below is the log generated</div><div class=""><br class=""></div><div class="">$ mpirun -np 60 ./mitgcmuv <br class="">forrtl: error (72): floating overflow<br class="">Image              PC                                    Routine                     Line          Source             <br class="">libifcoremt.so.5   00002AD443246555         for__signal_handl     Unknown     Unknown<br class="">libpthread-2.17.s  00002AD442DB35F0       Unknown               Unknown       Unknown<br class="">libnetcdf.so.15.2  00002AD44121C4B3       __libm_exp_e7         Unknown     Unknown<br class="">mitgcmuv           0000000000AC0FF7        exf_bulkformulae_        4598         exf_bulkformulae.f<br class="">mitgcmuv           0000000000B02334         exf_getforcing_          4430           exf_getforcing.f<br class="">mitgcmuv           000000000128726E         load_fields_drive        2141           load_fields_driver.f<br class="">mitgcmuv           0000000000C45A25        forward_step_            2340           forward_step.f<br class="">mitgcmuv           0000000001290200         main_do_loop_            2078         main_do_loop.f<br class="">mitgcmuv           0000000001C283F6         the_main_loop_           2097        the_main_loop.f<br class="">mitgcmuv           0000000001C28955         the_model_main_          2421      the_model_main.f<br class="">mitgcmuv           0000000001290615         MAIN__                   4286            main.f<br class="">mitgcmuv           0000000000403412         Unknown                Unknown      Unknown<br class=""><a href="http://libc-2.17.so/" target="_blank" class="">libc-2.17.so</a>       00002AD445C2A505        __libc_start_main    Unknown      Unknown<br class="">mitgcmuv           0000000000403319         Unknown                Unknown      Unknown</div><div class="">----------------------------------------------------------------------------------------------------------------------------------------------</div><div class=""><br class=""></div><div class="">The STDERR file reads:</div><div class="">(PID.TID 0030.0001) ** WARNING ** MNC_READPARMS: incomplete MNC pickup files implementation<br class="">(PID.TID 0030.0001) ** WARNING ** MNC_READPARMS: => pickup_write_mnc=T not recommanded<br class="">(PID.TID 0030.0001) ** WARNING ** MNC_READPARMS: => pickup_read_mnc=T not working for some set-up<br class="">(PID.TID 0030.0001) ** WARNING ** INI_MODEL_IO: globalFiles=TRUE is not safe in Multi-processors (MPI) run<br class="">(PID.TID 0030.0001) ** WARNING ** INI_MODEL_IO: use instead "useSingleCpuIO=.TRUE."<br class="">(PID.TID 0030.0001) ** WARNING ** INI_MODEL_IO: use tiled-files to write sections (for OBCS)<br class="">(PID.TID 0030.0001) ** WARNING ** EXF_CHECK: wind-stress position irrelevant</div><div class=""><br class=""></div><div class="">Attaching data, data.obcs, data.exf for your reference. I have set deltaTmom=120.0,</div><div class="">What I am understanding is that the model is blowing up due to overestimation of few values and not because of any error. Am I right?</div><div class=""><br class=""></div><div class="">Regards</div><div class="">Kunal<br class=""></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sun, Sep 13, 2020 at 6:47 AM Gus Correa <<a href="mailto:gus@ldeo.columbia.edu" target="_blank" class="">gus@ldeo.columbia.edu</a>> wrote:<br class=""></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" class=""><div class="">Hi Kunal</div><div class=""><br class=""></div><div class="">To try to nail down where, when, why it fails you could compile in debugging mode,</div><div class="">ie. start fresh ('make CLEAN' in the build directory or just wipe that directory off) <br class=""></div><div class="">and run gemake2 with the -devel flag (keep the other flags).</div><div class="">Then, to increase verbosity add:</div><div class="">debugLevel = 4,</div><div class="">to the "data" namelist &PARM01,</div><div class="">and increase the <br class=""></div><div class="">monitorFreq <br class=""></div><div class="">in &PARM03</div><div class="">to one or a few time steps.</div><div class="">The STDOUT.XXXX, and STDERR.XXXX files <br class=""></div><div class="">may give a hint of what is going on (when, where, wny it fails).</div><div class=""><br class=""></div><div class="">I hope this helps,</div><div class="">Gus Correa<br class=""></div></div><br class=""><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Sat, Sep 12, 2020 at 7:38 PM kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com" target="_blank" class="">kunal.madkaiker02@gmail.com</a>> wrote:<br class=""></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" class=""><div class="">Dear All,</div><div class=""><br class=""></div><div class="">I am trying to simulate U,V currents circulation along the West Coast of India.</div><div class="">I have a grid of 720 x 1560 with a high resolution of 1.45km x 1.45km, with 25 levels in the vertical from 0 to 2150m. I have set hFacMin=0.3 and hFacMinDz=10</div><div class=""><br class=""></div><div class="">But model blows up at the initial stage and I get the error: <br class=""></div><div class="">Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG<br class="">STOP MOM_IMPLICIT_R: error when solving 3-Diag problem.</div><div class=""><br class=""></div><div class="">I have tried changing viscAh from 1 to 1000 m2/s and viscAz from 0.02 to 0.001 m2/s. Also tried with viscAhgrid=0.1.</div><div class=""></div><div class="">I have defined the vertical levels keeping delZ(k+1)/delZ(k) < 1.4 ratio in mind. But the issue persists. Kindly advise. <br class=""></div><div class="">Please let me know if any additional information is required from my side.</div><div class=""><br class=""></div><div class="">Regards</div><div class="">Kunal<br class=""></div><div class=""><br class=""></div></div>
_______________________________________________<br class="">
MITgcm-support mailing list<br class="">
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank" class="">MITgcm-support@mitgcm.org</a><br class="">
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">
</blockquote></div>
_______________________________________________<br class="">
MITgcm-support mailing list<br class="">
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank" class="">MITgcm-support@mitgcm.org</a><br class="">
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank" class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br class="">
</blockquote></div>
<span id="cid:f_kf16sty10"><data></span><span id="cid:f_kf16stys2"><data.obcs></span><span id="cid:f_kf16styo1"><data.exf></span><span id="cid:f_kf16styv3"><STDOUT.0025></span>_______________________________________________<br class="">MITgcm-support mailing list<br class=""><a href="mailto:MITgcm-support@mitgcm.org" class="">MITgcm-support@mitgcm.org</a><br class="">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support<br class=""></div></blockquote></div><br class=""><div class="">
<div style="color: rgb(0, 0, 0); letter-spacing: normal; text-align: start; text-indent: 0px; text-transform: none; white-space: normal; word-spacing: 0px; -webkit-text-stroke-width: 0px; word-wrap: break-word; -webkit-nbsp-mode: space; line-break: after-white-space;" class=""><span class="Apple-style-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-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px;"><div style="word-wrap: break-word;" class=""><span class="Apple-style-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-ligatures: normal; font-variant-caps: normal; font-variant-east-asian: normal; font-variant-position: normal; font-weight: normal; letter-spacing: normal; line-height: normal; text-indent: 0px; text-transform: none; orphans: 2; white-space: normal; widows: 2; word-spacing: 0px;"><div class="">--</div><div class="">Jody Klymak    </div><div class=""><a href="http://ocean-physics.seos.uvic.ca/~jklymak/" class="">http://ocean-physics.seos.uvic.ca/~jklymak/</a></div><div class=""><br class="khtml-block-placeholder"></div><div class=""><br class="khtml-block-placeholder"></div><br class="Apple-interchange-newline"></span></div></span></div><br class="Apple-interchange-newline"><br class="Apple-interchange-newline">
</div>
<br class=""></div></body></html>