grads如何画大气静力稳定度的垂直剖面图

凉茶啦啦啦啦啦

想画一个横坐标为纬度，纵坐标为高度的大气静力温度度图，计算的时候，每一层的大气静力稳定度值为上下两层的假相当位温差除以高度差，但只能画出一层的平面图，不知道怎么画随纬度和高度变化的剖面图，用下面的gs文件只能画出一条线，完全不对，着急！！
'reinit'
'open K:\271_h\data0.ctl'
'set z 1 20'
'set t 1'
'set lon 190 220'
'set lat 70 80'
'define these=theta'
iz=1
while(iz<=20)
'set z 'iz
'these=theta(z=2)-theta(z=1)'
'define dh=height(z=2)-height(z=1)'
'these=these/dh'
iz=iz+1
endwhile
'set z 1 20'
'set lon 202'
'set gxout shaded'
'd these'