2种方法的资料合并出现错误

追风的阳光

我分别试了俩种方法试图将多个时次的hdf资料提取某些变量并合并，俩种方法都出现了错误
1
第一种我逐个打开文件，试图合并，grads显示下图错误，还能形成数据，但数据已发生改变

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''
tt=1
while(tt<=48)
'set t 'tt
zz=1
while(zz<=42)
'set z 'zz
'set x 108.125 130.625'
'set y 24.375 40.625'
'd U'
zz=zz+1
endwhile
tt=tt+1
endwhile
i=i+1
endwhile
'disable fwrite'
'reinit'
;
2
第二种我编写了一个批量描述ctl，来实现，grads出现下图错误


ctl：
dset d:/zssdata/merra2/MERRA300.prod.assim.inst3_3d_asm_Cp.201107%d2.hdf
options template
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


gs：
'reinit'
'open d:/zssdata/merra2/20110710test2.ctl'
'set fwrite d:/zssdata/merra2/U.dat'
'set gxout fwrite'
tt=1
while(tt<=48)
'set t 'tt
zz=1
while(zz<=42)
'set z 'zz
'set lon 108.125 130.625'
'set lat 24.375 40.625'
'd U'
zz=zz+1
endwhile
tt=tt+1
endwhile
'disable fwrite'
;

请问这俩个方法都错在了什么地方。。

river

先看一下基础知识吧

nunu18

grands 直接可以打开.hdf格式的文件吗？

一水天下

'set x 108.125 130.625'
'set y 24.375 40.625'里面应该改成lon，lat吧

追风的阳光

river 发表于 2014-3-21 17:32

                               
登录/注册后可看大图

先看一下基础知识吧

请问到底是哪儿有问题呢。。。还是没发现。。

阿木木

{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}{:eb502:}

river

追风的阳光 发表于 2014-3-24 20:55
请问到底是哪儿有问题呢。。。还是没发现。。

基础的错误，看一下手册就知道了

追风的阳光

river 发表于 2014-3-25 17:11
基础的错误，看一下手册就知道了

求问。。到底是哪儿的问题啊。。搞几天了。。

river

追风的阳光 发表于 2014-3-27 11:53
求问。。到底是哪儿的问题啊。。搞几天了。。

地板的回答，您老人家都看不见吗。

追风的阳光

river 发表于 2014-3-27 15:36
地板的回答，您老人家都看不见吗。

不是lon，lat的问题啊。。。原先写的lon，lat，后来试着写的x，y，还是一样的错误。。