- 积分
- 151
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2016-6-16
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
各位同仁,我想处理一批降雨数据197901-201512,格式为NC数据,单个文件我已经可以读出来了,现在关键是如果批处理,完成这444个文件的提取呢?写了一个CTL和GS文件,但是老是报错,请各位大神看看是什么意思?错在哪里?
ctl文件如下:
Dset F:/data/precipatation/prec_ITPCAS-CMFD_V0106_B-01_03hr_010deg_%y4%m2.nc
TITLE Data of precipation rate
OPTIONS template
UNDEF -9.99e+33
DTYPE netcdf
Xdef 700 Linear 70.05 0.1
Ydef 400 Linear 15.05 0.1
Tdef 248 Linear 00Z01JAN1979 180mn
Vars 1
Prec 0 -103 precipitation rate(mn)
endvaes
gs文件如下:
'c'
'reinit'
'xdfopen c:/GrADS20/data/test.ctl'
year=1979
month=1
t=1
while(year<=2015)
while(month<13)
if(month<10)
aaa=year%0%month *如果月小于10,则在月前加0,保证月是两位
else
aaa=year%month *如果月大于等于10,则原样输出,因为它已经是2位了
endif
'set lon 100 115'
'set lat 20 35'
'set t 't
'sstave=ave(prec,t=1,t=248)'
'enable print f:/data/prec'%aaa%'.gmf'
'print'
'disable print'
'c'
month=month+1
t=t+1
endwhile
month=1
year=year+1
endwhile
'reinit'
我写单个文件的时候,GS是这样的:
'reinit'
'sdfopen f:/data/precipatation/prec_ITPCAS-CMFD_V0106_B-01_03hr_010deg_197901.nc'
'set lon 100 115'
'set lat 20 35'
'sstave=ave(prec,t=1,t=248)'
'grads2ascii sstave c:/GrADS20/data/txt6.txt %g 1 1'
'grads2ascii lat c:/GrADS20/data/lat1.txt %g 1 1'
'grads2ascii lon c:/GrADS20/data/lon1.txt %g 1 1'
'reinit'
;
不懂如何把单个文件加入到批量处理中,红色的代码是我按照论坛的贴子来加上去的,用来做测试的,具体如何对每个NC文件进行批处理TXT文件,实在不懂,请各位大神指点指点。
|
|