[MITgcm-devel] selectP_inEOS_Zc

[Martin Losch]

Hi Jean-Michel,

Dustin and I are trying to debug his configuration with the shelfice package. We are struggling to make the initial conditions “balance” (i.e. not have large initial adjustments). 
I can do it for eosType=‘JMD95z’, but not for JMD95p and MDJWF. The latter two use selectP_inEOS_Zc = 2 and the former uses 0 (very simple pressure = g*rho*zc)
When I look at the results of phiHyd (i.e. totPhiHyd as stored in phiHyd.*.nc), then I see that for selectP_inEOS_Zc = 2, phiHyd = 0 within the ice shelf, but for selectP_inEOS_Zc = 0 is corresponds to the value that is prescribed in SHELFICEloadAnomalyFile. I am afraid, that as a consequence, the phiHyd at the bottom of the ice shelf is not correct for selectP_inEOS_Zc = 2 (at least it looks like that in timestep 0: there’s a lateral gradient that I want to avoid with prescribing SHELFICEloadAnomaly).

Can you remember, why there is this different behavior for the different selectP_inEOS_Zc?
I would expect the totPhiHyd to have the contribution of SHELFICEloadAnomalyFile (in the end it’s added as phi0surf in diags_phiHyd) in all cases.

Martin