求助 画等熵位涡的 剖面图

heyufei

本人画出了等熵位涡在设定等压面后的剖面图 现在想画等熵位涡南北向的剖面图 修改后出来的全是一条一条直线 求教啊
*等熵位涡分析
'reinit'
'open f:/zhang/fnl_20130710_18_00.00.ctl'
'set lon 80 130'
'set lat 15 45'
'set lev 1000 100'
"define t=TMPprs"
"define u=ugrdprs"
"define v=vgrdprs"
"define pp=lev"
"define coriol=2*7.29e-5*sin(lat*3.1415/180)"
"define dudy=cdiff(u,y)/(111177*cdiff(lat,y))"
"define dvdx=cdiff(v,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))"
"define dt=t(z-1)*pow(1000/PP(z-1),0.286)-t(z+1)*pow(1000/PP(z+1),0.286)"
"define dp=100*(PP(z-1)-PP(z+1))"
"define dtdp=dt/dp"
"define part1="isen(dvdx,t,PP,320)
"define part2="isen(dudy,t,PP,320)
"define part3="isen(dtdp,t,PP,320)
"define pv320=-9.8*(coriol+part1-part2)*part3"
"set z 0"
'set lon 80 130'
'set lat 15 40'
"d pv320*1000000"

'printim f:/zhang/1000ipv1018.jpg x1024 y768 white'

按这个可以出来图 可是修改SET Z 和 SET LON 之后 就变成了一条一条直线

river

没懂你的意思，怎么修改了？

大气太阳茄

请问楼主问题解决没？
我也遇到了同样问题

小预报员

问题解决了吗？

小预报员

很重要的问题啊，怎么没有人回答啊！

fzghui

没人来解答么。。。

夏朗的芒果

因为这个程序本来算出来的位涡就只有一层，画出的垂直剖面图自然各个层次都是一致的。

heyufei

回复 谢谢各位关心 我的问题还是没有解决

heyufei

成功地画时间高度图真的是好难好难啊

want_data

位涡和等熵位涡有什么区别啊