求助，困扰好久但感觉简单的问题

天道酬勤

楼主首先对四年6月份的500hpa高度场做了个平均，这是提取每年的gs文件：'reinit''sdfopen e:\osci\situation\hgt.mon.mean.nc'
'set fwrite e:\osci\situation\4.grd'
'set gxout fwrite'
'set lev 500'
'set t 654'
'd hgt'
'disable fwrite'
;
这个算出平均的ctl文件：
dset e:\osci\situation\han6.grd
title paomian
undef -9.99e+08
XDEF 145 linear 0 2.5
YDEF 73 linear 0 2.5
ZDEF 1 levels 500
TDEF 1 LINEAR jun1979 1mon
VARS 1
h 0 99 paomian
endvars


然后，想做个距平，意思是用这个去减30年平均。这个是gs文件
'reinit'
'sdfopen e:\osci\situation\hgt.mon.mean.nc'
'open e:\osci\situation\hgt.ctl'
'set lev 500'
'define pingjun=ave(hgt,t=390,t=750,12)'
'd h.2-pingjun'
;

感觉这是个简单的问题，可为什么解决不了，老出错这是错误原因。   请各位指点一下。谢谢了！！！

做个霸气的木头

你把我绕晕了，哪个文件对应哪个，到底是谁的距平。

mofangbao

'define pingjun=ave(hgt,t=390,t=750,12)'
问题应该是出在这句话上，这里的隐循环，你最后只会得到最后一个t出来的变量，没法得到你想要的效果，你自己考虑下

冰枫

提取数据也不对吧，就提取了一个时次的数据呀，没有循环怎么提取4年的呢？

天道酬勤

做个霸气的木头 发表于 2012-9-18 22:20

                               
登录/注册后可看大图

你把我绕晕了，哪个文件对应哪个，到底是谁的距平。

第一个是从nc中提取一年的位势高度场，因为我提取的四年没有规律，所以我一年一年的提取，然后求平均。
第二个是求出平均之后的ctl文件。
第三个是求距平的gs文件。

天道酬勤

冰枫 发表于 2012-9-18 23:04

                               
登录/注册后可看大图

提取数据也不对吧，就提取了一个时次的数据呀，没有循环怎么提取4年的呢？

不是，第一个是我提取一年的gs。总共提取了4次，然后求的平均。

天道酬勤

mofangbao 发表于 2012-9-18 22:26

                               
登录/注册后可看大图

'define pingjun=ave(hgt,t=390,t=750,12)'
问题应该是出在这句话上，这里的隐循环，你最后只会得到最后一 ...

想了一晚上，没理清头绪，之前没学过grads，都是现学的，可能基础不行吧。麻烦清风，给我说一下如何改吧。谢谢了！

mofangbao

天道酬勤 发表于 2012-9-19 10:12

                               
登录/注册后可看大图

想了一晚上，没理清头绪，之前没学过grads，都是现学的，可能基础不行吧。麻烦清风，给我说一下如何改吧。 ...

你本身的问题提的就不清楚，你需要自己多看书，多看例子，类似的帖子在论坛是有的，你多找找

做个霸气的木头

天道酬勤 发表于 2012-9-19 10:09

                               
登录/注册后可看大图

第一个是从nc中提取一年的位势高度场，因为我提取的四年没有规律，所以我一年一年的提取，然后求平均。
...

其实这些都可以在一个gs文件里完成的，既然你分步做了，那就这样吧。你先把你的
最后的那个gs改成这样试试
'reinit'
'sdfopen e:\osci\situation\hgt.mon.mean.nc'
'open e:\osci\situation\hgt.ctl'
'set t 1'
'set lev 500'
'define pingjun=ave(hgt,t=390,t=750,12)'
'set t 1'
'd h.2(t=1)-pingjun'
;
如果还是不对，你画画你求完平均的那个hgt文件，看能不能出图。步骤多就逐步检查吧~

天道酬勤

做个霸气的木头 发表于 2012-9-19 12:09

                               
登录/注册后可看大图

其实这些都可以在一个gs文件里完成的，既然你分步做了，那就这样吧。你先把你的
最后的那个gs改成这样试 ...

好了，太感谢了。。。可我想知道，为什么要有两个’ set t 1‘ 哈？