[MITgcm-devel] selectP_inEOS_Zc

[Jean-Michel Campin]

Hi Martin,

I will take a look at this soon, and will let you know.

Cheers,
Jean-Michel

On Thu, Jan 11, 2018 at 05:20:34PM +0100, Martin Losch wrote:
> 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
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mailman.mitgcm.org/mailman/listinfo/mitgcm-devel