[MITgcm-devel] diffusion in pkg/seaice

[Gael Forget]

Dear Martin and co,

with regard to seaice_diffusion.F, after discussing the issue with Jean Michel, I am preparing 
to add a diffKh argument, whose units will be m/s2 (as a standard Laplacian diffusion) 
and that will be used instead of DIFF1*MIN( _dxF(I,J,bi,bj), _dyF(I,J,bi,bj)).

Each tracer will be associated with a specific diffusivity (SEAICEdiffKhArea etc.) that will 
be specified alongside the advection scheme in data.seaice, and whose default will be 0.

In the single-dimensional advection scheme case, when advect.F and diffus.F are
used, I am not going to change anything for now. There are several options, and 
I am looking for opinions. We could for example add a call to seaice_diffusion.F
in seaice_advdiff.F with an "if diffKhXXX.NE.0 and diff1.EQ.0" condition. Would this work?

Cheers,
Gael