爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 4566|回复: 2

[分享资料] 求助,批处理.NC文件转为TXT文件

[复制链接]
发表于 2016-9-21 22:59:29 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册

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文件,实在不懂,请各位大神指点指点。
密码修改失败请联系微信:mofangbao
发表于 2016-9-22 09:59:52 | 显示全部楼层
我只会用Matlab转换
密码修改失败请联系微信:mofangbao
发表于 2018-7-17 21:30:10 | 显示全部楼层
15754881543 发表于 2016-9-22 09:59
我只会用Matlab转换

可以教教我吗?
密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表