- 积分
- 1831
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-5-6
- 最后登录
- 1970-1-1
|
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
有多个时次的hdf资料,我想将其合并成一个文件并提取其中变量,
下面是我的gs
'reinit'
'sdfopen d:/zssdata/merra2/MERRA300.prod.assim.inst3_3d_asm_Cp.20110710.hdf'
'sdfopen d:/zssdata/merra2/MERRA300.prod.assim.inst3_3d_asm_Cp.20110711.hdf'
'sdfopen d:/zssdata/merra2/MERRA300.prod.assim.inst3_3d_asm_Cp.20110712.hdf'
'sdfopen d:/zssdata/merra2/MERRA300.prod.assim.inst3_3d_asm_Cp.20110713.hdf'
'sdfopen d:/zssdata/merra2/MERRA300.prod.assim.inst3_3d_asm_Cp.20110714.hdf'
'sdfopen d:/zssdata/merra2/MERRA300.prod.assim.inst3_3d_asm_Cp.20110715.hdf'
'set fwrite d:/zssdata/merra2/U.grd'
'set gxout fwrite'
i=1
while(i<=6)
'set dfile 'i''
'set lon 108.125 130.625'
'set lat 24.375 40.625'
tt=1
while(tt<=48)
'set t 'tt
zz=1
while(zz<=42)
'set z 'zz
'd U'
zz=zz+1
endwhile
tt=tt+1
endwhile
i=i+1
endwhile
'disable fwrite'
'reinit'
;
上面是我提取合并的gs,显示‘request is completely outside limits’错误,并且提取出的数据数值有问题,与sdfopen打开的不同,gs中有什么错误吗?
上一个问题没有解决,我写了ctl文件,试图将多个文件合并:
dset d:/zssdata/merra2/MERRA300.prod.assim.inst3_3d_asm_Cp.201107%d2.hdf
title MERRA reanalysis. GEOS-5.2.0
dtype hdfsds
undef 1e+15
xdef 288 linear -179.375 1.25
ydef 144 linear -89.375 1.25
zdef 42 levels 1000 975 950 925 900 875 850 825 800 775 750 725 700 650 600 550
500 450 400 350 300 250 200 150 100 70 50 40 30 20 10 7 5 4 3 2 1 0.7 0.5 0.4 0.3
0.1
tdef 48 linear 00Z10Jul2011 180mn
vars 18
SLP=>slp 0 t,y,x Sea-level pressure
PS=>ps 0 t,y,x Surface pressure
PHIS=>phis 0 t,y,x Surface Geopotential
H=>h 42 t,z,y,x Geopotential height
O3=>o3 42 t,z,y,x Ozone Mixing Ratio
QV=>qv 42 t,z,y,x Specific humidity
QL=>ql 42 t,z,y,x Cloud liquid water mixing ratio
QI=>qi 42 t,z,y,x Cloud ice mixing ratio
RH=>rh 42 t,z,y,x Ralative humidity
T=>t 42 t,z,y,x Air Temperature
U=>u 42 t,z,y,x Eastward wind component
V=>v 42 t,z,y,x Northward wind component
EPV=>epv 42 t,z,y,x Ertel potential vorticity
OMEGA=>omega 42 t,z,y,x Vertical pressure velocity
XDim=>xdim 0 x XDim
YDim=>ydim 0 y YDim
Height=>height 42 z Height
Time=>time 0 t Time
endvars
ctl可以打开,但发现ctl打开后的同一时刻,同一高度的变量值与sdfopen打开的还是不一样,请问这是什么原因呢?
|
|