[MITgcm-devel] pkg/seaice OCEAN_drag in B-grid code

[Martin Losch]

Hi Jean-Michel,

I have to admit that I might have thought about that when writing the C-grid version seaice_dynsolver.F, but I cannot remember what my conclusion was (if there was one). Now this mixing of stress over (marginal) ice does not make sense to me at all. We should remove it (although nobody is *supposed* to use this part of the code anyway (o: ).

Martin


On Jan 16, 2012, at 5:32 PM, Jean-Michel Campin wrote:

> Hi Martin,
> 
> I looked to where OCEAN_drag is used in pkg/seaice: 
> In the C-grid code, assuming winds are provided, OCEAN_drag is
> used only when SEAICE_EXTERNAL_FLUXES is undef (which make sense
> to me, otherwise, EXF compute the wind-stress over open water).
> But in the B-grid code, in dynsolver.F after line 135
>> C now ice surface stress
> it looks like OCEAN_drag has some effect on the wind-stress 
> over seaice (since it's area weighted). Is it intentional ? 
> 
> Cheers,
> Jean-Michel
> 
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel