[MITgcm-devel] SHELFICEBoundaryLayer ustar calculation with hFacC>1

[Martin Losch]

Hi Dan,

clearly, this code was written with hFac<=1 in mind. I think your limiting drKp1 from below makes absolute sense, because this scheme was only meant to work for hFac<1, and as you can see from the code, the largest effect is larger the smaller hFacC is. For a fully wet cell, it should do anything. There are two more occurence of this in shelfice_forcing (there the variable is drLoc).

M.
> On 18 Apr 2015, at 21:56, Daniel Goldberg <dngoldberg at gmail.com> wrote:
> 
> Hi Martin
> 
> At the moment we are using ShelficeBoundaryLayer for the calculation of melt rates, but are allowing hFacC to be larger than 1.
> 
> But we find this allows the factor drKp1 to become negative (line 280 of shelfice_thermodynamics.F) and velocities from the cell below even when the top cell is thicker than drF(kTopC). 
> 
> Is this correct, or was this written with only hFacC(kTopC)<1 in mind? At the moment we are simply inserting the line
> 
> drKp1 = MAX(drKp1, 0.0)
> 
> after line 282 (although long term we hope to implement something like what J-M has suggested, using wet-point averaging to find ustar at the center of the cell).
> 
> Many thanks
> Dan
> 
> 
> 
> -- 
> 
> Daniel Goldberg, PhD
> Lecturer in Glaciology
> School of Geosciences, University of Edinburgh
> Geography Building, Drummond Street, Edinburgh EH8 9XP
> 
> 
> em: Dan.Goldberg at ed.ac.uk
> web: http://ocean.mit.edu/~dgoldberg
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel