Question on the correlation pattern plot ?

Hi, all

   I have a question when I want to get a correlation field between PDO and NH precipitation.

  1. I met an error when dealing with rtest of cross correlation ,
error info : "SLATEC/NCL: DBETAI: P AND/OR Q IS LE ZERO: NERR= 2: LEVEL= 2"
and codes were listed below :

  ccr1 = escorc(pdo({1951:2002}),ann_cru(:,:,{1951:2002})) ; CRU data
  nlat = dimsizes(ccr1(:,1)) ; nlat= 180 for North Hemisphere
  nlon = dimsizes(ccr1(1,:)) ; nlon = 720
  print(nlon)
  print(nlat)
  Nx1 = new((/nlat,nlon/),"integer")
  do ilat = 0,nlat-1
     do ilon = 0,nlon-1
        Nx1(ilat,ilon) = num(.not.ismissing(ann_cru(ilat,ilon,{1951:2002})))
     end do
  end do
  prob1 = rtest(ccr1,Nx1,0)

I believe that it is because Nx1 have many constant missing value in ocean area that cause such trouble.
But I can not find a way to resolve it.

Thanks

Lin

LEO_ARIES
2007-01-22

_______________________________________________
ncl-talk mailing list
ncl-talk_at_ucar.edu
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Mon Jan 22 2007 - 08:52:02 MST