ncl怎么进行批处理？

tianmeng

风之牧语 发表于 2014-10-28 17:09
rcm2rgrid插值后不能用plot = wrf_map_overlays(a,wks,(/contour /),pltres,mpres)这个语句画图了！用gsn ...

难道是我的版本太低，不能批处理吗？我用的是ncl6.0.0

风之牧语

tianmeng 发表于 2014-10-28 17:17
难道是我的版本太低，不能批处理吗？我用的是ncl6.0.0

   files=systemfunc("ls -1 wrfout_d01_2013-07*")
   f=addfiles(files,"r")
   ListSetType(f,"cat")这个是我在数据相同目录的写法！

风之牧语

tianmeng 发表于 2014-10-28 17:16
错误是：，fatal:_NclBuildConcatArray: Arrays of files are not yet supported
fatal:Execute: Error  ...

我看你程序主体和我以前发的帖子一模一样，应该没有什么错误吧！不过olon = fspan(114,123,144)
      olat = fspan(34,39,80)你这个格点插值的也太细了吧！还是要他说了才权威@longlivehj

tianmeng

风之牧语 发表于 2014-10-28 17:18
files=systemfunc("ls -1 wrfout_d01_2013-07*")
   f=addfiles(files,"r")
   ListSetType(f,"cat" ...

ls为什么要加-1，另外，我改成你的，还是不行，错误还是fatal:_NclBuildConcatArray: Arrays of files are not yet supported
fatal:Execute: Error occurred at or near line 45 in file wrf_Height.ncl

tianmeng

是不是ncl_ncarg-6.0.0.Linux_Debian_x86_64_nodap_gcc445.tar.gz,这个版本不提供批量处理呀！

tianmeng

@longlivehj

风之牧语

tianmeng 发表于 2014-10-28 17:29
@longlivehj

我都不知道你错的那一行是哪一行？

tianmeng

风之牧语 发表于 2014-10-28 17:30
我都不知道你错的那一行是哪一行？

times = wrf_user_getvar(f[:],"times",-1)  ; get all times in the file

tianmeng

风之牧语 发表于 2014-10-28 17:30
我都不知道你错的那一行是哪一行？

我觉得没错呀！我把文本放在我师兄的服务器转还是一样的错误

风之牧语

tianmeng 发表于 2014-10-28 17:33
times = wrf_user_getvar(f[:],"times",-1)  ; get all times in the file

你确定你是用的WRFOUT的数据？