<div dir="ltr"><div>Hi Gus,</div><div><br></div><div>As per your suggestion, I made the respective changes and tried to run the executable again. Below is the log generated</div><div><br></div><div>$ mpirun -np 60 ./mitgcmuv <br>forrtl: error (72): floating overflow<br>Image              PC                                    Routine                     Line          Source             <br>libifcoremt.so.5   00002AD443246555         for__signal_handl     Unknown     Unknown<br>libpthread-2.17.s  00002AD442DB35F0       Unknown               Unknown       Unknown<br>libnetcdf.so.15.2  00002AD44121C4B3       __libm_exp_e7         Unknown     Unknown<br>mitgcmuv           0000000000AC0FF7        exf_bulkformulae_        4598         exf_bulkformulae.f<br>mitgcmuv           0000000000B02334         exf_getforcing_          4430           exf_getforcing.f<br>mitgcmuv           000000000128726E         load_fields_drive        2141           load_fields_driver.f<br>mitgcmuv           0000000000C45A25        forward_step_            2340           forward_step.f<br>mitgcmuv           0000000001290200         main_do_loop_            2078         main_do_loop.f<br>mitgcmuv           0000000001C283F6         the_main_loop_           2097        the_main_loop.f<br>mitgcmuv           0000000001C28955         the_model_main_          2421      the_model_main.f<br>mitgcmuv           0000000001290615         MAIN__                   4286            main.f<br>mitgcmuv           0000000000403412         Unknown                Unknown      Unknown<br><a href="http://libc-2.17.so" target="_blank">libc-2.17.so</a>       00002AD445C2A505        __libc_start_main    Unknown      Unknown<br>mitgcmuv           0000000000403319         Unknown                Unknown      Unknown</div><div>----------------------------------------------------------------------------------------------------------------------------------------------</div><div><br></div><div>The STDERR file reads:</div><div>(PID.TID 0030.0001) ** WARNING ** MNC_READPARMS: incomplete MNC pickup files implementation<br>(PID.TID 0030.0001) ** WARNING ** MNC_READPARMS: => pickup_write_mnc=T not recommanded<br>(PID.TID 0030.0001) ** WARNING ** MNC_READPARMS: => pickup_read_mnc=T not working for some set-up<br>(PID.TID 0030.0001) ** WARNING ** INI_MODEL_IO: globalFiles=TRUE is not safe in Multi-processors (MPI) run<br>(PID.TID 0030.0001) ** WARNING ** INI_MODEL_IO: use instead "useSingleCpuIO=.TRUE."<br>(PID.TID 0030.0001) ** WARNING ** INI_MODEL_IO: use tiled-files to write sections (for OBCS)<br>(PID.TID 0030.0001) ** WARNING ** EXF_CHECK: wind-stress position irrelevant</div><div><br></div><div>Attaching data, data.obcs, data.exf for your reference. I have set deltaTmom=120.0,</div><div>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><br></div><div>Regards</div><div>Kunal<br></div></div><br><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">gus@ldeo.columbia.edu</a>> wrote:<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>Hi Kunal</div><div><br></div><div>To try to nail down where, when, why it fails you could compile in debugging mode,</div><div>ie. start fresh ('make CLEAN' in the build directory or just wipe that directory off) <br></div><div>and run gemake2 with the -devel flag (keep the other flags).</div><div>Then, to increase verbosity add:</div><div>debugLevel = 4,</div><div>to the "data" namelist &PARM01,</div><div>and increase the <br></div><div>monitorFreq <br></div><div>in &PARM03</div><div>to one or a few time steps.</div><div>The STDOUT.XXXX, and STDERR.XXXX files <br></div><div>may give a hint of what is going on (when, where, wny it fails).</div><div><br></div><div>I hope this helps,</div><div>Gus Correa<br></div></div><br><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">kunal.madkaiker02@gmail.com</a>> wrote:<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 All,</div><div><br></div><div>I am trying to simulate U,V currents circulation along the West Coast of India.</div><div>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><br></div><div>But model blows up at the initial stage and I get the error: <br></div><div>Note: The following floating-point exceptions are signalling: IEEE_INVALID_FLAG IEEE_DIVIDE_BY_ZERO IEEE_UNDERFLOW_FLAG<br>STOP MOM_IMPLICIT_R: error when solving 3-Diag problem.</div><div><br></div><div>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></div><div>I have defined the vertical levels keeping delZ(k+1)/delZ(k) < 1.4 ratio in mind. But the issue persists. Kindly advise. <br></div><div>Please let me know if any additional information is required from my side.</div><div><br></div><div>Regards</div><div>Kunal<br></div><div><br></div></div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>
_______________________________________________<br>
MITgcm-support mailing list<br>
<a href="mailto:MITgcm-support@mitgcm.org" target="_blank">MITgcm-support@mitgcm.org</a><br>
<a href="http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support" rel="noreferrer" target="_blank">http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support</a><br>
</blockquote></div>