用grads显示.grd格式的ctl文件的问题

东篱

用NCEP高度场1948-2013月平均资料分析青藏高原夏季风指数IPM，首先用一个 hgt-index.gs (479 Bytes, 下载次数: 10)

2013-5-12 22:40 上传

点击文件名下载附件 将NCEP资料hgt.mon.mean（文件太大没法上传）转化为.grd 格式文件 hgt_index.grd (3.05 KB, 下载次数: 6)

2013-5-12 22:42 上传

点击文件名下载附件 。再为.grd编写ctl文件为
dset f:\data\hgt_index.grd
title 1948--3013-PMI
undef -2.56e33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 linear 600 1
tdef 780 linear 00Z01jan1948 1mo

在grads中运行错误
求各位大神帮帮忙，检查下ctl文件有没什么问题，先在这里谢谢啦！

补充：
.gs文件
'reinit'
'sdfopen f:\data\hgt.mon.mean.nc'
'set gxout fwrite'
'set fwrite f:\data\tp\hgt-index.grd'
*'set t 1 780'
'set lev 600'

tt=1
while(tt<=780)
'set t 'tt''
'define a1=aave(hgt,x=33,x=33,y=50,y=50)'
'define a2=aave(hgt,x=37,x=37,y=47,y=47)'
'define a3=aave(hgt,x=41,x=41,y=50,y=50)'
'define a4=aave(hgt,x=37,x=37,y=53,y=53)'
'define a5=aave(hgt,x=37,x=37,y=50,y=50)'
'define temp=a1+a2+a3+a4-4*a5'
'd temp'

tt=tt+1
endwhile
'disable fwrite'
;

高原季风指数采用汤懋苍定义的指数，取夏季青藏高原四周的(32.5°N, 80°E)、(25°N, 90°E)、(32.5°N，100°E) 、(40°N, 90°E)四点和中心点 (32.5°N ，90°E) 共5 个点600 hPa高度值H1 、H2 、H3 、H4、H0 ，计算H1 ＋H2 ＋H3＋H4 - 4H0 ，取其月距平值作为反映高原夏季风强弱的高原夏季风指数IPM

lqouc

把ctl贴完整，再把gs直接贴上来。还有你的报错截图。
就只是你说的这些信息，没法帮你。

东篱

lqouc 发表于 2013-5-12 22:56

                               
登录/注册后可看大图

把ctl贴完整，再把gs直接贴上来。还有你的报错截图。
就只是你说的这些信息，没法帮你。

刚注册，不知道下载附件要减分，谢谢回复。

lqouc

东篱 发表于 2013-5-12 23:09

                               
登录/注册后可看大图

刚注册，不知道下载附件要减分，谢谢回复。

不时间不减分的问题，我的积分绝对够用，就是你提问要方便别人回答，这样才会有人帮你。
话说问题挺多的。
gs里面虽然是求平均，但是最好也把空间维度设置一下，这样比较标准。
ctl问题很大，你的水平空间既然求了平均，就是一个值，那就不要再写144*73了，z轴最好用level。你可以检查一下缺省值，如果你的数据有缺省值的话，缺省值应该也会变的。还有你真的贴完整了么？？？变量的描述哪里去了？？？

东篱

lqouc 发表于 2013-5-13 07:33

                               
登录/注册后可看大图

不时间不减分的问题，我的积分绝对够用，就是你提问要方便别人回答，这样才会有人帮你。
话说问题挺多的 ...

谢谢你的帮助，水平空间求平均，如果不写144*73该怎么改一下？谢谢！

新的ctl：

dset f:\data\hgt_index.grd
title 1948--3013-PMI
undef -2.56e33
xdef 144 linear 0 2.5
ydef 73 linear -90 2.5
zdef 1 levels 600 1
tdef 780 linear 00Z01jan1948 1mo

vars 6
a1 0 99 1th
a2 0 99 2th
a3 0 99 3th
a4 0 99 4th
a5 0 99 5th
temp 0 99 PMI
endvars

lqouc

因为你每一个时刻只有一个值，根本谈不上空间，所以就吧144和73 都改成1啊。数据是什么样子ctl就是什么样子。
xdef 1 linear 0 2.5
ydef 1 linear -90 2.5
话说新的ctl也不对，变量那里你只输出了temp，把其他的去掉。

东篱

lqouc 发表于 2013-5-14 11:42

                               
登录/注册后可看大图

因为你每一个时刻只有一个值，根本谈不上空间，所以就吧144和73 都改成1啊。数据是什么样子ctl就是什么样子 ...

懂了，实在是太感谢了

喵~

zdef  后面的1没用。。。。。。。。。。请无视我。

lqouc

喵~ 发表于 2013-5-15 16:16

                               
登录/注册后可看大图

zdef  后面的1没用。。。。。。。。。。请无视我。

这个不能无视，你说的很对，规范操作时很重要的

kakgqum80

话说问题挺多的。

                               
登录/注册后可看大图