求教：用fwrite 写多层次数据文件不能画图

nancy

各位： 问题如主题所示，使用era-40月平均多层grib文件
GS 文件如下：

'reinit'
'open E:\GrADS19\tropopause\t.ctl'

'set lon 72.5 135'
'set lat 15 55'
'set t 1 528'

'set fwrite E:\GrADS19\tropopause\temp.grd'
'set gxout fwrite'


i=1
while(i<=9)
say i
'set z 'i
'd tprs'
i=i+1
endwhile
'disable fwrite'

'reinit'
'open E:\GrADS19\tropopause\temp.ctl'
'enable print E:\GrADS19\tropopause\500.gmf'
'set grads off'
'set grid off'
'set lev 500'
'set t 1'
'set lon 72.5 135'
'set lat 15 55'
'd tprs'
'print'
'disable print'
;



CTL文件如下：

dset E:\GrADS19\tropopause\temp.grd
title Temperature
dtype grib
*options yrev
undef 9.999E+20
ydef 17 linear 15 2.5
xdef 26 linear 72.5 2.5
tdef 528 linear 00Z01jan1958 1mo
*  z has 9 levels, for prs
zdef 9 levels
500 400 300 250 200 150 100 70 50
vars 1
Tprs 9 99 Temperature
ENDVARS


不知怎么回事，请求您的热心帮助！谢谢！

传说中的谁

请给出你的出错信息，大家可以一起探讨。

mofangbao

恩 给出错误截图或者提示 最好给你的gs加上注释和你对该程序的理解，这样便于帮您解决问题，大家一起交流进步，

凡不凡

正在学GRADS