垂直螺旋度

松子茶

下面是我编辑的gs，之前画平面图可以出图，稍作修改画沿纬向的剖面就一直出现问题，求大神指点
*垂直螺旋度沿纬度剖面
'reinit'
'sdfopen f:\nc\421.nc'
'enable print f:\nc\Hp_N.gmf'
'set lon 110 120'
'set lat 35
'set lev 1000 100'
*'set mpdset cnworld'
'set grid off'
'set display color white'
'c'
'set grads off'
'set ylevs 1000 925 850 700 600 500 400 300 200 100'
'set xlopts 1 5 0.2'
'set ylopts 1 5 0.2'
'set csmooth on'
day=20
while(day<=22)
hour=0
while(hour<=18)
tt=(day-20)*4+hour/6+1
'set t 'tt''
'set xlint 1'
'set ylint 1'
* 绘制垂直螺旋度
'set gxout contour'
'define wd=hcurl(u,v)'
'define lx=-(wd*w*(-1000))/12.64'
'set ccolor 1'
'd lx*1e+6'
'draw title Hp for 2018Apr'day' 'math_mod(hour,24)''
'print'
'c'
hour=hour+6
endwhile
day=day+1
endwhile
'disable print'
'reinit'
；

                               
登录/注册后可看大图

松子茶

                               
登录/注册后可看大图

subtropical

楼主你好 请问你的相对螺旋度公式是什么 我最近要画风暴相对螺旋度，我搜了公式 h=V*(dU/dZ)-U*(dV/dZ) 请问你的公式也是这个吗？这个要怎么计算呢？