[MITgcm-devel] exf wind-stress masking

[Martin Losch]

Hi Jean-Michel,

I am not sure, but I think that for B-grid seaice, the wind (stress) needs to be available a Z-points and this interpolation is done in pkg/seaice/dynsolver.F, so that when masking in exf removes wind values over land, you underestimate the wind over boundary values of seaice.
Having said that, the boundary points of B-grid seaice have (uice,vice)=0 as a boundary conditions, so that this masking does not make any difference, so on second thought you are probably right, and this differentiation can be removed safely.

M.

On Jan 3, 2012, at 5:49 PM, Jean-Michel Campin wrote:

> Hi, 
> 
> I am going to check-in some unfinished changes in exf:
> I realized that the masking of the wind-stress is wrong (when not on 
> C-grid, which happen when interpolated or read from A-grid file).
> I fixed this problem and will check-in this (involves also EXF_PARAMS.h
> & exf_init_fixed.F) and update output of global_with_exf exp.
> 
> I also started to add some params in EXF_PARAMS.h + exf_init_fixed.F
> for 2-components vector interpolation (I did this before the stress
> mask fix), will be part of the changes but they are not yet used.
> 
> Regarding masking, I don't see a clear reason why it's only
> applied when not using SEAICE. If exf works fine without
> masking (when using SEAICE), then why do we care about applying
> masks ?
> 
> Also, noticed that exf_init_fixed.F does not compile with default 
> pkg/icefront settings (might have been better to put all this
> icefront - exf code in pkg/icefront and still call some utility 
> exf & cal S/R).
> 
> Cheers,
> Jean-Michel
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel