ncl批量读取txt文件问题

fealbai

ncl程序：
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_code.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/gsn_csm.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/contributed.ncl"
load "/usr/local/lib/ncarg/nclscripts/csm/shea_util.ncl"
begin

diri = "/home/Administrator/160/"
fils = systemfunc ("ls "+ diri + "16*")
print(fils)
f   = addfiles(fils,"r")  
data = asciiread("fils",(/488,20/),"float")
print(data)
end
：-------------------------------------------------------------------------------------
12个降水资料，txt格式，想批量读取
之前测试过 可以批量读取*.nc资料
是不是addfiles的问题
望各位大侠指点，不胜感激！

侠骨柔情

如果每个文件的大小不一样该如何处理？
begin
files=systemfunc("ls ./*.storm.track")
print(dimsizes(files))
nrow=new(13,"integer")
do i=0,dimsizes(files)-1
nrow(i)=numAsciiRow(files(i))
end do

do i=0,dimsizes(files)-1
data=asciiread(files(i),(/nrow(i),4/),"float")
print(data)
end do

然后报错：
fatal:Dimension sizes of left hand side and right hand side of assignment do not match
fatal:["Execute.c":8578]:Execute: Error occurred at or near line 21 in file storm.track.ncl

hzwjy

addfiles不是给TXT文件用的。

fils = systemfunc ("ls "+ diri + "16*")
print(fils)
do i=0, dimsizes(fils)-1
data = asciiread(fils(i),(/488,20/),"float")
print(data)
end do

fealbai

hzwjy 发表于 2012-9-27 20:11

                               
登录/注册后可看大图

addfiles不是给TXT文件用的。

fils = systemfunc ("ls "+ diri + "16*")

问题已解决，非常感谢大侠的帮忙，

realmadridhz

fealbai 发表于 2012-9-27 21:01

                               
登录/注册后可看大图

问题已解决，非常感谢大侠的帮忙，

怎么解决我也想知道，我在处理JTWC的best-track资料要用到这个

maliyasunny

同样的问题已经解决，谢谢啊~~~

xiaomo

hzwjy 发表于 2012-9-27 20:11
addfiles不是给TXT文件用的。

fils = systemfunc ("ls "+ diri + "16*")

赞一个，也为同样问题苦恼

Cash

学习了！{:eb502:}{:eb502:}{:eb502:}{:eb502:}

shawnSS

能有什么方法不用循环读么？ 类似于addfile批量读的方法

白兔糖0217

极好！

miniphoenix

以后要常来这里啊，有人一起讨论真好