利用grads2.1制作散度、垂直速度、水汽通量散度剖面图的gs

qiuba

经过一段时间学习，自己调试绘制了一些有关水汽通量散度的gs，现给大家分享一下，不足之处，请多提宝贵意见。
以2018年4月20-23日华北地区大范围低云降水为背景，以首都国际机场为例绘制了三张剖面图，主要是想通过利用剖面图这个工具进一步分析水汽通量散度在日常天气预报中的应用问题。脚本gs如下：
'reinit'
'sdfopen d:\grads\data\uwnd.2018.nc'
'sdfopen d:\grads\data\vwnd.2018.nc'
'sdfopen d:\grads\data\rhum.2018.nc'
'sdfopen d:\grads\data\shum.2018.nc'
'sdfopen d:\grads\data\omega.2018.nc'
'sdfopen d:\grads\data\air.2018.nc'
'set t 441 452'
'set cachesf 1219392 bytes'
*-----水平散度-----
'set lev 1000 100'
'define div=hdivg(uwnd.1,vwnd.2)*1e5'
*-----水汽通量散度-----
'set lev 1000 300'
'define qdiv=hdivg(uwnd.1*shum.4*1e3/9.8,vwnd.2*shum.4*1e3/9.8)*1e5'
*--------GRADS设置-----------
'set grads off'
'set grid off'
'set csmooth on'
'set map 15 1 9'
'set parea 1 10.2 1.5 8'
'set xlopts 1 6 0.14'
'set ylopts 1 6 0.14'
'set font 0'
*---------输出图形1------------
'set lon 116.06'
'set lat 40.08'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300'
'set xlabs 08|12|16|20|00|04|08|12|16|20|00|04|08|12|16|20|00|04|08'
'set gxout shaded'
'set cint 0.3'
'set cmax 2'
'd qdiv'
'run d:\grads\cbarn.gs'
'set ylpos 0 l'
'set gxout contour'
'set cint 0.3'
'set cthick 5'
'set clskip 1'
'd qdiv'
'writeslhz 0.5 4.3 百帕 FSJ 3.0 1 1.0 0.5  90'
'writeslhz 0.6 0.9 4月20日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 2.65 0.9 4月21日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 5.7 0.9 4月22日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 8.75 0.9 4月23日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 5.2 0.04 时间 FSJ 3.0 1 1.0 0.5  0'
'writeslhz 0.6 8 水汽通量散度 FSJ 3.0 1 1.0 0.5  0'
'gxprint d:/grads/1/qdiv.png white'
'c'
*---------输出图形2------------
'set grads off'
'set lon 116.06'
'set lat 40.08'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
'set xlabs 08|12|16|20|00|04|08|12|16|20|00|04|08|12|16|20|00|04|08'
'set gxout shaded'
'set cint 0.1'
'set cmax 3'
'd div'
'run d:\grads\cbarn.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.2'
'set cthick 5'
'set clskip 1'
'd div'
'writeslhz 0.5 4.3 百帕 FSJ 3.0 1 1.0 0.5  90'
'writeslhz 0.6 0.9 4月20日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 2.65 0.9 4月21日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 5.7 0.9 4月22日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 8.75 0.9 4月23日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 5.2 0.04 时间 FSJ 3.0 1 1.0 0.5  0'
'writeslhz 0.6 8 散度 FSJ 3.0 1 1.0 0.5  0'
'gxprint d:/grads/1/div.png white'
'c'
*---------输出图形3------------
'set lev 1000 100'
'set grads off'
'set lon 116.06'
'set lat 40.08'
'set zlog on'
'set ylevs 1000 925 850 700 600 500 400 300 250 200 150 100'
'set xlabs 08|12|16|20|00|04|08|12|16|20|00|04|08|12|16|20|00|04|08'
'set gxout shaded'
'set cint 0.05'
'set cmax 5'
'd omega.5'
'run d:\grads\cbarn.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.05'
'set cthick 11'
'd omega.5'
'writeslhz 0.5 4.3 百帕 FSJ 3.0 1 1.0 0.5  90'
'writeslhz 0.6 0.9 4月20日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 2.65 0.9 4月21日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 5.7 0.9 4月22日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 8.75 0.9 4月23日 HTJ 2.0 1 0.8 0.5  0'
'writeslhz 5.2 0.04 时间 FSJ 3.0 1 1.0 0.5  0'
'writeslhz 0.6 8 垂直速度 FSJ 3.0 1 1.0 0.5  0'
'gxprint d:/grads/1/omega.png white'
;
效果图如图所示
从水汽通量散度图中可以看出，4月20日16时前机场上空水汽通量散度值为负值，表明水汽聚集，为低云降水提供有利水汽条件；从16时开始至22日零时为止，低层850百帕以下水汽通量散度为正值，表明水汽通量散度值流失明显，其原因正是由于低云降水造成的；从22日零时开始至16时止，主要是为负值的水汽聚集，此阶段主要是维持低云天气，降水不明显（仅08时前出现的一些分散性降水）；再之后由于500百帕以下水汽通量散度值为正值，大气中整层水汽流失严重，低云天气结束，云层逐渐抬高，天气也逐渐转好，与实际地面观测非常吻合。
另外，可能是由于版本不兼容问题，绘图中出现数字与汉字不在一条线上，影响美观，以后我再学习和改正。

