[Mitgcm-support] Re: horizontal viscosity

[mitgcm-support at dev.mitgcm.org]

Hi Martin,

I agree with Chris.

Also, i think that 2D arrays make sense (like the Visbeck K), and do not 
require too much space so that, even if someone is not using it, 
it does not matter.
If its 3 D arrays (viscAh viscA4 + Temp, Salin, Tracers, ...) it would
be good to use a CPP_OPTION, so that the size of the code do not 
increase when you know you will not use it.

Cheers,

Jean-Michel

> Hi,
>
> are there any plans to replace the so-far constant horizontal
> viscosities, viscAh and viscA4, by 2D or 3D fields? So that one could
> e.g., scale the coefficients with variable grid size, or implement
> something like a Smagorinski-type viscosity? If I were to
> code this up,
> would it make sense to go for (global) 2D fields
> Ah(i,j,bi,bj) or should
> I immediately introduce 3D fields Ah(i,j,k,bi,bj)? Is this
> anything of
> interest to anyone else or is it just me (in which case I
> would do as I
> please and not check anything in (o: )?
>
> Martin