I've discovered a couple more interesting facts about
this process this morning. It turns out that it is
mostly Triangulate that is giving me problems. I've
found I do NOT have to exclude any values to produce
the proper triangulation, and that the "co-linear"
problem occurs on my Windows box, but not my Linux
box. On Windows, setting the TOLERANCE keyword to 1
appears to solve the problem.
Also interesting is that there is a small gap (mostly
camouflaged in my web article) where the longitude vector
wraps around from 257.5 to 0 degrees. This is especially
apparent in the filled contour method, and less apparent
in the NSIDC regrid method. It appears to disappear
completely in the GridData method, perhaps justifying
my confidence in its power, if we can just learn to
harness it. :-)
I'll update my web page article sometime soon. But
managing to do this with GridData opens up a path
I have been searching for for at least the last
David Fanning, Ph.D.
Fanning Software Consulting, Inc.
Coyote's Guide to IDL Programming: http://www.yqcomputer.com/
Sepore ma de ni thue. ("Perhaps thos speakest truth.")