correlation between data from txt file and gridded sst data

Hi all, I want to plot a spatial distribution of the correlation between the sst data and my station data. I try by the following code;   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" begin fili = "Jeddah_Rain_1982_2013.txt" f1  = addfile("sst.mnmean.anom_1982-2013.cdf","r") sst = f1->ssta lat = f1->Y lon = f1->X ; initialize site and year count ; note, set to work for six years of data only mncnt   = 0 yrcnt   = 0 sitecnt = 0 ntime = 380 ;---Read in file as array of strings so we can parse each line.   lines  = asciiread(fili,-1,"string")   nlines = dimsizes(lines)-1   ; First line is a header   delim     = ","   rr_obs     = str_get_field(lines(1:),3,delim)   rr_obs@_FillValue = -999 copy_VarCoords(sst,rr_obs) printVarSummary(rr_obs) ;print(rr_obs+"  ") ssta = sst(Y|:,X|:,time|:)   corel  = esccr(rr_obs,ssta,1) ;  copy_VarCoords(sst(0,:,:), corel) print(corel+"  ") end  and I had got the error; fatal:Argument type mismatch on argument (0) of (esccr) can not coerce fatal:["Execute.c":8128]:Execute: Error occurred at or near line 34 in file statistics_tmp.ncl the ascii file has three columns year, month and the rainfall. So, can anyone help me? Thanks Regards   Mohamed Ibrahim El-Sayed OMAR Egyptian Meteorological Authority Scientific Research Department Climate Studies Centre Kobry El-Quobba, P.O.Box. 11784 Cairo, EGYPT Mobil phone  : 002 01060802805  Home Phone : 002 0482370261

_______________________________________________
ncl-talk mailing list
List instructions, subscriber options, unsubscribe:
http://mailman.ucar.edu/mailman/listinfo/ncl-talk
Received on Thu Oct 31 04:59:20 2013