[MITgcm-devel] viscC4Leith on the cube

[Jean-Michel Campin]

Hi Dimitri,

I am looking at this right now.

But seeing what has been changed in dome/input/data,
might be worth to try:
> sideDragFactor=0.,
> useAreaViscLength=.TRUE.,
> viscC4Leith=1.458198138065,     (old value was viscC4Leith=0.1)

and to use the magic forlula (from Baylor) to get the new value 
of viscC4Leith from the old value (viscC4Leith=1.0):
> viscC?leith?_new=pi*(viscC?leith?_old)^(1/3)

Cheers,

Jean-Michel

On Tue, Dec 13, 2005 at 04:36:47PM -0800, Dimitris Menemenlis wrote:
> Jean-Michel and Baylor, here is part II of message regarding problem 
> with viscC4Leith on the cube.
> 
> You can reproduce the problem using the experiment in 
> verification/global_ocean.cs32x15/input
> by adding "viscC4Leith=1.0," to the data file.
> 
> For checkpoint57s_post, this works fine (after setting OLx=OLy=3).
> 
> For checkpoint57t_post, the integration bombs during the first time 
> step, even for deltat=1.
> 
> Dimitris
> 
> -- 
> Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
> Jet Propulsion Lab, California Institute of Technology
> MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
> tel: 818-354-1656;  fax: 818-393-6720
> _______________________________________________
> MITgcm-devel mailing list
> MITgcm-devel at mitgcm.org
> http://mitgcm.org/mailman/listinfo/mitgcm-devel