Hi There,
I would like to do the correlation of two 2d datasets (MODIS & WRF-output) which are on different grids and trying something like:
;;modis data
lat1 = f1->Latitude
lon1 = f1->Longitude
var1 = f1->var1
var1@lat2d = lat1
var1@lon2d= lon1
;;wrf output
lat = a1->XLAT(0,:,:)
lon = a1->XLONG(0,:,:)
var2 = a1->var2
lat1d = ndtooned(lat) ; (N*M)
lon1d = ndtooned(lon)
ndim = dimsizes(lat)
N = ndim(0)
M = ndim(1)
;;calculate correlation
FO = linint2_points (lon1(0,:), lat1(:,0), var1, False, lon1d, lat1d, 0)
fo = onedtond( ndtooned(FO), (/N,M/) )
rr = pattern_cor(var1,var2,1.0,1)
and get:
-9999
the _FillValue is in both cases -9999, things get wrong before correlation take place FO = -9999
is that the right way to do this and what is wrong ?
Thanks!
Izo
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Mar 26 08:37:23 2012
This archive was generated by hypermail 2.1.8 : Mon Apr 09 2012 - 13:43:03 MDT