[MITgcm-devel] cube-sphere question

[Dimitris Menemenlis]

I may have found a possible cause for mysterious crashes
near cube-face edges as viscosity is decreased.  I think
there is a problem with grid definition both for the
32x32x6 case as well as for the hi-res case.

Shouldn't DXC(1,:,:) equal DXC(cx,:,:) ?
The grid as defined now has DXC(1,:,:) < DXC(cx,:,:)
(see attached figure fro 32x32x6 case).

Conversely, shouldn't DXG(1,:,:) < DXG(cx,:,:) ?
The grid as defined now has DXG(1,:,:) = DXG(cx,:,:)

Please confirm whether this is an issue or not.

Cheers, D.

-- 
Dimitris Menemenlis <menemenlis at jpl.nasa.gov>
Jet Propulsion Lab, California Insitute of Technology
MS 300-323, 4800 Oak Grove Dr, Pasadena CA 91109-8099
tel: 818-354-1656;  fax: 818-393-6720
-------------- next part --------------
A non-text attachment was scrubbed...
Name: FIG.ps.gz
Type: application/x-gzip
Size: 2870 bytes
Desc: not available
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20040221/12524f87/attachment.gz>