fwrite提取数据遇到问题，求教

kelly

wrf数据的ctl
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  
1000.00000
975.00000
950.00000
925.00000
900.00000
850.00000
800.00000
750.00000
700.00000
650.00000
600.00000
550.00000
500.00000
400.00000
300.00000
200.00000
100.00000
tdef   73 linear 12Z17JUN2012      60MN      
VARS   26
U             17  0  x-wind component (m s-1)
V             17  0  y-wind component (m s-1)
W             17  0  z-wind component (m s-1)
Q2             1  0  QV at 2 M (kg kg-1)
T2             1  0  TEMP at 2 M (K)
QVAPOR        17  0  Water vapor mixing ratio (kg kg-1)
HGT            1  0  Terrain Height (m)
TSK            1  0  SURFACE SKIN TEMPERATURE (K)
RAINC          1  0  ACCUMULATED TOTAL CUMULUS PRECIPITATION (mm)
RAINNC         1  0  ACCUMULATED TOTAL GRID SCALE PRECIPITATION (mm)
XLAT           1  0  LATITUDE, SOUTH IS NEGATIVE (degree_north)
XLONG          1  0  LONGITUDE, WEST IS NEGATIVE (degree_east)
height        17  0  Model height (km)
tc            17  0  Temperature (C)
theta         17  0  Potential Temperature (K)
td            17  0  Dewpoint Temperature (C)
rh            17  0  Relative Humidity (%)
clflo          1  0  Low Cloud Fraction (%)
clfmi          1  0  Mid Cloud Fraction (%)
clfhi          1  0  High Cloud Fraction (%)
rh2            1  0  Relative Humidity at 2m (%)
u10m           1  0  Rotated wind component (m s-1)
v10m           1  0  Rotated wind component (m s-1)
slp            1  0  Sea Levelp Pressure (hPa)
dbz           17  0  Reflectivity (-)
cape          17  0  CAPE (J/kg)
ENDVARS


我想提取850这层的各时次的温度（tc）并写入一个文件中，gs如下：
'reinit'
'open D:\sucai\wrf\wrfout.ctl'
****************shuchu T/TH/WIND
'set fwrite  D:\sucai\wrf\wrf_1.grd'
'set gxout fwrite'
i=1
while(i<=73)
'set t 'i
'set z 13'
'set lon 112.06458 122.10512'
'set lat 28.11144 36.28711'
'd tc'
i=i+1
endwhile
'disable fwrite'
*'quit'
;

我对提取出的数据wrf_1.grd 配上了ctl
dset D:\sucai\wrf\wrf_1.grd  
undef 1.e30
title  OUTPUT FROM WRF V3.4 MODEL
xdef  744 linear  112.06458   0.01351351
ydef  606 linear   28.11144   0.01351351
zdef   1 levels  1000
tdef   73 linear 12Z17JUN2012      60MN      
VARS   1
tc            0  99  Temperature (C)
ENDVARS

画图gs
'reinit'
'open D:\sucai\wrf\wrf_1.ctl'
'set t '1
'set grads off'
'set grid off'
'set frame off'
'set lon 114.7 119.8'
'set lat 29.3 34.8'
'set xlopts  0 0 0'
'set ylopts  0 0 0'
'set gxout contour'
'set cint 0.1'
'd tc'
'printim D:\sucai\wrf\shishi1.png white'
'c'
*'quit'
;

kelly

请问图上缺的是咋回事？  t=1,2,3缺的地方还不同呢

cassienju

还有个pdef，是不是写出数据的时候忽略了？

aa22556688

导出数据的时候，不要用 set lon 和 set lat，容易发生数据漂移，最好用 set x 和 set y，

风儿飘飘

提取850hPa的数值，那么设置z 时应为set z 6或set lev 850

风儿飘飘

提取850hPa的数值，那么设置z 时应为set z 6或set lev 850

18393810639

我是特别想把这个数据中的不同变量写进不同的文件里。