[MITgcm-devel] negative H in calc_r_star.F

[Dimitris Menemenlis]

Jean-Michel, I notice your CVS check-ins this morning.
Shall I take the code out for a spin on the CS510 - or wait a little longer?

Looking again at the figure that I sent last night, I notice that it 
includes an
example of a pathological case that neither (1) nor (2) can deal with.
See bathymetry and ice shelf topography detail circled in red in attached.
I guess this is a case of "careless preparation of topographies" and the
only fix is to remove this "pocket" of water?

Thanks

Dimitris Menemenlis


On 05/31/2012 02:06 PM, Jean-Michel Campin wrote:
> Hi,
>
> I will make the little changes for method (1) (will be needed anyway for method (2))
> in calc_surf_dr.F and within IF (useShelfIce) THEN / ENDIF in ini_masks_etc.F
> so that results will not change (machine trucation) if not using shelfice.
> Will see later how method (2) can be implemented.
>
> There was also the idea of allowing hFac>  1 near the bottom instead
> of rounding off the depth if hFac below is too small.
> Might do both at the same time.
>
> Cheers,
> Jean-Michel
-------------- next part --------------
A non-text attachment was scrubbed...
Name: NegativeDepth.pdf
Type: application/pdf
Size: 11731 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20120601/bc6e45fe/attachment-0001.pdf>