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
This archive was generated by hypermail 2.1.8 : Fri Nov 01 2013 - 08:58:14 MDT