【已解决】求日平均的gs问题

muggle

我从ECMWF网上下的2009年850 和 500 hPA环流场nc资料 如图所示


（不知道为什么这里q ctlinfo显示出来的z是850 和 -350 但是set z 2的时候确实是500 所以暂时不看这个小问题）


因为下载的是每6小时一次的资料，求日平均的gs我是这么写的



得到grd文件后，我用ctl出图检验的时候发现，在某个时间之前能正常出图，而后面的图就画不出来，好像是因为缺测？？还是我的gs写错了？各位大大帮我看一下，麻烦了

看云淡风清

你这个资料我没用过，不过看着有点奇怪，变量名和高度场名是一样的，还有你这个是两层循环，按照下面的格式试试看吧：
i=1
while(i<=1457)
'set t 'i''
'set lev 850'
'define ght=ave（）'
'd ght'
'set lev 500'
'define ght=ave（）'
'd ght'

lqouc

sdfopen放在循环的外面去

关于那个小问题：
”不知道为什么这里q ctlinfo显示出来的z是850 和 -350 但是set z 2的时候确实是500 所以暂时不看这个小问题”
人家写的zdef是linear，后面跟的是初值和步长，不是平时常用的levs。这种数据一般不会出错的。

muggle

lqouc 发表于 2017-2-27 18:23
sdfopen放在循环的外面去

关于那个小问题：

恩恩。level那个是我自己太粗心了，没注意到。谢谢你的提醒~
但是你说的吧sdfopen放在循环外面，运行得到的文件还是出现了同样的问题

这是我修改之后的gs 以及出错的截图 麻烦你帮我看看是什么原因

muggle

看云淡风清 发表于 2017-2-27 17:52
你这个资料我没用过，不过看着有点奇怪，变量名和高度场名是一样的，还有你这个是两层循环，按照下面的格式 ...

嗯，我明白你的意思。变量名重复我确实没考虑过这方面，但是你说的修改方法改的是两层循环，我觉得这部分我没出错

lqouc

循环写的不对
应为:
ii=1
while(ii<=2)
i=1
while(i<=1457)
之前的写法ii=2的时候，i不是从1开始的而是直接就是1461。

ps：楼主以后提问请直接粘贴代码，不要贴截图，否则别人帮你改还要重新输入代码

muggle

lqouc 发表于 2017-2-28 09:31
循环写的不对
应为:
ii=1

成功了！谢谢你的指导！

这种细节上的错误在自己的思维里面很难改正出来。。

下次我会改正的，不贴图了。
谢谢~

lqouc

muggle 发表于 2017-2-28 10:29
成功了！谢谢你的指导！

这种细节上的错误在自己的思维里面很难改正出来。。

熟练了就好了。
欢迎常来论坛~