gs运行错误

爱高斯拉普拉斯

各位大神、前辈，我刚接触GRADS，我想实现一年数据（grib1格式）的批量读取，按照兰溪的神贴生成了ctl和对应的idx，也都匹配成功了，我编写了一个gs，运行提示第二行有误（如上图），第二行是读取ctl，我单独读取又可以，不知道问题出在哪儿。下面附上gs和ctl，烦请大神们帮忙看看

爱高斯拉普拉斯

gs:
reinit
open E:/sj/grib1/2000/fnl_20000101_00.ctl
set gxout fwrite
set fwrite E:/sj/grib1/2000/fnl_20000101_00.dat
set lat 32 42
set lon 110 120
t=1
while(t<=1460)
set t
z=1
while(z<=26)
setz z
HPBLsfc
z=z+1
endwhile
z=1
while(z<=26)
setz z
RHprs
z=z+1
endwhile
t=t+1
endwhile
disable fwrite

爱高斯拉普拉斯

ctl:
dset E:/sj/grib1/2000/fnl_2000%m2%d2_%h2_00.grib1.spasub.dong152823
options template
index E:/sj/grib1/2000/fnl_20000101_00_00.grib1.spasub.dong152823.idx
undef 9.999E+20
title E:/sj/grib1/2000/fnl_20000101_00_00.grib1.spasub.dong152823
*  produced by grib2ctl v0.9.12.5p33k
dtype grib 255
options yrev
ydef 16 linear 30.000000 1
xdef 21 linear 105.000000 1.000000
tdef 1460 linear 00Z01jan2000 6hr
zdef 26 levels
1000 975 950 925 900 850 800 750 700 650 600 550 500 450 400 350 300 250 200 150 100 70 50 30 20 10
vars 33
HPBLsfc  0 221,1,0  ** surface Planetary boundary layer height [m]
RHprs 21 52,100,0 ** Relative humidity [%]
RH2m  0 52,105,2 ** 2 m Relative humidity [%]
RHsig995  0 52,107,9950 ** sig=.995 Relative humidity [%]
RHsg44_100  0 52,108,11364 ** sigma=0.44-1.00 layer Relative humidity [%]
RHsg72_94  0 52,108,18526 ** sigma=0.72-0.94 layer Relative humidity [%]
RHsg44_72  0 52,108,11336 ** sigma=0.44-0.72 layer Relative humidity [%]
RHsg33_100  0 52,108,8548 ** sigma=0.33-1.00 layer Relative humidity [%]
RH30_0mb  0 52,116,7680 ** 30-0 mb above gnd Relative humidity [%]
RHclm  0 52,200,0 ** atmos column Relative humidity [%]
TMPsfc  0 11,1,0  ** surface Temp. [K]
TMPprs 26 11,100,0 ** Temp. [K]
TMP1829m  0 11,103,1829 ** 1829 m Temp. [K]
TMP2743m  0 11,103,2743 ** 2743 m Temp. [K]
TMP3658m  0 11,103,3658 ** 3658 m Temp. [K]
TMP2m  0 11,105,2 ** 2 m Temp. [K]
TMPsig995  0 11,107,9950 ** sig=.995 Temp. [K]
TMPSoilT  0 11,112,10 ** 0-10 cm undergnd Temp. [K]
TMPSoilM  0 11,112,2760 ** 10-200 cm undergnd Temp. [K]
TMP30_0mb  0 11,116,7680 ** 30-0 mb above gnd Temp. [K]
TMPmwl  0 11,6,0 ** max wind level Temp. [K]
TMPtrp  0 11,7,0 ** tropopause Temp. [K]
VGRDprs 26 34,100,0 ** v wind [m/s]
VGRD1829m  0 34,103,1829 ** 1829 m v wind [m/s]
VGRD2743m  0 34,103,2743 ** 2743 m v wind [m/s]
VGRD3658m  0 34,103,3658 ** 3658 m v wind [m/s]
VGRD10m  0 34,105,10 ** 10 m v wind [m/s]
VGRDsig995  0 34,107,9950 ** sig=.995 v wind [m/s]
VGRD30_0mb  0 34,116,7680 ** 30-0 mb above gnd v wind [m/s]
VGRDmwl  0 34,6,0 ** max wind level v wind [m/s]
VGRDtrp  0 34,7,0 ** tropopause v wind [m/s]
VVELprs 21 39,100,0 ** Pressure vertical velocity [Pa/s]
VVELsig995  0 39,107,9950 ** sig=.995 Pressure vertical velocity [Pa/s]
ENDVARS

fresh_fish

{:5_235:}