[MITgcm-devel] JMD95Z

[Martin Losch]

Patrick,

I am probably not any longer up to date with the code, but as far as I 
know, there are (at least there were) two JMD95 options, 'JMD95Z' for 
the Jacket&McDougall(1995,JAOT)-EOS with pressure=g*rho0*depth, and 
'JMD95P' for the same EOS with exact pressure (lagged by one time 
step). Both work, and if you have a look again at ini_eos.F: JMD95Z may 
not be used with pressure coordinates, because in that case pressure is 
just 'rC', the vertical coordinate. All if-statements, except for the 
one that excludes joint use of JMD95Z and p-coordinate, say 
(equationOfState(1:5) .eq. 'JMD95'), which allows for both JMD95Z and 
P, and that's correct.

Martin

On Tuesday, October 7, 2003, at 11:47  PM, Patrick Heimbach wrote:

>
>
>     I looked into
>     verification/global_ocean.cubed32x32x30/input/data
>     and saw eosType = 'JMD95Z'
>     Is this correct?
>     I see this appearing only in
>     ini_eos.F as a stopping criteria,
>     nowhere else in the code...
>     -p.
>
>     --
>
>     Patrick Heimbach ........................... M.I.T
>     FON: +1/617/253-5259 .......... EAPS, Room 54-1518
>     FAX: +1/617/253-4464 ..... 77 Massachusetts Avenue
>     mailto:heimbach at mit.edu ....... Cambridge MA 02139
>     http://www.mit.edu/~heimbach/ .............. U.S.A
>
>
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://dev.mitgcm.org/mailman/listinfo/mitgcm-devel