<p>
        Hi, 
</p>
<p>
        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? 
</p>
<p>
        <br>
</p>
<p>
        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?
</p>
<p>
        <br>
</p>
<p>
        Thanks in advance!
</p>
<p>
        <br>
</p>
<p>
        Yanhong 
</p>