看云淡风清

2.1画的效果比2.0不能好再多

qiuba

其实opengrads2.0画出的效果更好一些，主要是虚线更美观一些，但opengrads2.0不识别鲁局的writeslhz命令，
图注部分无法显示，目前尚未解决，有些遗憾

river

我的是数字和汉字是在同一水平线上的，只是和你的一样数字和汉字的粗细不一致。
建议去鲁局的帖子里反馈http://bbs.06climate.com/forum.p ... p;extra=&page=1，或者看脚本里鲁局的QQ，然后联系反馈

qiuba

谢谢river，确实很漂亮，是用opengrads画的吧，虚线更美观，数值标注也很规范，我也画过，只是它不识别鲁局的writeslhz命令，我也求教过鲁局，前几天鲁局说可能是我用的版本高，与writeslhz不兼容吧，也问了一些我画图的环境、系统之类的问题，也在解决这个问题。方便的话能附上您的脚本程序吗？我想对照一下，看能不能找到问题症结。

qiuba

river 发表于 2018-7-2 22:03
我的是数字和汉字是在同一水平线上的，只是和你的一样数字和汉字的粗细不一致。
建议去鲁局的帖子里反馈ht ...

这是我用opengrads2.0气象家园整合版本画的图，方便的话，能否帮我试一下，看问题在哪里，是不是我的脚本路径设置问题？

qiuba

这是opengrads的运行提示

qiuba

river 发表于 2018-7-2 22:03
我的是数字和汉字是在同一水平线上的，只是和你的一样数字和汉字的粗细不一致。
建议去鲁局的帖子里反馈ht ...

我看您的图更漂亮，再想问一下，是用哪个版本绘制的？

river

qiuba 发表于 2018-7-3 07:09
谢谢river，确实很漂亮，是用opengrads画的吧，虚线更美观，数值标注也很规范，我也画过，只是它不识别鲁局 ...

我就是直接复制你的脚本，资料用的是2012年的，我手头上就这个资料。opengrads的版本是2.1.a2.oga.1，windows版的。估计以后咱们都用不上grads最新windows版了，好久没有更新过了。
矢量汉字如果2.1可以用的话，那么2.0肯定可以用的，肯定还是你放的位置有问题。而且两个版本grads不能同时安装在一台电脑上吧，环境路径会有冲突······

四叶草

还是习惯用2.0，2.1用着很纠结！楼主图不错，谢谢分享。