时间高度剖面出错

酒国 · 发表于 2015-12-20 16:59:01

登录后查看更多精彩内容~

您需要登录才可以下载或查看，没有帐号？立即注册

x

本帖最后由酒国于 2015-12-20 17:01 编辑

大家好，我想做一物理量的时间剖面图，可是画出来的有问题。如果画一个时次的经向和纬向高度剖面是可以的，请问哪里出错了。

'reinit'
'open d:\test\20110923\fnl_20110923_00_00.ctl'
'set lon 102 125'
'set lat 12 27'
'set t 23 26'
'set lev 1000 100'
'set grads off'
'set grid off'
'define h1=(1+0.86*qs)*1004*TMPprs'
'define dtx=cdiff(h1,x)'
'define dty=cdiff(h1,y)'
'define dx=cdiff(lon,x)*3.1416/180'
'define dy=cdiff(lat,y)*3.1416/180'
'set gxout contour'
'set cthick 6'
'set ccolor 1'
'set xlopts 1 6 0.16'
'set ylopts 1 6 0.16'
'set clopts 1 6 0.10'
'set csmooth on'
'set lat 19.5'
'set t 23 26'
'set lev 1000 100'
'define aa=-1*(VGRDprs*dty/dy)/6.37e6*100000'
'define aaa=ave(aa,lon=102,lon=124)'
'd aaa'

四叶草 · 发表于 2015-12-20 19:27:34

本帖最后由四叶草于 2015-12-20 19:28 编辑

时间剖面图应该是固定经纬度，高度和时间变化，而你的经度没有固定，请楼主仔细看看！

river · 发表于 2015-12-20 19:38:05

四叶草发表于 2015-12-20 19:27
时间剖面图应该是固定经纬度，高度和时间变化，而你的经度没有固定，请楼主仔细看看！

你看出来问题了，我就不说了，就占个楼层吧

酒国 · 发表于 2015-12-20 19:39:53

四叶草发表于 2015-12-20 19:27
时间剖面图应该是固定经纬度，高度和时间变化，而你的经度没有固定，请楼主仔细看看！

谢谢！已经可以出图了

ypliu2015 · 发表于 2015-12-20 20:19:18

略看不懂

{:5_196:}

付亚男 · 发表于 2015-12-20 22:40:39

兴安盟的刘大神么……

酒国 · 发表于 2015-12-21 10:31:03

river 发表于 2015-12-20 19:38
你看出来问题了，我就不说了，就占个楼层吧

我还有点不明白，希望river大神能够指点一下。我看的论文上面他做的图是沿32.5N（110.5-125.5E)做的时间高度剖面，如下图

他这个在编写程序时经度是变化的吗？我现在的程序需要做什么调整？

酒国 · 发表于 2015-12-21 10:42:33

酒国发表于 2015-12-21 10:31
我还有点不明白，希望river大神能够指点一下。我看的论文上面他做的图是沿32.5N（110.5-125.5E)做的时 ...

'set lon 111.1'
'set t 21 31'
'set lev 1000 300'
'set cint 500'
'set xaxis -5 5 1'
'set clab forced'
'define aa=-1*((UGRDprs*dtx)/(cos(lat*3.1416/180)*dx)+VGRDprs*dty/dy)/6.37e6*100000'
'define aaa=ave(aa,lat=12,lat=27)'
'set lat 19.7'
这种思路对不对：我把纬度先求个平均，还有在给出一点纬度，这样做出来的图是不是就是纬向平流的时间高度剖面？

river · 发表于 2015-12-21 11:24:24

酒国发表于 2015-12-21 10:42
'set lon 111.1'
'set t 21 31'
'set lev 1000 300'

没看到原文，所以不太好理解···

酒国 · 发表于 2015-12-21 14:12:16

river 发表于 2015-12-21 11:24
没看到原文，所以不太好理解···

能量平流.pdf (740.88 KB, 下载次数: 7)

		自动登录	找回密码
密码			立即注册

[分享资料] 时间高度剖面出错

登录后查看更多精彩内容~