point-correlation

Hi Everyone,

I am trying to plot a cross co-relation map of point value with gridded data. The following error is troubling.

fatal:Subscript out of range, error in subscript #2
fatal:An error occurred reading y
fatal:Execute: Error occurred at or near line 32 in file point_corr.ncl

and the script is as follows:

load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/contributed.ncl"
load "$NCARG_ROOT/lib/ncarg/nclscripts/csm/shea_util.ncl"

begin
;**********************************************
; Reading Data and correlation function
;**********************************************
        f1 = addfile("6104_01.nc","r")
        f2 = addfile("tll_6104_01.nc","r")
        x1 = f1->a
        x2 = f2->a
        print(x1)
        print(x2)
        y = x1(lon|:,lat|:,time|:)
        c = escorc(y(:,:,0),x2)
        printVarSummary(c)

;**********************************************
; create plot
;*********************************************
  wks = gsn_open_wks ("ps","cor_prcp_01") ; open ps file

  gsn_define_colormap(wks,"radar")
  res = True
  res_at_cnFillOn = True
  res_at_cnLinesOn = False

  res_at_gsnSpreadColors = True
  res_at_lbAutoLabelStride = True

    res_at_cnLevelselectionMode = "ManualLevels" ; manually set cn levels
  res_at_cnMinLevelValF = -1. ; min level
  res_at_cnMaxLevelValF = 1. ; max level
  res_at_tiMainString = ":F25:Prcp Point Corr"

  res_at_gsnLeftString = " "
  res_at_tiMainFontHeightF = 0.03
  res_at_gsnAddCyclic = False
  plot = gsn_csm_contour_map_ce(wks,c,res)
end

It will be really great if anyone can help me out.
Thanks

-Ibraheem
_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Tue Jun 23 2009 - 21:29:18 MDT