求助：关于批处理nc数据的问题

MacTaylor

有一批nc数据，一个月一个文件，现将文件名改为：197901.nc，197902.nc，依次类推，对二进制的数据知道写ctl文件来进行批处理运算，那么nc数据应该怎么做呢？试过了写ctl文件，sdfopen ***\%y4%m2.nc，
year=1978
while(year<=2011)
month=01
  while(month<=12)
'sdfopen ***\'year''month'.nc'
...
    month=month+1
  endwhile
year=year+1
endwhile
都没有成功，请教大家这个应该怎么做啊？

小虹

批量读取nc数据，用你的方法成功了，谢谢！！！

MacTaylor

kc121 发表于 2012-6-1 16:18

                               
登录/注册后可看大图

具体方法是什么样的，楼主给我们解释一下吧。不要只顾着自己解决问题了就可以了啥。还是分享一下吧。谢谢！

其实我也是糊里糊涂的解决了。。。ctl文件如下：
dset  ^%y4%m2.nc
undef 1e+15
options template
title MERRA data
dtype NetCDF
ydef  144 linear -90 1.25
xdef  288 linear -180 1.25
zdef   21 levels 1000 975 950 925 900 875 850 825 800 775 750 725 700 650 600 550 500 450 400 350 300
tdef  396 linear 00Z01JAN1979 1mo
vars 3
qv 21 t,z,y,x Specific humidity
u 21 t,z,y,x Eastward wind component
v 21 t,z,y,x Northward wind component
endvars
然后open ***.ctl就行了，之前的问题是打不开ctl文件，怎么改也不行，后来换了台机子就好了。。。所以我说我也不知道怎么回事，希望对你有帮助。

易小凯

写成%y4%m2.nc的形式就行

mofangbao

你把nc文件配个ctl批量描述下

MacTaylor

mofangbao 发表于 2012-5-6 16:16

                               
登录/注册后可看大图

你把nc文件配个ctl批量描述下

我之前就这样做的，能打开ctl文件，但是d之后，都显示all undefined values，我的ctl如下，麻烦帮我看看哪里错了？
dset  ^%y4%m2.nc
undef 1e+15
options template
title MERRA data
dtype netcdf
xdef  288 linear -179.375 1.25
ydef  144 linear -89.375 1.25
zdef   21 levels 1000 975 950 925 900 875 850 825 800 775 750 725 700 650 600 550 500 450 400 350 300
tdef  396 linear 00Z01JAN1979 1mn
vars 3
qv 21 t,z,y,x Specific humidity
u 21 t,z,y,x Eastward wind component
v 21 t,z,y,x Northward wind component
endvars

MacTaylor

易小凯 发表于 2012-5-6 15:59

                               
登录/注册后可看大图

写成%y4%m2.nc的形式就行

就是这样写的，能帮我看看哪错了吗？

mofangbao

打开后你q files看看

MacTaylor

mofangbao 发表于 2012-5-7 09:20

                               
登录/注册后可看大图

打开后你q files看看

能打开，q file也有用，能显示信息，但是就是画图d 变量的时候没值

易小凯

MacTaylor 发表于 2012-5-7 06:27

                               
登录/注册后可看大图

我之前就这样做的，能打开ctl文件，但是d之后，都显示all undefined values，我的ctl如下，麻烦帮我看看哪 ...

dset  ^%y4%m2.nc
换成绝对路径

MacTaylor

易小凯 发表于 2012-5-7 12:47

                               
登录/注册后可看大图

dset  ^%y4%m2.nc
换成绝对路径

换成绝对路径也试过，还是一样的问题。。。。。

易小凯

那个nc资料貌似不需要ctl文件，要连续打开可以用'sdfopen d:\data\ncep\vwnd.1999.nc d:\data\ncep\vwnd.%y4.nc 4383'这个格式
如果那个是二进制格式的，你试一下
DSET  D:\data\40yr_month_pressure\era40_monthly_pressure_%y4.dat
TITLE monthly ERA-40, from 9/1957 to 8/2002
UNDEF  9.999E+20
OPTIONS  SEQUENTIAL template big_endian
看看这种格式行不