[ncl-talk] warning:["FileSupport.c":3894]:NCL does not know anything about file suffix <csv>.
grace
313695096 at qq.com
Sun Jun 19 21:53:53 MDT 2016
Hi:
All,I have write a script to read and calculate .csv data according to the "Reading csv"example on the NCL website,but it appears warning and fatal:
warning:["FileSupport.c":3894]:NCL does not know anything about file suffix <csv>.
But NCL will try its best to figure out the file format.
fatal:An internal error has occurred. The file format requested does not appear to be supported, could not open (/public/home/huanglei/cmaqdata/20160615/Caotan.csv)
warning:["FileSupport.c":3894]:NCL does not know anything about file suffix <csv>.
But NCL will try its best to figure out the file format.
........
.........
........
This is my script:
;
begin
filename1=systemfunc("ls /public/home/huanglei/cmaqdata/2016*/Caotan.csv")
fin1=addfiles(filename1,"r")
time_al= new((/29,24/),float,"No_FillValue")
oth_al = new((/29,24/),float,"No_FillValue")
xa_al = new((/29,24/),float,"No_FillValue")
tc_al = new((/29,24/),float,"No_FillValue")
bj_al = new((/29,24/),float,"No_FillValue")
xy_al = new((/29,24/),float,"No_FillValue")
wn_al = new((/29,24/),float,"No_FillValue")
ya_al = new((/29,24/),float,"No_FillValue")
hz_al = new((/29,24/),float,"No_FillValue")
yl_al = new((/29,24/),float,"No_FillValue")
ak_al = new((/29,24/),float,"No_FillValue")
sl_al = new((/29,24/),float,"No_FillValue")
do i=0,28
;---Read in file as array of strings so we can parse each line.
lines = asciiread(fin1[i],-1,"string")
nlines = dimsizes(lines)-1 ; First line is a header
;
; Start reading in the fields we care about:
; The first line is a header, so skip this.
;
delim = ","
time = str_get_field(lines(1:25),1,delim)
oth = str_get_field(lines(1:25),2,delim)
xa = str_get_field(lines(1:25),3,delim)
tc = str_get_field(lines(1:25),4,delim)
bj = str_get_field(lines(1:25),5,delim)
xy = str_get_field(lines(1:25),6,delim)
wn = str_get_field(lines(1:25),7,delim)
ya = str_get_field(lines(1:25),8,delim)
hz = str_get_field(lines(1:25),9,delim)
yl = str_get_field(lines(1:25),10,delim)
ak = str_get_field(lines(1:25),11,delim)
sl = str_get_field(lines(1:25),12,delim)
time_al(i,:)= time
oth_al(i,:) = oth
xa_al(i,:) = xa
tc_al(i,:) = tc
bj_al(i,:) = bj
xy_al(i,:) = xy
wn_al(i,:) = wn
ya_al(i,:) = ya
hz_al(i,:) = tc
yl_al(i,:) = bj
ak_al(i,:) = xy
sl_al(i,:) = wn
end do
time_all=ndtooned(time_al)
oth_all=ndtooned(oth_al)
xa_all=ndtooned(xa_al)
tc_all=ndtooned(tc_al)
bj_all=ndtooned(bj_al)
xy_all=ndtooned(xy_al)
wn_all=ndtooned(wn_al)
ya_all=ndtooned(xy_al)
hz_all=ndtooned(wn_al)
yl_all=ndtooned(ya_al)
ak_all=ndtooned(ak_al)
sl_all=ndtooned(sl_al)
print(time_all)
; calculate the avg
oth_avg=avg(oth_al)
xa_avg=avg(xa_all)
tc_avg=avg(tc_all)
bj_avg=avg(bj_all)
xy_avg=avg(xy_all)
wn_avg=avg(wn_all)
ya_avg=avg(xy_all)
hz_avg=avg(wn_all)
ya_avg=avg(ya_all)
ak_avg=avg(ak_all)
sl_avg=avg(sl_all)
ssx = hz_avg + ak_avg + sl_avg
nsx = ya_avg + ya_avg
end
My NCL version is 6.2.0.
How can I slove the problem?
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/ncl-talk/attachments/20160620/b9f9ccb1/attachment.html
More information about the ncl-talk
mailing list