Hi all,
I am trying to put irregularly spaced geo data onto a grid. Doesn't
seem to work with the test data below, at least not in spherical
gridding mode. The range of values is spread massively, cmp. original
data with results. This is what I get
IDL> print,min(data,max=md),md
2.01690 3.04770
IDL> print,min(result,max=mr),mr
-35.4423 14.1332
How can that happen? Help very much appreciated
Carsten
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
; Paste this into your idl comand line
; This is some test data (part of an actual ozone data set):
data = [2.4373, 2.7141, 2.3265, 2.8369, 2.7715, 2.3347, 2.3997, $
2.0169, 3.0362, 2.7579, 2.8260, 2.1901, 3.0477, 2.8023]
lat = [40.9423, 43.8911, 43.3905, 46.0351, 47.9382, 49.0647, 48.9252, $
52.8167, 52.0945, 55.3319, 56.3629, 58.4239, 59.8312, 59.0832]
lon = [94.3195, 82.7758, 91.3826, 87.9544, 80.6488, 81.2267, 85.4324, $
83.4327, 84.4105, 84.7905, 85.4792, 86.5243, 87.5318, 97.5873]
limits = [80.00, 40.00,100.00, 60.00]
spacing = [2.0,2.0]
; Do a spherical triangulation:
triangulate,lon,lat,tr,/degree,sphere=sphere,fvalue=data
result=trigrid(data,spacing,limits,sphere=sphere,/degree, $
xgrid=gx,ygrid=gy)
; look at the data / result:
print,min(data,max=md),md
print,min(result,max=mr),mr
; end of test script
;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;;
|
comp.lang.idl-pvwave archive
Members
Search
Help
Login
Home
