Problem of Incompatable grids
Hui Du
huidu at UCDAVIS.EDU
Fri Sep 30 14:21:48 EDT 2005
Dear GrADs users:
I want to calculate the correlation between two variables which come
from two data sets. One is model output, the other is reanalysis data.
The domain of the model is a little bigger than the reanalysis, But part
of them are overlapped exactly. I have written the grid set at the last.
We can see that they are exactly the same except the blue part. My
script is
function wp ( args)
tb=subwrd(args,1)
te=subwrd(args,2)
'define avet2a=ave(t2.2,t='tb',t='te')'
'define avet20=ave(t2,t='tb',t='te')'
'define
cort2a=scorr(avet2a,avet20,lon=230.505,lon=252.651,lat=28.4976,lat=46.6943)'
'd cort2a'
return
When I ran the above script, I got:
ga-> run b.gs 613 736
Error from SCORR: Incompatable grids
Operation Error: Error from scorr function
Error ocurred at column 1
DEFINE error: Invalid expression.
I wonder where my problem was, when I calculate two model sets which
have the exactly same grid set. It was ok
Thanks lot
**********************************************************************************************************
Control
longitude = 223.6738, 224.011, 224.3488, 224.6871, 225.0258, 225.3651,
225.7048, 226.045, 226.3856, 226.7267, 227.0682, 227.4102, 227.7525,
228.0952, 228.4384, 228.7819, 229.1257, 229.4699, 229.8145, 230.1593,
230.5045, 230.85, 231.1957, 231.5418, 231.888, 232.2346, 232.5814,
232.9284, 233.2756, 233.623, 233.9706, 234.3183, 234.6662, 235.0143,
235.3625, 235.7108, 236.0592, 236.4077, 236.7563, 237.105, 237.4537,
237.8024, 238.1512, 238.5, 238.8488, 239.1976, 239.5463, 239.895,
240.2437, 240.5923, 240.9408, 241.2892, 241.6375, 241.9857, 242.3338,
242.6817, 243.0294, 243.377, 243.7244, 244.0716, 244.4186, 244.7654,
245.1119, 245.4583, 245.8043, 246.15, 246.4955, 246.8407, 247.1855,
247.5301, 247.8743, 248.2181, 248.5616, 248.9048, 249.2475, 249.5898,
249.9318, 250.2733, 250.6144, 250.955, 251.2952, 251.6349, 251.9742,
252.313, 252.6512, 252.989, 253.3262, 253.6629 ;
latitude = 28.4976, 28.76577, 29.03428, 29.30313, 29.57229, 29.84179,
30.11159, 30.38172, 30.65216, 30.9229, 31.19394, 31.46528, 31.73691,
32.00883, 32.28103, 32.55351, 32.82626, 33.09929, 33.37257, 33.64611,
33.9199, 34.19395, 34.46824, 34.74277, 35.01752, 35.29251, 35.56772,
35.84315, 36.11879, 36.39464, 36.67069, 36.94694, 37.22338, 37.5,
37.7768, 38.05378, 38.33092, 38.60824, 38.8857, 39.16333, 39.44109,
39.71899, 39.99703, 40.2752, 40.55349, 40.83189, 41.1104, 41.38903,
41.66775, 41.94656, 42.22545, 42.50443, 42.78348, 43.0626, 43.34177,
43.62101, 43.90029, 44.17961, 44.45897, 44.73835, 45.01777, 45.29719,
45.57663, 45.85606, 46.1355, 46.41492, 46.69433, 46.9737 ;
Observation
longitude = 230.5045, 230.85, 231.1957, 231.5418, 231.888, 232.2346,
232.5814, 232.9284, 233.2756, 233.623, 233.9706, 234.3183, 234.6662,
235.0143, 235.3625, 235.7108, 236.0592, 236.4077, 236.7563, 237.105,
237.4537, 237.8024, 238.1512, 238.5, 238.8488, 239.1976, 239.5463,
239.895, 240.2437, 240.5923, 240.9408, 241.2892, 241.6375, 241.9857,
242.3338, 242.6817, 243.0294, 243.377, 243.7244, 244.0716, 244.4186,
244.7654, 245.1119, 245.4583, 245.8043, 246.15, 246.4955, 246.8407,
247.1855, 247.5301, 247.8743, 248.2181, 248.5616, 248.9048, 249.2475,
249.5898, 249.9318, 250.2733, 250.6144, 250.955, 251.2952, 251.6349,
251.9742, 252.313, 252.6512 ;
latitude = 28.4976, 28.76577, 29.03428, 29.30313, 29.57229, 29.84179,
30.11159, 30.38172, 30.65216, 30.9229, 31.19394, 31.46528, 31.73691,
32.00883, 32.28103, 32.55351, 32.82626, 33.09929, 33.37257, 33.64611,
33.9199, 34.19395, 34.46824, 34.74277, 35.01752, 35.29251, 35.56772,
35.84315, 36.11879, 36.39464, 36.67069, 36.94694, 37.22338, 37.5,
37.7768, 38.05378, 38.33092, 38.60824, 38.8857, 39.16333, 39.44109,
39.71899, 39.99703, 40.2752, 40.55349, 40.83189, 41.1104, 41.38903,
41.66775, 41.94656, 42.22545, 42.50443, 42.78348, 43.0626, 43.34177,
43.62101, 43.90029, 44.17961, 44.45897, 44.73835, 45.01777, 45.29719,
45.57663, 45.85606, 46.1355, 46.41492, 46.69433 ;
More information about the gradsusr
mailing list