[MITgcm-devel] exf: u/v-wind interpolation and cubed sphere

[Dimitris Menemenlis]

Martin et al., I am back on e-mail contact.  I have verified that, for both arctic40km and cs510, the North Pole is a corner (XG,YG) point.  In that case the version of exf_set_uv.F <= Revision 1.11 should work fine, no?

The reason for introducing the offset is that the pre-1.11 algorithm cannot accomodate situations where the North Pole is not a corner point and because there was an anomaly in the North-Pole grid definition of the arctic40km grid, now fixed.

For time being, until Ed's interpolation algorithm is in place, shall we go back to exf_set_uv.F Revision 1.11 with a warning if  North Pole is not a corner point (not sure how to check for this in an elegant and foolproof way though)?

A third option is to specify AngleCS and AngleSN a priori and to avoid online agle computations alltogether.  These arrays did not exist at time I wrote exf_set_uv and are not currently specified for cs510 nor for arctic40km.

Dimitris

-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://mitgcm.org/pipermail/mitgcm-devel/attachments/20060623/e807dbe7/attachment.htm>