爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

新浪微博登陆

只需一步, 快速开始

搜索
查看: 3704|回复: 0

[分享资料] 关于hdfsds资料的问题

[复制链接]

新浪微博达人勋

发表于 2014-3-18 11:30:04 | 显示全部楼层 |阅读模式

登录后查看更多精彩内容~

您需要 登录 才可以下载或查看,没有帐号?立即注册 新浪微博登陆

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打开的还是不一样,请问这是什么原因呢?

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册 新浪微博登陆

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表