求助：如何用fwrite提取wrfout中单个变量，指定层次上的数据？

∝冷暖如墨≈

我想从wrfout_d03.dat中提取850hPa上geopt的值，fwrite的gs和ctl等文件如下：
1.对wrfout_d03.ctl做了如下修改

dset e:\wrf\data\wrfout_d03.dat

options  byteswapped

undef 1.e30

title  OUTPUTFROM WRF V3.5.1 MODEL

*pdef  666 399 lcc 29.591  117.547  333.500 200.000  60.00000  30.00000 118.50000   1500.000  1500.000

*xdef 1602 linear  112.09885  0.00675676

*ydef  841 linear  26.67278   0.00675676

xdef 666 linear 112.09885   0.00675676

ydef 399 linear  26.67278   0.00675676

zdef   19levels  

1000.00000

950.00000

900.00000

850.00000

2.提取的gs文件：

'reinit'

'open e:\wrf\data\wrfout_d03.ctl'

'set fwrite e:\wrf\data\geopt\850.grd'

'set gxout fwrite'

'set x 1 666'

'set y 1 399'

'set t 1 31'

'set z 4'

'd geopt'

'disable fwrite'

;

3.对提取出的数据的ctl文件

dset E:\WRF\data\geopt\850.grd

undef 1.e30

title OUTPUT FROM WRF V3.5.1 MODEL

pdef  666 399lcc  29.591  117.547 333.500  200.000  60.00000 30.00000 118.50000   1500.000   1500.000

xdef 1602 linear 112.09885   0.00675676

ydef  841linear   26.67278   0.00675676

*xdef 666 linear 112.09885   0.00675676

*ydef 399 linear  26.67278   0.00675676

zdef 1 levels 850.000000

tdef 31 linear 18Z13JUN2011      60MN     

vars 1

geopt 1  0  Geopotential (m2/s2)

endvars


4.画出的图：




麻烦各位大神帮我看一下，问题出现在哪边？万分感激！

∝冷暖如墨≈

{:5_275:}{:5_275:}