求助！处理区域平均的脚本变量问题

呆萌萱

我先用月平均资料计算生成年平均资料，脚本如下：

'reinit'
'sdfopen F:\benke\xiaobo\1979_2015_3.nc'
'set fwrite F:\benke\xiaobo\100ave.grd'
'set gxout fwrite'
'set lev 100'

it=13
'set t 'it''
while(it<=444)
'define tave=ave(t,t='it',t='it+11')'
'd tave'
it=it+12
endwhile


'disable fwrite'
'reinit'
;


之后我编写了100ave.grd的描述文件，脚本如下：

dset F:\benke\xiaobo\100ave.grd
title ave of year
undef -9.99e+33
*dtype netcdf
xdef 144 linear 0 2.5
ydef 73  linear -90 2.5
zdef 1   levels 100
tdef 36 linear JAN1980 1yr
vars 1
tave 36  0  Temperature
endvars



我再用gs算区域平均，脚本如下：

'reinit'
'open F:\benke\xiaobo\100ave.ctl'
'set fwrite F:\benke\xiaobo\100quyu.grd'
'set gxout fwrite'
it=1
while(it<=36)
'set t 'it''
'define  a=aave(tave,lon=75,lon=105,lat=25,lat=40)'
'd tloop(a)'
it=it+1
endwhile

'disable fwrite'
'print'
'reinit'

；


但是最后出现了如图的错误，这怎么解决啊？？？

呆萌萱

真的好困惑啊 为什么一直变量有问题

liyf

虽说不忍吐槽了，但还是提醒一下，自己再看看fwrite、ctl如何编怎么写，

liyf

才发现，怎么又是你

278803532

，错了就自己Open以下ctl, d一下变量啊，自己找错啊。。。。

呆萌萱

liyf 发表于 2016-4-7 22:23
才发现，怎么又是你

我是个初学者，谢谢指点

呆萌萱

liyf 发表于 2016-4-7 22:23
才发现，怎么又是你

我是个初学者，谢谢指点

呆萌萱

278803532 发表于 2016-4-8 08:27
，错了就自己Open以下ctl, d一下变量啊，自己找错啊。。。。

恩恩  我是初学者 没经验