【求助】grads中关于多个sum函数的应用

Iree_Aurora

请走过路过的大神们帮帮忙！grads小白对于这个问题想不明白了……
编写求位势高度场的多年平均场的gs文件时，因为所涉及的时间是没有规律的，所以就先sum，再平均，gs如下：
'reinit'
'xdfopen f:\ivf\annul_temp\hgt1961-2018\hgt.nc.1961-2018.ctl' *将1961-2018年的位势高度场合成一个ctl文件，便于下面根据不同时间求和

'set lev 850'
'set mpdset cnworld'
'define s=sum(hgt,t=700,t=708)+sum(hgt,t=731,t=737)+sum(hgt,t=2188,t=2212)+sum(hgt,t=5100,t=5106)+sum(hgt,t=5457,t=5471)+sum(hgt,t=5838,t=5845)+sum(hgt,t=5873,t=5879)+sum(hgt,t=6214,t=6221)+sum(hgt,t=6225,t=6233)+sum(hgt,t=8753,t=8767)+sum(hgt,t=15338,t=15344)+sum(hgt,t=16441,t=16447)+sum(hgt,t=17192,t=17202)+sum(hgt,t=18266,t=18275)+sum(hgt,t=18648,t=18654)+sum(hgt,t=18991,t=19004)'
*是的，这里就很长，因为需要求这么多段时间
'd s/165'
'gxprint f:\ivf\annul_temp\hgt1961-2018\hgt850.png png white'
;

然后运行，是可以画出图来的：

但是！！问题来了！！！
我用同样的gs文件（ctl提前运行好的），
ctl:
dset F:/IVF/annul_temp/air1961-2018/air.%y4.nc
title mean Daily Air Reanalysis
options template
options yrev
undef -9.9621e+36
xdef lon 144 linear 0 2.5  
ydef lat 73 linear -90 2.5
zdef level 17 levels 1000 925 850 700 600 500 400 300 250 200 150 100 70 50 30 20 10
tdef time 21184 linear 00Z01Jan1961 1440mn
vars 1
air=>air 17  t,z,y,x  mean Daily Air Temperture
endvars
gs:
'reinit'
'xdfopen f:\ivf\annul_temp\air1961-2018\air.nc.1961-2018.ctl'
'set lev 1000'
'set mpdset cnworld'
'define s=sum(air,t=700,t=708)+sum(air,t=731,t=737)+sum(air,t=2188,t=2212)+sum(air,t=5100,t=5106)+sum(air,t=5457,t=5471)+sum(air,t=5838,t=5845)+sum(air,t=5873,t=5879)+sum(air,t=6214,t=6221)+sum(air,t=6225,t=6233)+sum(air,t=8753,t=8767)+sum(air,t=15338,t=15344)+sum(air,t=16441,t=16447)+sum(air,t=17192,t=17202)+sum(air,t=18266,t=18275)+sum(air,t=18648,t=18654)+sum(air,t=18991,t=19004)'
'd s/165'
'gxprint f:\ivf\annul_temp\air1961-2018\air1000.png png white'
;

画多年平均温度场，就是这样:(

PS:试过把sum那一堆减少一些（减到了sum(air,t=8753,t=8767)，可以正常作图）

Question:

1、为什么位势高度场可以计算这么多sum，而温度场却不可以呢？
2、如果温度场必须要计算这么多时段的均温，请问该如何解决这个问题或者有没有更好的计算方案呢？

麻烦各位大佬帮忙解决一下这个疑问，谢谢大噶！！！

river

应该不存在这样的问题才对，位势高度可以，温度也必然是可以的······

Iree_Aurora

river 发表于 2020-9-12 20:14
应该不存在这样的问题才对，位势高度可以，温度也必然是可以的······

对啊，我也这么觉得，可是它输出来就是这样，请问您还有求平均温度场其他的思路吗（上面ctl也是今天看到您的帖子现学现用的，在此对您表示一下感谢哈~）

river

Iree_Aurora 发表于 2020-9-12 20:55
对啊，我也这么觉得，可是它输出来就是这样，请问您还有求平均温度场其他的思路吗（上面ctl也是今 ...

有没有什么报错信息吗？

Iree_Aurora

river 发表于 2020-9-12 20:58
有没有什么报错信息吗？

您看，没有报错

river

Iree_Aurora 发表于 2020-9-12 21:07
您看，没有报错

你这样的，你先试试把ctl里缺测值那一项用*注释掉。
然后还不行就直接用打开ctl的方式在GrADS里面输入命令，随机设置时间，看看原始资料都能不能出图，有的时候原始资料可能有问题

Iree_Aurora

river 发表于 2020-9-15 22:33
你这样的，你先试试把ctl里缺测值那一项用*注释掉。
然后还不行就直接用打开ctl的方式在GrADS里面输入命 ...

您看，我在命令窗口直接输入，8000之前的数据好像没有问题，但是加到后面就出不了图了，您看这是不是数据的问题，后面他有个Contouring: 1e+07 to 1.2e+08 interval 1e+07 的提示，您看是不是这个问题

Iree_Aurora

river 发表于 2020-9-15 22:33
你这样的，你先试试把ctl里缺测值那一项用*注释掉。
然后还不行就直接用打开ctl的方式在GrADS里面输入命 ...

我好像找到原因了！！谢谢您提供的思路！我想起来我1980年下的是ncep1的数据，1980年后是ncep2的数据，ncep2的数据似乎不能直接与ncep1用一个ctl打包起来，我可以再请教一下您吗，就是如果不重新下ncep11980年之后的数据，请问我该如何实现，这种1961-2018年位势高度场求平均呢（时间没有规律)

river

Iree_Aurora 发表于 2020-9-19 10:40
我好像找到原因了！！谢谢您提供的思路！我想起来我1980年下的是ncep1的数据，1980年后是ncep2的数据，nc ...

这种情况的话，我觉得思路就是分成两段分别批量描述分别求，前面一段加起来除以165，后面一段加起来除以165，最后再把这两个数字加起来

Iree_Aurora

river 发表于 2020-9-19 12:09
这种情况的话，我觉得思路就是分成两段分别批量描述分别求，前面一段加起来除以165，后面一段加起来除以1 ...

好的！谢谢您！我再去试试！