[Mitgcm-support] Re: horizontal viscosity

mitgcm-support at dev.mitgcm.org mitgcm-support at dev.mitgcm.org
Wed Jul 9 15:28:07 EDT 2003


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            



More information about the MITgcm-support mailing list