一页四图之不同高度天气图，叠加高原地形、台风和季风低压中心

听海落雪

        最近因为要写个例总结，用 NCEP 资料画了几张天气图，跟大家分享一下。
（1）使用了 page.gs 的一页四图功能，分别绘制了 200 hPa、500 hPa、700 hPa 以及 地面 天气图。
（2）使用  tibetan.gs 绘制了东亚 700 hPa 以上的地形（具体文件和用法见附件）。
（3）台风路径数据 来自佛山天气網http://www.fs121.com/typhoon.aspx.
（4）大家要注意 tibetan命令的位置和用法。在 200 hPa 和 500 hPa 图上，等高线等浮于地形之上；700 hPa 和 地面图上，等高线等置于底层。
（5）本地市州边界地图使用 MeteoInfo C++ 版（下载地址：http://meteothink.org）制作，在此对王老师表示感谢。

'reinit'

'sdfopen F:\Data\GrADS\2019\slp.2019.nc'
*'set dfile 1'
'sdfopen F:\Data\GrADS\2019\hgt.2019.nc'
*'set dfile 2'
'sdfopen F:\Data\GrADS\2019\air.2019.nc'
*'set dfile 3'
'sdfopen F:\Data\GrADS\2019\uwnd.2019.nc'
*'set dfile 4'
'sdfopen F:\Data\GrADS\2019\vwnd.2019.nc'
*'set dfile 5'

'set grads off'
'set grid off'

'set lon 70 120'
'set lat 15 50'
'set t 735'

'set rgb 27    0   0 150'

********* 一页四图 *********

****** 左上图（200hPa） ******
'run page q1'

*'set parea 1.2 10.2   0.2 8.2'
'set xlopts 1 5 0.28'
'set ylopts 1 5 0.28'
'set xlab on'
'set ylab on'
'set xlabs 70`3.|80`3.|90`3.|100`3.|110`3.|120`3.`1E'
'set ylabs 15`3.|20`3.|25`3.|30`3.|35`3.|40`3.|45`3.|50`3.`1N'

'set lev 200'
*********** 高原地形 ***********
'set gxout contour'
'set csmooth on'
'set cint 4'
'set cstyle 1'
'set ccolor 27'
'set clopts 27  0.5  0.12'
'set clab forced'
'set clab masked'
'd hgt.2/10'
'tibetan 700'
*********************************
*********** 等高线 ***********
'set cint 4'
'set cstyle 1'
'set ccolor 27'
'set clopts 27  0.5  0.12'
'set clab forced'
'set clab masked'
'd hgt.2/10'
******************************
*********** 等温线 ***********
'set cint 4'
'set cstyle 2'
'set ccolor 2'
'set clopts 2  0.5  0.12'
'set clab forced'
'set clab masked'
'd air.3-273.15'
******************************
*********** 画风羽 ***********
'set gxout barb'
'set ccolor 1'
'set digsize 0.06'
'd uwnd.4*2.5;vwnd.5*2.5'
******************************
'set mpdset cnworld'
'set map 15 1 1'
'draw map'
'set mpdset pxdq_dj'
'set map 13 1 2'
'draw map'
******************************
'query w2xy 71.2 48'
x100=subwrd(result,3)
y100=subwrd(result,6)
'set string 1 c 6 0'
'set strsiz 0.45'
'draw string 'x100' 'y100'  `4  (a)'

'query w2xy 111.3 16.2'
x200=subwrd(result,3)
y200=subwrd(result,6)
'set string 1 c 6 0'
'set strsiz 0.30'
'draw string 'x200' 'y200'  `2  200 hPa'
******************************


****** 右上图（500hPa） ******
'run page q2'

*'set parea 1.2 10.2   0.2 8.2'
'set xlopts 1 5 0.28'
'set ylopts 1 5 0.28'
'set xlab on'
'set ylab on'
'set xlabs 70`3.|80`3.|90`3.|100`3.|110`3.|120`3.`1E'
'set ylabs 15`3.|20`3.|25`3.|30`3.|35`3.|40`3.|45`3.|50`3.`1N'

'set lev 500'
*********** 高原地形 ***********
'set gxout contour'
'set csmooth on'
'set cint 4'
'set cstyle 1'
'set ccolor 27'
'set clopts 27  0.5  0.12'
'set clab forced'
'set clab masked'
'd hgt.2/10'
'tibetan 700'
*********************************
*********** 等高线 ***********
'set cint 4'
'set cstyle 1'
'set ccolor 27'
'set clopts 27  0.5  0.12'
'set clab forced'
'set clab masked'
'd hgt.2/10'
******************************
*********** 等温线 ***********
'set cint 4'
'set cstyle 2'
'set ccolor 2'
'set clopts 2  0.5  0.12'
'set clab forced'
'set clab masked'
'd air.3-273.15'
******************************
*********** 画风羽 ***********
'set gxout barb'
'set ccolor 1'
'set digsize 0.06'
'd uwnd.4*2.5;vwnd.5*2.5'
******************************
'set mpdset cnworld'
'set map 15 1 1'
'draw map'
'set mpdset pxdq_dj'
'set map 13 1 2'
'draw map'
******************************
'query w2xy 71.2 48'
x100=subwrd(result,3)
y100=subwrd(result,6)
'set string 1 c 6 0'
'set strsiz 0.45'
'draw string 'x100' 'y100'  `4  (b)'

'query w2xy 111.3 16.2'
x200=subwrd(result,3)
y200=subwrd(result,6)
'set string 1 c 6 0'
'set strsiz 0.30'
'draw string 'x200' 'y200'  `2  500 hPa'
******************************
*********** 画台风和季风低压 ***********
'query w2xy 107.5 19.3'
x300=subwrd(result,3)
y300=subwrd(result,6)
'draw wxsym 41 'x300' 'y300' 0.8 2 5'

