GRADS怎么处理多个NC文件？谢谢

duan2008

'reinit'
'sdfopen D:\TDDOWNLOAD\hgt.1981.nc'
'set dfile 1'
'set lat 0 60'
'set lon 70 180'
'set lev 500'
'set t 1'
'sdfopen D:\TDDOWNLOAD\hgt.1982.nc'
'set dfile 2'
'set lat 0 60'
'set lon 70 180'
'set lev 500'
'set t 1'
'define hgt12=hgt.1+hgt.2'
'set t 1'
'd hgt12'
'enable print D:\TDDOWNLOAD\1981-1982HGT.gmf'
请问高手  我的程序那里错误了？我想合并这两年的HGT。程序无法求取两年的和，有可能是第一年的时间抽和第二年的不一样，谢谢

river

'reinit'
'sdfopen D:\TDDOWNLOAD\hgt.1981.nc'
'sdfopen D:\TDDOWNLOAD\hgt.1982.nc'
'set lat 0 60'
'set lon 70 180'
'set lev 500'
'set t 1'
'define hgt12=hgt.1+hgt.2'
'd hgt12'
'enable print D:\TDDOWNLOAD\1981-1982HGT.gmf'
'print'
'disable print'
;
还不行就把错误提示发上来啊，这样才能对症下药啊

duan2008

river 发表于 2012-11-30 20:11

                               
登录/注册后可看大图

'reinit'
'sdfopen D:\TDDOWNLOAD\hgt.1981.nc'
'sdfopen D:\TDDOWNLOAD\hgt.1982.nc'

谢谢啊，还是不行。错误如下


我在想是时间的设置问题吗？因为第一个T=1是1981年1月1日，而第二个文件是1982年1月1日，而且有可能是365天和366天，在求和的时候有数据缺失？ 谢谢啊。

river

duan2008 发表于 2012-11-30 20:30

                               
登录/注册后可看大图

谢谢啊，还是不行。错误如下

把那句set t 1去掉，把define那句改成define hgt12=hgt.1(t=1)+hgt.2(t=1)试试

duan2008

river 发表于 2012-11-30 21:18

                               
登录/注册后可看大图

把那句set t 1去掉，把define那句改成define hgt12=hgt.1(t=1)+hgt.2(t=1)试试

这样可以计算出一个值，但是对于366天的值怎么计算了，把hgt12存入一个NC文件，但是T是366天的，谢谢。

river

duan2008 发表于 2012-11-30 21:57

                               
登录/注册后可看大图

这样可以计算出一个值，但是对于366天的值怎么计算了，把hgt12存入一个NC文件，但是T是366天的，谢谢。

用循环可以吧，加几句
'reinit'
'sdfopen D:\TDDOWNLOAD\hgt.1981.nc'
'sdfopen D:\TDDOWNLOAD\hgt.1982.nc'
'enable print D:\TDDOWNLOAD\1981-1982HGT.gmf'
'set lat 0 60'
'set lon 70 180'
'set lev 500'
it=1
while（it<=365）
'set t 'it''
'define hgt12=hgt.1(t='it')+hgt.2(t='it')
'draw title 'it''
'print'
'c'
it=it+1
'disable print'
;
但是多出来一天就不好办了啊

xjliutao

学习了，我一个都没有处理好 你都可以处理多个了

linikiku

支持一下！！！！

djcj520

river 发表于 2012-11-30 20:11
'reinit'
'sdfopen D:\TDDOWNLOAD\hgt.1981.nc'
'sdfopen D:\TDDOWNLOAD\hgt.1982.nc'

我也遇到过

yaotianci

谢谢，正好在学习grads，受益匪浅。