辐散风和速度势的图

chuanhao

我通过论坛一个帖子画了一张辐散风和速度势的图，查书弄不懂速度势的物理意义后百度，上面说速度势为正，是辐散运动，为负的时候，是辐合运动。但是我画出来的图辐合中心是正值。求各位大神解答，不甚感激！

centaurus

我也出现这种问题，直接用的NCEP的再分析的逐月CHI画的。想问楼主解决了吗？

zhouling

速度势为正，是辐散运动，为负的时候，是辐合运动。可能是你编的程序出错了

听雨轩

我觉得速度为正的时候是幅合 负值表示辐散   具体参加文献Influence of El Niño Modoki on spring rainfall over south China （Juan Feng1 and Jianping Li1）图2a

chuanhao

听雨轩 发表于 2014-3-21 22:11

                               
登录/注册后可看大图

我觉得速度为正的时候是幅合 负值表示辐散   具体参加文献Influence of El Niño Modoki on spring rai ...

十分感谢，

chuanhao

zhouling 发表于 2014-3-21 23:21

                               
登录/注册后可看大图

速度势为正，是辐散运动，为负的时候，是辐合运动。可能是你编的程序出错了

恩，应该是程序编写的符号问题

nuistajun

有木有速度势的gs？？

chuanhao

nuistajun 发表于 2014-3-28 21:44
有木有速度势的gs？？

'reinit'
'open F:\xin\sqtl\qu.ctl'
'open F:\xin\sqtl\qv.ctl'
'open F:\xin\sqtl\avequ.ctl'
'open F:\xin\sqtl\aveqv.ctl'
'define UU=qu.1-qu.3'
'define VV=qv.2-qv.4'
'set ylint 10'
'SET XLINT 20'
'set xlopts 1 5 0.16';'set ylopts 1 5 0.16'
*'set mpdset cnworld'
'set map 1 1 1'
'set grads off'
'set grid off'
***************Computing Velocity Potential from Wind Components**********
'define chi=fish_chi(uu,vv)'   
'chi=chi-aave(chi,global)'                                                                                                            
'set gxout contour'
'set lat 0 50'
'set lon 50 160'
*'set cthick 5'
'd chi/1e6'  
****Calculating the Divergent Wind****                                    
'define one=1+0*lat'
'uchi = - muadv(one,chi)'
'vchi = - mvadv(one,chi)'
'set gxout vector'
*'set ccolor 15'
'set lat 0 50'
'set lon 50 160'
'set cmin 2'
'set cmax 20'
'set arrscl 0.5 40'
'd uchi;vchi'
'q time'
xx=subwrd(result,3)
'set line 1 1 5'
'draw rec 4.15 4.18  5.92 5.43'
'draw title Divergence wind and velocity potential of EXP'
'enable print F:\xin\dsdb\shiyan\shuiqifusan1.gmf'
'print'
'disable print'

;

nuistajun

chuanhao 发表于 2014-3-30 12:09
'reinit'
'open F:\xin\sqtl\qu.ctl'
'open F:\xin\sqtl\qv.ctl'

可不可以画成x轴是经度，Y轴是高度这样的速度势。

'reinit'
'open C:\Users\Lenovo\Desktop\wind\positive\u1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\v1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\w1.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\u3.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\v3.ctl'
'open C:\Users\Lenovo\Desktop\wind\positive\w3.ctl'
'set x 1 144'
'set y 1 73 '
'set z 1 12'
'set annot 8 8'
'set mpdraw on'
'set map auto'
'set mpdset cnworld'
'set grads off'
'set grid off'
'define chi=fish_chi(uwnd.1-uwnd.4,vwnd.2-vwnd.5)'
*pull dummy
'chi=chi-aave(chi,global)'
'set gxout contour'
'set cthick 9'
'set lon 50 300'
'set lat 0'
'set z 1 12'
'set zlog on'
'd ave(chi/1e6,lat=-10,lat=20)'
*pull dummy
'set gxout vector'
'set lon 50 300'
'set lat 0'
'set z 1 12'
'set zlog on'
'd ave(uwnd.1-uwnd.4,lat=-10,lat=20); ave(omega.3-omega.6,lat=-10,lat=20)'
*pull dummy
'printim C:\Users\Lenovo\Desktop\wind\positive\grid('a').png white'
;

chuanhao

nuistajun 发表于 2014-3-30 14:34
可不可以画成x轴是经度，Y轴是高度这样的速度势。

'reinit'

应该可以吧，只是换了一维而已