'query w2xy 81.0 24.5'
x200=subwrd(result,3)
y200=subwrd(result,6)
'set string 2 c 7 0'
'set strsiz 0.4'
'draw string 'x200' 'y200'  `4  D'
******************************


****** 左下图（700hPa） ******
'run page q3'

*'set parea 1.2 10.2   0.2 8.2'
'set xlopts 1 5 0.28'
'set ylopts 1 5 0.28'
'set xlab on'
'set ylab on'
'set xlabs 70`3.|80`3.|90`3.|100`3.|110`3.|120`3.`1E'
'set ylabs 15`3.|20`3.|25`3.|30`3.|35`3.|40`3.|45`3.|50`3.`1N'

'set lev 700'
*********** 等高线 ***********
'set cint 4'
'set cstyle 1'
'set ccolor 27'
'set clopts 27  0.5  0.12'
'set clab forced'
'set clab masked'
'd hgt.2/10'
******************************
*********** 等温线 ***********
'set cint 4'
'set cstyle 2'
'set ccolor 2'
'set clopts 2  0.5  0.12'
'set clab forced'
'set clab masked'
'd air.3-273.15'
******************************
*********** 画风羽 ***********
'set gxout barb'
'set ccolor 1'
'set digsize 0.06'
'd uwnd.4*2.5;vwnd.5*2.5'
******************************
*********** 高原地形 ***********
'set gxout contour'
'set csmooth on'
'set cint 4'
'set cstyle 1'
'set ccolor 27'
'set clopts 27  0.5  0.12'
'set clab forced'
'set clab masked'
'd hgt.2/10'
'tibetan 700'
*********************************
'set mpdset cnworld'
'set map 15 1 1'
'draw map'
'set mpdset pxdq_dj'
'set map 13 1 2'
'draw map'
******************************
'query w2xy 71.2 48'
x100=subwrd(result,3)
y100=subwrd(result,6)
'set string 1 c 6 0'
'set strsiz 0.45'
'draw string 'x100' 'y100'  `4  (c)'

'query w2xy 111.3 16.2'
x200=subwrd(result,3)
y200=subwrd(result,6)
'set string 1 c 6 0'
'set strsiz 0.30'
'draw string 'x200' 'y200'  `2  700 hPa'
******************************
*********** 画台风和季风低压 ***********
'query w2xy 107.5 19.3'
x300=subwrd(result,3)
y300=subwrd(result,6)
'draw wxsym 41 'x300' 'y300' 0.8 2 5'

'query w2xy 79.6 23.9'
x200=subwrd(result,3)
y200=subwrd(result,6)
'set string 2 c 7 0'
'set strsiz 0.4'
'draw string 'x200' 'y200'  `4  D'
******************************


****** 右下图（地面） ******
'run page q4'

*'set parea 1.2 10.2   0.2 8.2'
'set xlopts 1 5 0.28'
'set ylopts 1 5 0.28'
'set xlab on'
'set ylab on'
'set xlabs 70`3.|80`3.|90`3.|100`3.|110`3.|120`3.`1E'
'set ylabs 15`3.|20`3.|25`3.|30`3.|35`3.|40`3.|45`3.|50`3.`1N'

'set z 1'
*********** 等高线 ***********
'set cint 4'
'set cstyle 1'
'set ccolor 27'
'set clopts 27  0.5  0.12'
'set clab forced'
'set clab masked'
'd slp.1/100'
******************************
*********** 高原地形 ***********
'tibetan 700'
*********************************
'set mpdset cnworld'
'set map 15 1 1'
'draw map'
'set mpdset pxdq_dj'
'set map 13 1 2'
'draw map'
******************************
'query w2xy 71.2 48'
x100=subwrd(result,3)
y100=subwrd(result,6)
'set string 1 c 6 0'
'set strsiz 0.45'
'draw string 'x100' 'y100'  `4  (d)'

'query w2xy 112.2 16.2'
x200=subwrd(result,3)
y200=subwrd(result,6)
'set string 1 c 6 0'
'set strsiz 0.30'
'draw string 'x200' 'y200'  `2  surface'
******************************
*********** 画台风和季风低压 ***********
'query w2xy 107.5 19.3'
x300=subwrd(result,3)
y300=subwrd(result,6)
'draw wxsym 41 'x300' 'y300' 0.8 2 5'

'query w2xy 79.6 23.9'
x200=subwrd(result,3)
y200=subwrd(result,6)
'set string 2 c 7 0'
'set strsiz 0.4'
'draw string 'x200' 'y200'  `4  D'
******************************

'printim F:\Data\GrADS\2019\two_MD\pic\high_2019070312.png white'

*'enable print F:\Data\GrADS\2019\two_MD\pic\high_2019070312.gmf'
*'print'
*'disable print'

四叶草

谢谢分享！

fanyuxiang

写论文正好用得上，感谢楼主分享

lidd1130

river

真好看，支持一下！

听海落雪

四叶草 发表于 2019-8-18 08:35
谢谢分享！

听海落雪

fanyuxiang 发表于 2019-8-18 09:37
写论文正好用得上，感谢楼主分享

不用谢，如果有问题，随时交流。

听海落雪

river 发表于 2019-8-18 14:38
真好看，支持一下！

谢谢支持！！！！！！！！

river

听海落雪 发表于 2019-8-18 23:23
谢谢支持！！！！！！！！

建议可以升级grads 2.1或者更高的2.2版本了，出图质量更高了

听海落雪

river 发表于 2019-8-19 00:34
建议可以升级grads 2.1或者更高的2.2版本了，出图质量更高了

好的，马上去下载