- 积分
 - 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文件,实在不懂,请各位大神指点指点。 
 |   
 
 
 
 |