帮我看看这个gs怎么运行不成功？？

夜深沉

想求不同年份不同日期的高度平均场，可是运行总出错，附我的gs和出错信息，求大侠指点一二。
'sdfopen e:\hgt.1983.nc'
'sdfopen e:\hgt.1988.nc'
'enable print e:\500hgt.gmf'
'set gxout contour'
'set cterp on'
'set csmooth on'
'set lon 20 120'
'set lat 20 90'
'set z 6'
'define aveH1=ave(hgt.1,t=735,t=738)'
'define aveH2=ave(hgt.2,t=801,t=804)'
'define anom=(aveH1+aveH2)/2'
'd anom/10'
'print'
'disable print'

愉悦

时间问题
grads没办法统一那个t，你把其中一个年份的转化成dat，然后写CTL的时候把时间和另一个年份写成一样就可以了

夜深沉

没明白意思喔~~~

mofangbao

define之前加上set t 1

牧场物语之主

我觉得也是时次的问题吧，grads不知道你   'define aveH1=ave(hgt.1,t=735,t=738)'     和    'define aveH2=ave(hgt.2,t=801,t=804)'     分别处理的是哪个文件。

river

define语句之前加上set t 1
里面的1其实可以是随便一个固定的数字的，就是为了表示平均完只剩一个时次了

夜深沉

mofangbao 发表于 2013-10-30 12:39

                               
登录/注册后可看大图

define之前加上set t 1

'sdfopen e:\hgt.1983.nc'
'sdfopen e:\hgt.1988.nc'
'enable print e:\500hgt.gmf'
'set gxout contour'
'set cterp on'
'set mpdset cnworld'
'set csmooth on'
'set lon 20 120'
'set lat 20 90'
'set z 6'
'set t 1'
'define aveH1=ave(hgt.1,t=735,t=800)'
'define aveH2=ave(hgt.2,t=801,t=824)'
'define anom=(aveH1+aveH2)/2'
'd anom/10'
'print'
'disable print'
结果提示还是跟以前一样啊~~~