如何画温度距平剖面图？？

爱喝柠檬汁

我觉得可能是我的垂直层设置的不对，@river

爱喝柠檬汁

river 发表于 2017-12-17 09:02
怎么说呢，单就那个函数看，aave求的是区域平均。你把它放在grads里面如果设置了层次那就是那一层的那个 ...

大神帮我看看呀，
具体内容在上面几楼，

river

爱喝柠檬汁 发表于 2017-12-17 16:53
我用楼上的方法画出来的效果如图，我觉得似乎不是我想要的文献中的效果，gs如下：

'reinit'

文献里没有说是区域平均啊，你自己想一下区域平均算出来的是一个值，他说的应该是每个格点上的平均值，而且你也没搞清楚他说的是时间的平均还是纬度方向的平均······

爱喝柠檬汁

请问每个格点上的平均值怎么算呢？
还有我觉得是纬度上的平均，如下图，时间是固定某一时次的，我之前截图少截了一块

爱喝柠檬汁

我弄出来了，不知道对不对，请@river
您帮我看看呗，还有没有什么错误，嘻嘻
在楼下

爱喝柠檬汁

'reinit'
'sdfopen e:\aaa.nc'
'set gxout fwrite'
'set fwrite e:\tzjuping.dat'
'set lon 123'
'set lat 45 60'
'set t 109'


i=1000
while(i>=100)
'set lev 'i

'define avez=ave(z,lat=45,lat=60)'
'define avet=ave(t,lat=45,lat=60)'

'd z/9.8-avez/9.8'
'd t-avet'

i=i-100
endwhile

'disable fwrite'
'set gxout contour'
；
用这个gs输出一个dat文件，里面的值就是距平值
下面是对应的ctl
dset e:\tzjuping.dat
undef -9.99E+33
title juping
xdef   1 linear  123  0.5
ydef  31 linear   45  0.5
zdef  10 levels  1000 900 800 700 600 500 400 300 200 100
tdef   1 linear  00z28aug2017 6hr
vars  2
z     10  99 ** <z>     hpa
t     10  99 ** <t>     k
endvars

然后是画图的gs
'reinit'
'open E:/tzjuping.ctl'
'set grads off'
'set mpdset hires'
'set map 15 1 1'

'set t 1'
'set x 1'
'set lat 45 60'
'set lev 1000 100'
'set zlog on'
'set gxout shaded'
'd z'
'run cbarn.gs'

'set gxout contour'
'set ccolor 15'
'set cstyle 1'
'set cthick 2'
'd t'

'printim E:/11.png white'
;
然后是出图

爱喝柠檬汁

river 发表于 2017-12-17 22:23
文献里没有说是区域平均啊，你自己想一下区域平均算出来的是一个值，他说的应该是每个格点上的平均值，而 ...

您能帮我看看有什么错误吗？我弄得是不是文献中的意思呢？嘻嘻

爱喝柠檬汁

@river @river @river