[MITgcm-support] Format of curvilinear grid files

[Christopher L. Wolfe]

Hi all,

I've been trying to set up an experiment on a curvilinear grid with pre-computed grid spacings, areas, and angles. I thought I had written out all the grid descriptors in the proper format (18 descriptors, big endian, 64-bit), but I kept having the model try to read past the end of my grid file and crashing. Some poking around on the MITgcm website brought me to the Held-Suarez-on-a-cube-sphere tutorial (http://mitgcm.org/public/r2_manual/latest/online_documents/node137.html) where it mentions that the grid faces are 32x32 but the grid files contain 18 33x33 fields. This made me think that I totally don't understand what's supposed to be in the grid files.

Can anyone confirm that the grid files are indeed supposed to contain one extra point in each direction? If so, what defines the values of the grid descriptors on that extra point?

In any event, if anyone can point me to a matlab script or other such thing with an example of writing out curvilinear grid descriptors, I'd very much appreciate it.

Thanks!
Christopher

-----------------------------------------------------------
Dr. Christopher L. P. Wolfe
Research Oceanographer
Climate, Atmospheric Science, and Physical Oceanography
Scripps Institution of Oceanography, UCSD
clwolfe at ucsd.edu                           858-534-4560
-----------------------------------------------------------





-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-support/attachments/20130702/2400dcc7/attachment.htm>