[MITgcm-devel] negative H in calc_r_star.F

[Martin Losch]

I think you can still use this code if it works. I just could never make it work for my simple experiments. And you Eta looks OK to me (Eta is not the issue anyway). I don't know anything about the negative H.

Martin

On May 25, 2012, at 9:12 PM, Menemenlis, Dimitris (3248) wrote:

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