[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