[MITgcm-support] STOP ABNORMAL END: S/R INI_PRESSURE

[Martin Losch]

Hi Kunal,

this iteration solves the hydrostatic pressure equation (dp = - g rho(z) dz) and usually converges very quickly.

I think you need to make sure that your input fields do not contain any NaN nor funny values like -999.9

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.

Martin

> On 23. Feb 2020, at 13:34, kunal madkaiker <kunal.madkaiker02 at gmail.com> wrote:
> 
> Dear MITgcm Users,
> 
> I am trying to run MITgcm along West Coast of India, using a Curvilinear grid.
> I have given climatological profile of Temperature, Salinity and U,V surface wind stress for Dec.
> 
> Model is unable to compute RMS-difference for hydrostatic pressure. I tried to run model with 
> tRef and sRef values and it gave me result. So, I know problem may lie in my temperature and 
> salinity files, but unable to figure out where exactly, since I am able to list and plot these 
> netcdf files correctly.
> 
> I have also tried converting NaN values into -999.9 in the input files but isnt helping.
> 
> STDERR.0059 states,
> (PID.TID 0059.0001) ** WARNING ** MNC_READPARMS: incomplete MNC pickup files implementation
> (PID.TID 0059.0001) ** WARNING ** MNC_READPARMS: => pickup_write_mnc=T not recommanded
> (PID.TID 0059.0001) ** WARNING ** MNC_READPARMS: => pickup_read_mnc=T not working for some set-up
> (PID.TID 0059.0001) ** Warning ** KPP_READPARMS: ignores "data.kpp" file since
> (PID.TID 0059.0001) ** Warning ** KPP_READPARMS: useKPP= F (set from "data.pkg")
> (PID.TID 0059.0001) ** WARNING ** INI_MODEL_IO: globalFiles=TRUE is not safe in Multi-processors (MPI) run
> (PID.TID 0059.0001) ** WARNING ** INI_MODEL_IO: use instead "useSingleCpuIO=.TRUE."
> (PID.TID 0059.0001) ** WARNING ** INI_MODEL_IO: use tiled-files to write sections (for OBCS)
> (PID.TID 0059.0001) *** ERROR *** Initial hydrostatic pressure did not converge after 15 steps
> 
> I edited the ini_pressure.F source code and changed steps from 15 to 30 but error persists.
> 
> STDERR.0059 states
> (PID.TID 0059.0001) // =======================================================
> (PID.TID 0059.0001) // Check Model config. (CONFIG_CHECK):
> (PID.TID 0059.0001) // CONFIG_CHECK : Normal End
> (PID.TID 0059.0001) // =======================================================
> (PID.TID 0059.0001) 
> (PID.TID 0059.0001)  MDS_READ_FIELD: opening global file: new_trial_temp_nov.bin
> (PID.TID 0059.0001)  MDS_READ_FIELD: opening global file: new_trial_salt_nov.bin
> (PID.TID 0059.0001) Start initial hydrostatic pressure computation
> (PID.TID 0059.0001) Iteration   1, RMS-difference =                 NaN
> 
> Kindly assist with the error. Thanks in advance!
> Please let me know if any additional information is required.
> 
> Regards
> Kunal Madkaiker
> Research Scholar,
> Centre for Atmospheric Science, IIT Delhi
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support