GrADS副高脊线

帅蛋fly

'reinit'
'sdfopen E:\shuai\moonson\500hPa\hgt.mon.mean.nc'
'sdfopen E:\shuai\moonson\500hPa\uwnd.mon.mean.nc'
'sdfopen E:\shuai\moonson\500hPa\vwnd.mon.mean.nc'
'sdfopen E:\shuai\moonson\500hPa\air.mon.mean.nc'
'enable print  E:\shuai\moonson\500hPa\ave-2016-1.gmf white'
'set grid off'
'set grads off'

'set map 1 1 8'
'set mpdset hires'
'set mpdset cnworld'
'set lat 0 50'
'set lon 70 180'
'set lev 500'

'set xlopts 1 5 0.15'
'set ylopts 1 5 0.15'
*'set xlint 30'
*'set ylint 5'
'set cint 5'

*******求某年的夏季平均
'define aver=ave(hgt.1,time=jun2016,time=aug2016)'
'define uu = ave(uwnd.2,time=jun2016,time=aug2016)'
'define vv = ave(vwnd.3,time=jun2016,time=aug2016)'
'define tt = ave(air.4,time=jun2016,time=aug2016)+274.15'

'set gxout contour'
'set cint 4'
'set cthick 8'
'set ccolor 1'
'set cterp on'
'set csmooth on'
'd aver/10'

******画588线并突出
'set gxout contour'
'set clevs 588'
'set  cthick 11'
'set clab forced'
'd aver/10'

'set gxout contour'
'set clevs 0'
'set cthick 11'
'set ccolor 8'
'set clab off'
'd maskout(maskout(maskout(maskout(uu,-hcurl(uu,vv)),tt-250),lon-103),lat-10)'
*500hPa副高脊线（在负涡度区和温度大于250K的热带地区，u=0的线）

'set string 1 c 5'
'set strsiz 0.18'
'draw string 10 1.3 (l)'
'set line 1 1 5'
'draw rec 9.8 1 10.2 1.5'

'set strsiz 0.18 0.18
'
'draw title 2016 500hPa Summer Average Height Field'

'print'
'disable print'
结果如图：

river

不错哦，好看，也学到一个新方法

xuedeqiang

发一个500hpa+850hpa图

海豚

感谢群主分享，学习了