批处理生成ctl后，提取数据时绝大部分数值显示缺测

气象学渣

我正在用wrf模式的预报资料，数据格式是grib2，我现在批量处理1-10月份每天的预报时效为24小时的数据，已经生成了描述文件ctl和映射文件idx，我想要提取出每天的温度，但写好gs运行之后几乎都是缺测值-9.99e+08这样的，只有五六个正常数值，但不做批处理单独提取某一天的温度却可以显示正常值。请问各位大神该怎么解决？

附ctl：
dset D:\fjRAP_wrf_d01_2017%m2%d200f024.grib2
index D:\fjRAP_wrf_d01_2017010100f024.grib2.idx
undef 9.999E+20
title D:\fjRAP_wrf_d01_2017%m2%d200f024.grib2
* produced by g2ctl v0.1.1
* command line options: D:\20180115Dates\fj\201708\fjRAP_wrf_d01_2017%m2%d200f024.grib2
* griddef=1:0:(261 x 261):grid_template=0:winds(N/S): lat-lon grid:(261 x 261) units 1e-06 input WE:SN output WE:SN res 48 lat 13.000000 to 36.400000 by 0.090000 lon 107.600000 to 131.000000 by 0.090000 #points=68121:winds(N/S)

dtype grib2
options template
ydef 261 linear 13.000000 0.09
xdef 261 linear 107.600000 0.090000
tdef 304 linear 00Z02jan2017 1dy
* PROFILE hPa
zdef 24 levels 100000 97500 95000 92500 90000 87500 85000 82500 80000 75000 70000 65000 60000 55000 50000 45000 40000 35000 30000 25000 20000 15000 10000 200
options pascals
vars 131
（变量太多就略去了，我想要的变量为TMPsfc）
ENDVARS


附gs：
'reinit'
'open D:\fjRAP_wrf_d01_2017.ctl'
n=1
while(n<=304)
'set lon 118.04'
'set lat 24.25 24.34'
'set t 'n''
'define a=ave(TMPsfc,lat=24.25,lat=24.34)'
'd a-273.15'
out=subwrd(result,4)
say out
n=n+1
endwhile
;


运行结果（部分截图）：



求各位指点！！！