fwrite命令截取数据出现变量错位

swallow19810923

各位大侠，本人在使用fwrite命令截取数据的时候出现变量错位（比如层次错位）的情况，请高手们给予指点，不胜感激：）
gs文件：
'reinit'
'open F:\ncep\201001\201001.ctl'
'set mpdset cnworld cnriver'
'set lon 80 130'
'set lat 10 50'
'set fwrite F:\dat\201001.dat'
'set gxout fwrite'
lev.1=850;lev.2=700;lev.3=500;lev.4=200
t=1
while(t<=124)
'set t 't
l=1
while(l<=4)
lev=lev.l
'set lev 'lev
'd hgtprs'
'd RHprs'
'd TMPprs'
'd UGRDprs'
'd VGRDprs'
'd VVELprs'
l=l+1
endwhile
'd HGTsfc'
'd PRESsfc'
'd RH2m'
'd SPFH2m'
'd TMPsfc'
'd TMP2m'
t=t+1
endwhile
'disable fwrite'
'c'

ctl文件：
dset F:/dat/201001.dat
undef 9.999E+20
title 201001grib data
xdef 51 linear 80.000000 1.000000
ydef 41 linear 10.000000 1.000000
tdef 124 linear 00Z01jan2010 6hr
zdef 4 levels 850 700 500 200
vars 12
HGTsfc  0 99 surface Geopotential height [gpm]
HGTprs 4 99 Geopotential height [gpm]
PRESsfc  0 99 surface Pressure [Pa]
RHprs 4 99 Relative humidity [%]
RH2m  0 99 2 m Relative humidity [%]
SPFH2m  0 99 2 m Specific humidity [kg/kg]
TMPsfc  0 99 surface Temp. [K]
TMPprs 4 99 Temp. [K]
TMP2m  0 99 2 m Temp. [K]
UGRDprs 4 99 u wind [m/s]
VGRDprs 4 99 v wind [m/s]
VVELprs 4 99 Pressure vertical velocity [Pa/s]
ENDVARS

男紫汗

没取过这么多而复杂的数据。。。变量可以同时d吗

river

你仔细研究下grads要求的数据排列方式，再看看你的gs就知道错在哪了

swallow19810923

谢谢“男紫汗”大侠， 我试过只截取一个多层的，一个单层的，是可以的。d的变量一多久错位了

swallow19810923

谢谢“river”大侠， 我再看看吧：）

jolincai

变量太多  只做过简单的 循环对吗 lev.1 lev.2这种表达对吗..

sosbbxc

3L说的对，你仔细看一下grads里面对维数排列顺序的要求，从内到外依次是经度，纬度，高度，变量和时间，时间维在最外面，如果按照你这种d法，变量这一维就跑到最外面那一维了，那d出来的数据肯定全乱了。所以我一般在提取有好多变量的多维数组的时候，要么用fortran，要么每次只d一个变量，需要几个变量就生成几个dat文件

18393810639

男紫汗 发表于 2016-10-11 10:14
没取过这么多而复杂的数据。。。变量可以同时d吗

我是特别想把这个数据中的不同变量写进不同的文件里。
大神有什么好办法吗？不胜感激啊。

男紫汗

18393810639 发表于 2017-5-9 08:47
我是特别想把这个数据中的不同变量写进不同的文件里。
大神有什么好办法吗？不胜感激啊。

这个似乎没有简单的方法，老老实实写