等熵位涡的剖面

流浪的心

如何求等熵面位涡

jingshen

等熵面上位涡怎么画呢，求ncl脚本

依恋。

mofangbao 发表于 2012-7-19 09:05
剖面图就是固定的维度不一样而已，把其中的一个维度固定，原来高度层变为变化的，画图的语句还是一样的

老师您好，不太理解您说的把其中一个维度固定，为了画剖面，算涡度是没有办法把经度和维度固定的，算好用等压面转成等熵面，等压面也不能固定，我能想到的思路就是各个等熵面分别输出，于是ctl中zdef设置就不会了。麻烦您看一下，没有循环是可以出图的，设置pv为循环变量，相当于输出不同等熵面的物理量，这样的话，ctl中zdef就不知道怎么写了，因为纵坐标已经不是常规z轴了
'reinit'
'open D:\data\fnl\data\typhoon.ctl'
'set gxout fwrite'
'set fwrite D:\data\fnl\dengshang.grd'
pv=325
while(pv<327)
'set lon 90 140'
'set lat 10 60'
'set lev 1000 100'
'set t 2'
"define tt=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=tt(z-1)*pow(1000/PP(z-1),0.286)-tt(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,tt,PP,pv)
"define part2="isen(dudy,tt,PP,pv)
"define part3="isen(dtdp,tt,PP,pv)

"define pv320=-9.8*(coriol+part1-part2)*part3"
'set z 1'
"display pv320*1000000"
pv=pv+1
endwhile

'disable print'
function isen(field,tgrid,pgrid,tlev)
*----------------------------------------------------------------------
* Bob Hart (hart@ems.psu.edu) /  PSU Meteorology
* 2/26/1999

依恋。

mofangbao 发表于 2012-7-19 09:05
剖面图就是固定的维度不一样而已，把其中的一个维度固定，原来高度层变为变化的，画图的语句还是一样的

万分期待您的回复。

依恋。

mofangbao 发表于 2012-7-19 09:05
剖面图就是固定的维度不一样而已，把其中的一个维度固定，原来高度层变为变化的，画图的语句还是一样的

找到原因了，打扰您了，还是谢谢你的提点

zzgf

请问楼主画出垂直剖面图了吗