[MITgcm-support] How to add a new equation of state in mitgcm

[赖燕红]

Hi, 

I want to add a new equation of state in MITgcm. Particularly, this equation of state is assumed to be a piecewise function, that means the density has a different dependence on pressure and temperature under different temperatures. I found some files in /models/ that might relate to the settings of the equation of state, for example, /models/src/ini_eos.F & Find_rho.F & find_alpha.F, /models/inc/EOS.h & PARMS.h. But I'm still very confused on how to realize a piecewise equation of state that depends on temperature. Does anyone know or ever try how to realize this? 




Also, if I want to add a new parameter in /input/data PARM01, where do I need to define this new parameter? Or do I need to modify other files at the same time?




Thanks in advance!




Yanhong 
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mailman.mitgcm.org/pipermail/mitgcm-support/attachments/20211217/1bb72e23/attachment.html>