想做散度沿高度的垂直剖面图，现阶段只会写单层散度图，请问如何加入高度循环

蘑菇酱

查了气象家园很多帖子，只找到了关于时间的垂直剖面，我的gs脚本如下
'reinit'
'sdfopen G:\2017\uwnd2017.nc'
'sdfopen G:\2017\vwnd2017.nc'
'sdfopen G:\2017\rhum2017.nc'
'sdfopen G:\2017\shum2017.nc'
'sdfopen G:\2017\omega2017.nc'
'sdfopen G:\2017\air2017.nc'
'set t 735 748'
'set lev  850'
'define qdiv=hdivg(uwnd.1*shum.4*1e3/9.8,vwnd.2*shum.4*1e3/9.8)*(1e5)'
'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'

'set gxout shaded'
'set cint 0.3'

'd qdiv'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.3'
'set cthick 11'

'd qdiv'

蘑菇酱

好像帖脚本的时候把经度纬度设置删掉了

蘑菇酱

G:\2017\tu\水汽通量散度\7月5日06时\500.png

蘑菇酱

这是单层图

蘑菇酱

重新查找了以后明白了，就是在开始定义的时候加上‘set lev 1000 300’，然后出图前固定经度/纬度。以后还是把问题想清楚以后多看看，贴一下改了以后出图的gs
'reinit'
'sdfopen G:\2017\uwnd2017.nc'
'sdfopen G:\2017\vwnd2017.nc'
'sdfopen G:\2017\rhum2017.nc'
'sdfopen G:\2017\shum2017.nc'
'sdfopen G:\2017\omega2017.nc'
'sdfopen G:\2017\air2017.nc'
'set t 735'

'set lev 1000 300'
'define qdiv=hdivg(uwnd.1*shum.4*1e3/9.8,vwnd.2*shum.4*1e3/9.8)*(1e5)'

'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'

'enable print G:\2017\1\qdiv.gmf'
'set lon 104 126'
'set lat 31'
'set zlog on'

'set ylevs 1000 925 850 700 600 500 400 300'
'set gxout shaded'
'set cint 0.3'

'd qdiv'
'run cbar.gs'
'set ylpos 0 r'
'set gxout contour'
'set cint 0.3'
'set cthick 11'

'd qdiv'

阳光暴雨

‘set z 1 8’
'set zlog on'

'set ylevs 1000 925 850 700 600 500 400 300'