Spherical gridding at a pole [message #4787] |
Tue, 25 July 1995 00:00  |
ECSPRS
Messages: 6 Registered: November 1994
|
Junior Member |
|
|
Hi.
I'm having some problems with the TRIANGULATE function when trying to grid some
irregular data over the south pole. I have data points defined for each hour
of the day (hence each 15 deg in longitude) and for latitudes -90 to -50. So
the data looks like the spokes of a wheel radiating out from the pole. Now I
wanted to grid this and make a contour plot on top of a map of Antarctica so I
used the spherical implementation of TRIANGULATE vis;
TRIANGULATE,lon,lat,tr,SPHERE=s,FVALUE=epot,/DEGREES
Where lon, lat and epot are vector arrays of longitude, latitude and data.
This statement causes either IDL to crash or Windows to crash (wow...how
surprising). Any ideas? It just occured to me that the point -90 latitude
will appear 24 times in the array, once for each longitude. Would this help
mess things up? (I'll try taking the extra ones out out but I'd rather send
this first before disaster strikes).
I also noticed that while in the doco the GS, Limits parameters are in brackets
as optional for TRIGRID (spherical gridding mode), the command won't work
unless you put them in there.
I'm using IDL 4.0
Cheers,
Paul
|
|
|