[MITgcm-devel] negative H in calc_r_star.F

[Martin Losch]

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?
> 
> Dimitris Menemenlis
> 
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel