NCL读取nc文件简单问题求教

妞胖胖

大侠们好：
我才开始学NCL，有一个简单问题想请教大家，大家帮我看看错在哪里吧，谢谢啊。
经过WRF之后，得到了很多output，比如wrfout_d03*这类的，但是他们没有自带后缀。
我就用NCL最简单的命令去读取这些东西，script如下：

begin

   f = addfile("/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_01:00:00.nc","r")
  T = f->T
   printVarSummary(T)

end
这样是可以的，会显示如下:Variable: T
Type: float
Total Size: 1154736 bytes
            288684 values
Number of Dimensions: 4
Dimensions and sizes:   [Time | 1] x [bottom_top | 27] x [south_north | 99] x [west_east | 108]
Coordinates:
Number Of Attributes: 6
  FieldType :   104
  MemoryOrder : XYZ
  description : perturbation potential temperature (theta-t0)
  units :       K
  stagger :
  coordinates : XLONG XLAT

然后我又进行修改，改成
begin

  diri = "/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/"
  fils = systemfunc ("ls " +diri+ "wrfout_d03_2008-11-01_1*.nc")
  f = addfiles(fils,"r")
   T = f->T
   printVarSummary(T)

end
结果就不行了。
错误如下：ls: /nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_1*.nc: No such file or directory
warning:Argument 0 of the current function or procedure was coerced to the appropriate type and thus will not change if the function or procedure modifies its value
fatal:(-2147483647) has no file extension, can't determine type of file to open
fatal:["NclVar.c":1376]:Assignment type mismatch, right hand side can't be coerced to type of left hand side
fatal:Execute: Error occurred at or near line 12 in file file.ncl

或者我把  fils = systemfunc ("ls " +diri+ "wrfout_d03_2008-11-01_1*.nc") 改成 fils = systemfunc ("ls " +diri+ "wrfout_d03_2008-11-01_1*") ,也不行，错误如下：
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_19:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_18:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_17:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_16:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_15:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_14:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_13:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_12:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_11:00:00) has no file extension, can't determine type of file to open
fatal:(/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout_d03_2008-11-01_10:00:00) has no file extension, can't determine type of file to open
fatal:(f) does not reference a file
fatal:Execute: Error occurred at or near line 13 in file file.ncl

我要怎么办呢？那些文件加上后缀说没有该文件，不加后缀又说没有扩展名，我究竟如何解决这个问题呢？
谢谢你们啊。

xo2001

妞胖胖 发表于 2011-12-1 09:28

                               
登录/注册后可看大图

谢谢兰溪。我按照你的第一个方法做了，加上了.nc之后，会有另外的一个提示：(f)does not prefer to a fil ...

T = n[:]->T
就这个冒号~加了才行

言深深

好吧，看起来有点儿复制，深深也不懂

兰溪之水

①把你的wrfout文件后面都加上.nc，然后
fils = systemfunc ("ls /nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout*.nc")
   f    = addfiles (fils, "r")  
或者
②循环读取处理
fs = systemfunc("ls /nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout*00")
  nfs= dimsizes(fs)  
  do ifs=0,nfs-1
    f= addfile

兰溪之水

①把你的wrfout文件后面都加上.nc，然后
fils = systemfunc ("ls /nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout*.nc")
   f    = addfiles (fils, "r")  
或者
②循环读取处理
fs = systemfunc("ls /nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/wrfout*00")
  nfs= dimsizes(fs)  
  do ifs=0,nfs-1
    f= addfile(fs(ifs)+".nc","r")
  end do

妞胖胖

兰溪之水 发表于 2011-12-1 00:26

                               
登录/注册后可看大图

①把你的wrfout文件后面都加上.nc，然后
fils = systemfunc ("ls /nuwa_cluster/home/rcecstu03/RunLiu/W ...

谢谢兰溪。我按照你的第一个方法做了，加上了.nc之后，会有另外的一个提示：(f)does not prefer to a file，后来我又修改了一下，就可以了。
diri = "/nuwa_cluster/home/rcecstu03/RunLiu/WRF/run/PRD/"
   fils = systemfunc("ls "+diri+"wrfout_d03_2008-11-01_0*") + ".nc"
   print(fils)
   n = addfiles(fils,"r")
   T = n[:]->T
   printVarSummary(T)
这样就可以了。不用把wrfout都加上nc了。
谢谢你啊，大半夜的给我解答疑惑。

兰溪之水

妞胖胖 发表于 2011-12-1 09:28

                               
登录/注册后可看大图

谢谢兰溪。我按照你的第一个方法做了，加上了.nc之后，会有另外的一个提示：(f)does not prefer to a fil ...

哦，解决就好~

妞胖胖

兰溪之水 发表于 2011-12-1 09:35

                               
登录/注册后可看大图

哦，解决就好~

{:2_25:}哈哈，谢谢。

da达

蓝溪之水。。。大神啊！

llh

{:5_213:}{:5_213:}{:5_213:}