[MITgcm-devel] viscC4Leith on the cube

Dimitris Menemenlis menemenlis at sbcglobal.net
Wed Dec 14 00:05:33 EST 2005


Jean-Michel merci beaucoup!  This:

>  useAreaViscLength=.TRUE.,

did the trick and the cubed-sphere config does not blow up anymore.
I have now started a new forward sensitivity integration that includes

>  SadournyCoriolis=.TRUE.,
>  useAbsVorticity=.TRUE.,

as you suggested in an earlier message.

The "sideDragFactor" could be tested later, if you think it might be 
important.  I will also try to contact Baylor about this see if he has 
suggestions.  Right now I have "no_slip_sides=.FALSE.".

I also would like to try the highOrderVorticity scheme.  Have you 
implemented your suggested fix:

> A fix could be to consider the average of the 2,
> Z(2,33)_f2 + Z(1,32)_f2 / 2 for both ?

Dimitris




More information about the MITgcm-devel mailing list