[MITgcm-support] 2D with rotation

[Martin Losch]

Hi Nicolas,

in principle, you should get away with (s)Ny=1, and there are a few  
examples that use this number (e.g. tutorial_plume_on_slope).

We have found that for some reason Ny=1 can lead to numerical  
instabilities after some time on our Vector-computer, that only went  
away, when we increased sNy = 8 or 16. I don't have an explanation,  
except that the overlabs, might be "overlapping" from both sides, if  
the sNy<=2*Oly, but that is not really plausible. Another explanation  
is some screw up due to optimization/vectorization in conjunction with  
the above explanation. The code should work fine for just sNy=1.

Martin

On Dec 7, 2009, at 3:13 PM, Nicolas Grisouard wrote:

> Hi everyone,
>
> I am starting to introduce rotation in my simulations which are  
> still 2D in the zonal plane. What is the minimal value for sNy (# of  
> points in the meridional direction) one has to set for the code to  
> work properly?
>
> I have done a few trials with only 1 point and it seems to work  
> reasonalbly well... once the noice has been removed by a biharmonic  
> operator.
>
> Thanks,
> Nicolas
> -- 
> Nicolas Grisouard
> PhD Student - ERES research group
> Laboratoire des Ecoulements Geophysiques et Industriels
> BP53
> GRENOBLE CEDEX 9 FRANCE
> tel : +33 (0) 476 825 037 - fax : +33 (0) 476 827 022
> http://nicolas.grisouard.free.fr
>
> _______________________________________________
> MITgcm-support mailing list
> MITgcm-support at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-support