[MITgcm-devel] negative H in calc_r_star.F

[Jean-Michel Campin]

Hi Dimitris,

Yes, I would give it a try. I think we have soltion (1) implemented
(but I did not really tested it) so it will be interesting to see how it goes.

I might still make little changes in calc_surf_dr.F but with rStar you 
are not using it, so should not matter.

And regarding the grid-point in the red circle in the pdf,
I imagine this plot as beeing a section, so the pocket of water might 
still be connected to open ocean in the other direction ? 
Or is this too optimistic ?

Cheers,
Jean-Michel

On Fri, Jun 01, 2012 at 08:27:42AM -0700, Dimitris Menemenlis wrote:
> 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


> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel