怎样从wrf模式预报结果中用fwrite提取所需数据

kelly

dset D:\sucai\wrf\wrf03d2014092512.dat
options big_endian   
undef 1.e30
title  OUTPUT FROM WRF V3.4 MODEL
pdef  288 288 lcc  32.236  117.107  144.500  144.500  60.00000  30.00000  117.30000   3000.000   3000.000
xdef  744 linear  112.06458   0.01351351
ydef  606 linear   28.11144   0.01351351
zdef   17 levels  
……
这样出来的数据有问题啊，画的图跟从原数据画出来的就是不一样
'reinit'
'open D:\sucai\wrf\wrfout.ctl'
'set fwrite  D:\sucai\wrf\wrf_1.grd'
'set gxout fwrite'
i=1
while(i<=73)
'set t 'i
'set z 13'
'set lon 114.7 119.8'
'set lat 29.3 34.8'
'd tc'
*'d rh'
*'d u'
*'d v'
i=i+1
endwhile
'disable fwrite'
*'quit'