grads怎样读取四五个ctl数据，然后显示出这些数据中其中一个变量的平均值

wuyang0074

大家好：
我想问一下，grads怎样读取四五个ctl数据，然后显示出这些数据中其中一个变量的平均值？

谢谢各位！

柚子皮

grads可以同时读取多个文件，所用文件的变量后面缀上文件号就可以了，型如: define H=(hgt.1+hgt.2+hgt.3+...)/n。其中hgt.1的.1可以省略。大概就是这样了

wuyang0074

柚子皮 发表于 2017-10-15 08:44
grads可以同时读取多个文件，所用文件的变量后面缀上文件号就可以了，型如: define H=(hgt.1+hgt.2+hgt.3+. ...

您能帮我看一下吗？这是哪里出错了？
'clear'
'reinit'

'open 2000_07.ctl'
'open 2005_07.ctl'


'set lon 86 127'
'set lat 37 54'
'set grads off'
'set timelab off'
'set grid off'
'set lev 1'
'set mpdraw off'

'set xlint 10'
'set ylint 5'

'set gxout shaded'
'set csmooth on'
'set clevs 0.01 0.05 0.09 0.13 0.17 0.21 0.25 0.29 0.33 0.37 0.41 0.45 0.49 0.53 0.57 0.61 0.65 0.69 0.74 0.8 0.9 1'

'define a=ave(smois,t=1,t=117)'
'define b=ave(smois.2,t=1,t=117)'
'd (a+b)/2'

'draw shp MGGY.shp'

'set strsiz 0.2 0.2'
'draw string 0.7 7 (2000-07)'
'cbarn 0.7 0 5.5 1.2'

柚子皮

'define a=ave(smois,t=1,t=117)'  定义之前，一般对lon,lat,lev,time作规定，define之前补充一个'set t 1'

柚子皮

柚子皮 发表于 2017-10-15 12:51
'define a=ave(smois,t=1,t=117)'  定义之前，一般对lon,lat,lev,time作规定，define之前补充一个'set t 1'

提问的时候，最好附上你的程序和错误提示，方便大家帮你分析

wuyang0074

柚子皮 发表于 2017-10-15 12:52
提问的时候，最好附上你的程序和错误提示，方便大家帮你分析

谢谢您，我改了，可是还是不行，您能不能贴出一段您读取多个数据的小程序，我对照着学习一下呀。这样能快一些，谢谢您了。