[Thread Prev][Thread Next][Index]
Correlation w/Two diff Datasets
Howdy ferreteers,
I have been converting the variance script I found from the initial
install for use with my data sets. Unfortunately, I have a problem with
the data sets being unequal in size as well as not on the same gridpts
and different resolutions:
************************************************************
./TRMM_SST_concat.nc (resolution of 0.25km or about a quarter degree)
name title I J
K L
SST SST 1:1440 1:320 ...
1:762
X=0E:0E(360) Y=-40:40
Scattereometer.nc (1 degree resolution)
name title I J K
L
U zonal pseudostress at 10m 1:360 1:171 ...
1:839
X=0.5W(-0.5):0.5W Y=85.5S:85.5N
V meridional pseudostress at 10m 1:360 1:171 ...
1:839
X=0.5W(-0.5):0.5W Y=85.5S:85.5N
precip97_01.nc (1 degree resolution)
name title I J
K L
PRECIP Daily Precip 1:360 1:180 ...
1:1704
X=0E:0E(360) Y=90S:90N
************************************************************
I thought the answer would be to specify my regions with x, y and
tlimits in my script:
***********************************************************
Set mode interpolate
use "/rlp/scocks/qscat_99_01mask.nc"
use "/rlp/scocks/precip97_01.nc"
set window/clear/size=0.75/aspect=0.8
Set Region/x=120e:70w/y=60s:60n
let magnitude=(U[d=1]^2+V[d=1]^2)^0.25
let uwnd=U[d=1]/magnitude; let vwnd=V[d=1]/magnitude
let r=6.38e6 ; let omega=7.292e-5 ; let d2r=atan(1.0)/45
let Div=(uwnd[x=@ddc]+vwnd[y=@ddc]-(vwnd/r)*tan(y[d=1,g=u]*d2r))
LET p = Precip[d=2, x=120e:70w, y=35s:35n]; LET q = Div[d=1, x=120e:70w,
y=35s:35n]
*****************************************************************
But I got the following error from the last step:
l
**ERROR: inconsistent sizes of data regions: X axis
Variable 2 is I=171:181
expression is I=171:180
I specifically looked for this type of problem on the archive and didn't
see it. Does anyone know how to solve this problem?
SCOCKS
TAMU
[Thread Prev][Thread Next][Index]
Dept of Commerce /
NOAA /
OAR /
PMEL /
TMAP
Contact Us | Privacy Policy | Disclaimer | Accessibility Statement