求grads脚本，用FNL文件计算绝对角动量的

hymn514

想做出类似图片的那种剖面图，查了很多文献一头雾水不知道应该怎么画出来，求大神赐教，如果有现成的脚本那就更好啦。

_whp

角动量是什么变量，还是计算量？

_whp

'Reinit'
'Open d:\fnl.ctl'　＊是删除＊＊后内容的CTL文件
'set grid on'
'set grads off'
'set lon 70 150'
'set lat 10 60'
'set lev 1000 100'
'set display color white'
'c'
'set t 1'
'define u = UGRDprs'
'define v = VGRDprs'
'define tm= Tmpprs'
'define rh= Rhprs'
'set rgb 31   0 168 255'
'set rgb 32   0 103 255'
'set rgb 33  47 154   0'
'set rgb 34  38 255   0'
'set rgb 35 255 228   3'
'set rgb 36 255 203   3'
'set rgb 37 255 178   3'
'set rgb 38 255 154   0'
'set rgb 39 255  77   0'
'set rgb 40 255   0   0'
'set rgb 41 198   2   4'
'set rgb 42 157   0   0'
'set rgb 43 160   0 160'
'set rgb 44 213   0 217'
'set rgb 45 239   0 255'
'set rgb 46  20  20  20'
'define vi =mag(u,v)'
'set clevs 5.5 8.0 10.8 13.9 17.2 20.8 24.5 28.5 32.7 37.0 41.5 46.2 51.0 56.1 61.3'
'set ccols 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 46'
'set gxout Contour'
'set csmooth on'
'set cmin 5'
'set lev 1000 300'
'set lon 114.5'
'set lat 33 45'
'set zlog on'
'set csmooth on'
'd vi'

'define es=(6.112*exp(17.67*(tm-273.15)/(tm-29.65)))'
'define qs=rh*(0.62197*es/(lev-0.378*es))/100'
'define vi=tm*pow((1000/lev),(0.2854*(1.0-0.28*qs)))'
'set cthick 5'
'set cint 5'
'set gxout ccolor rainbow'
'set lev 1000 300'
'set lon 114.5'
'set lat 33 45'
'set zlog on'
'set csmooth on'
'd vi'

*图上标题
'draw title 2019-08-11 06Z Mag(u,v)(114.5`ao`nE)'

'printim D:\FNL\GIF\20190811_06_Mag.gif white'
'c'
'quit'

_whp

Open d:\fnl.ctl'　＊是删除＊＊后内容的CTL文件，CTL文件中有很多注释，请将每行的注释删除。另外图中的M是不是全风速，我的GS上的是全风速

_whp

                               
登录/注册后可看大图

_whp

是否这个结果？

hymn514

_whp 发表于 2019-11-7 11:57
是否这个结果？

不太像啊~~~~~~~~