[MITgcm-support] ​Initial conditions ​At higher initial values for temperature and salinity, model simulation yields nan values.​

[Martin Losch]

The implemented equations of state are empirical functions for the tyical ocean ranges of temperature and salinity, so extreme salinity values may lead to weird densities. But what’s make the model blow up are large density gradients (so large salinity gradients). Are you initial fields smooth?

You can try the availalbe, implemented EOS: eosType = ‘LINEAR’, ‘JMD95Z’, ‘MDJWF’, ’TEOS10’ and see if one of them works, but except for the ‘LINEAR’ one, I expect salinity of 80 to lead to funny density.

Then I would try to find EOSs for high salinity ranges in the literature. I am sure that this has been done in the context of Dead Sea simulations where salinities can go up to 400. Once you have a suitable EOS, you’ll probably have to implement it (see e.g. find_rho.F and related routines).

Martin

> On 23. Feb 2023, at 13:31, Killampalli Rajesh <rajesh70530 at gmail.com> wrote:
> 
> Hi,
> 
> I'm trying to simulate the T,S,U,V in the hypersaline lake, where i'm using field observation of Temperature and Salinity as initial conditions (hydrogThetaFile=temp.bin and hydroSaltFile=salt.bin ). In this study area temp and salinity observation are high values are recorded for supposed salinity ranges from 30-80ppm. The model is blowing at the second time step using this data (error is nan values ). 
> 
> I Think as similar for forcing fields  useExfCheckRange=True/False (check range of input fields and stop if out of range) has to be set for hydrogThetafile and hydroSaltFile. But, for PARM05 to check the range of initial conditions is not available. for this case how to set for higher values. 
> 
> The model runs without any error found, while the field observation data is modified and maintains the range for temp (25-29) and salinity(33-38).
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-support

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20230223/e7d38454/attachment.html>