（已解决）关于将多层数据写成一个文件的求助

MacTaylor · 发表于 2012-7-20 23:53:52

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

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由 MacTaylor 于 2012-7-21 00:05 编辑

处理JRA25的gr数据，想把hgt多层数据写到一个文件里面，gs脚本如下：

'reinit'
'open F:\JRA25\data\anl_p25.monthly.ctl'
'set gxout fwrite'
'set fwrite F:\JRA25\hgt.bin'
i=1
while(i<=396)
j=1
  while(j<=23)
'set x 1 144'
'set y 1 73'
'set t 'i''
'set z 'j''
'd HGTprs'
j=j+1
endwhile
i=i+1
endwhile
'reinit'

生成文件的配套ctl文件如下：
dset  F:\JRA25\hgt.bin
title  Geopotential height
undef 9.999E+20
ydef 73  linear -90.0 2.5
xdef  144  linear 0.0 2.5
tdef 396  linear  01jan1979  1mo
zdef 23  levels
1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10 7 5 3 2 1 0.4
tdef  396 linear 01jan1979 1mo
vars 1
hgt  1  99  Geopotential height          [gpm]
endvars

但是用grads打开后，只能画出第一层的，其他层都不行，错误提示：
Data request warning: request is completely outside file limits
  Entire grid contents are set to missing data
  Grid limits of file:             X=1 144  Y=1 73  Z=1 23  T=1 396  E=1 1
  Grid limits of request:    X=1 144  Y=1 73  Z=1 23  T=1 396  E=1 1
  Warning issued for variable = hgt
Cannot contour grid - all undefined values

应该是加的那个层数循环有问题，请教大家应该怎么写啊？

MacTaylor · 发表于 2012-7-21 00:06:16

搞定了，一个低级错误。。。

royalosyin · 发表于 2012-7-21 04:22:21

怎么搞定的？顺便把修正的地方也贴出来貝

MacTaylor · 发表于 2012-7-21 04:59:32

royalosyin 发表于 2012-7-21 04:22

登录/注册后可看大图

怎么搞定的？顺便把修正的地方也贴出来貝

很低级的错误。。。

hgt 1 99 Geopotential height [gpm]改成
hgt 23 99 Geopotential height [gpm]

blui · 发表于 2013-4-8 20:30:15

为什么要改成23呢？

MacTaylor · 发表于 2013-4-9 03:57:38

blui 发表于 2013-4-8 20:30

登录/注册后可看大图

为什么要改成23呢？

因为是23层

methouse · 发表于 2017-2-17 10:09:52

谢谢分享，学习学习。

cibaouc · 发表于 2017-11-17 19:11:44

感谢楼主的分享，遇到同样的情况，顺利解决

		自动登录	找回密码
密码			立即注册

[分享资料] （已解决）关于将多层数据写成一个文件的求助

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