[MITgcm-devel] negative H in calc_r_star.F

[Dimitris Menemenlis]

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?

Dimitris Menemenlis