<div dir="ltr"><div>Hello Martin and Kaveh,</div><div><br></div><div>As mentioned by both of you, I had NaN values in my input files (over land region), <br></div><div>which I replaced with zero. Now the model is running successfully.</div><div><br></div><div>Thanks to you both.</div><div><br></div><div>Regards</div><div>Kunal<br></div></div><br><div class="gmail_quote"><div dir="ltr" class="gmail_attr">On Mon, Feb 24, 2020 at 2:46 PM Martin Losch <<a href="mailto:Martin.Losch@awi.de">Martin.Losch@awi.de</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">Hi Kunal,<br>
<br>
this iteration solves the hydrostatic pressure equation (dp = - g rho(z) dz) and usually converges very quickly.<br>
<br>
I think you need to make sure that your input fields do not contain any NaN nor funny values like -999.9<br>
<br>
The input fields should have useful values where there is ocean, and probably zero over land (although that will be masked by the model anyway. If you get a NaN in the first iteration then this is a clear indication that there are still NaNs left somewhere in the input fields.<br>
<br>
Martin<br>
<br>
> On 23. Feb 2020, at 13:34, kunal madkaiker <<a href="mailto:kunal.madkaiker02@gmail.com" target="_blank">kunal.madkaiker02@gmail.com</a>> wrote:<br>
> <br>
> Dear MITgcm Users,<br>
> <br>
> I am trying to run MITgcm along West Coast of India, using a Curvilinear grid.<br>
> I have given climatological profile of Temperature, Salinity and U,V surface wind stress for Dec.<br>
> <br>
> Model is unable to compute RMS-difference for hydrostatic pressure. I tried to run model with <br>
> tRef and sRef values and it gave me result. So, I know problem may lie in my temperature and <br>
> salinity files, but unable to figure out where exactly, since I am able to list and plot these <br>
> netcdf files correctly.<br>
> <br>
> I have also tried converting NaN values into -999.9 in the input files but isnt helping.<br>
> <br>
> STDERR.0059 states,<br>
> (PID.TID 0059.0001) ** WARNING ** MNC_READPARMS: incomplete MNC pickup files implementation<br>
> (PID.TID 0059.0001) ** WARNING ** MNC_READPARMS: => pickup_write_mnc=T not recommanded<br>
> (PID.TID 0059.0001) ** WARNING ** MNC_READPARMS: => pickup_read_mnc=T not working for some set-up<br>
> (PID.TID 0059.0001) ** Warning ** KPP_READPARMS: ignores "data.kpp" file since<br>
> (PID.TID 0059.0001) ** Warning ** KPP_READPARMS: useKPP= F (set from "data.pkg")<br>
> (PID.TID 0059.0001) ** WARNING ** INI_MODEL_IO: globalFiles=TRUE is not safe in Multi-processors (MPI) run<br>
> (PID.TID 0059.0001) ** WARNING ** INI_MODEL_IO: use instead "useSingleCpuIO=.TRUE."<br>
> (PID.TID 0059.0001) ** WARNING ** INI_MODEL_IO: use tiled-files to write sections (for OBCS)<br>
> (PID.TID 0059.0001) *** ERROR *** Initial hydrostatic pressure did not converge after 15 steps<br>
> <br>
> I edited the ini_pressure.F source code and changed steps from 15 to 30 but error persists.<br>
> <br>
> STDERR.0059 states<br>
> (PID.TID 0059.0001) // =======================================================<br>
> (PID.TID 0059.0001) // Check Model config. (CONFIG_CHECK):<br>
> (PID.TID 0059.0001) // CONFIG_CHECK : Normal End<br>
> (PID.TID 0059.0001) // =======================================================<br>
> (PID.TID 0059.0001) <br>
> (PID.TID 0059.0001)  MDS_READ_FIELD: opening global file: new_trial_temp_nov.bin<br>
> (PID.TID 0059.0001)  MDS_READ_FIELD: opening global file: new_trial_salt_nov.bin<br>
> (PID.TID 0059.0001) Start initial hydrostatic pressure computation<br>
> (PID.TID 0059.0001) Iteration   1, RMS-difference =                 NaN<br>
> <br>
> Kindly assist with the error. Thanks in advance!<br>
> Please let me know if any additional information is required.<br>
> <br>
> Regards<br>
> Kunal Madkaiker<br>
> Research Scholar,<br>
> Centre for Atmospheric Science, IIT Delhi<br>
> _______________________________________________<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>
<br>
_______________________________________________<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>