[MITgcm-devel] negative H in calc_r_star.F

[Menemenlis, Dimitris (3248)]

Ooops!

This means that for past three years we have boldly been going where we had no business of going :-(

http://ecco2.jpl.nasa.gov/data1/cube/cube86/input/ETANjan08.jpg
is a figure of ETAN on the Antarctic face of the cube86 global CS510 configuration after 16 years of integration.

Apart from fact that one can end up with negative H at edges of grid cells, which was not the case for cube86,
what other aspects of the existing "r_star + shelfice"  solutions are likely to be problematic?

Shall we add a stop or a warning to shelfice_check.F when r_star is enabled?

Dimitris Menemenlis

On May 25, 2012, at 12:52 AM, Martin Losch wrote:

> I am not aware, that shelfice can be used with r-star! 
> 
> M.
> 
> On May 24, 2012, at 9:30 PM, Dimitris Menemenlis wrote:
> 
>> Jean-Michel and Martin, is it allowable to have negative H in calc_r_star.F  when using pkg/shelfice ?
>> 
>> When we went from hFacMinDr=50 to hFacMinDr=0 in order to better represent the ice shelf cavity,
>> we end up with following situation:
>> 
>> WARNING: r*FacW < hFacInf at       1 pts : bi,bj,Thid,Iter=   1   1   1        10
>> e.g. at i,j=   9  99 ; rStarFac,H,eta =  0.019415 -1.499939E-01  1.689013E-01  1.252815E-01
>> STOP in CALC_R_STAR : too SMALL rStarFacW !
>> 
>> I am pretty sure that H is non-zero in center of cell - Hong is double checking.
>> ... but shouldn't H at the edges of the cells be set to zero when
>> MIN( Ro_surf(i-1,j), Ro_surf(i,j) ) - MAX( R_low(i-1,j), R_low(i,j) )
>> is negative?