请教水汽通量散度的时间剖面图怎么画

璧夏

我会画一层的水汽通量散度，又在论坛里学习其他同胞的时间剖面图然后写了下面的程序，
'reinit'
'sdfopen D:\lw\zl\uwnd.2015.nc'
'sdfopen D:\lw\zl\vwnd.2015.nc'
'sdfopen D:\lw\zl\rhum.2015.nc'
'sdfopen D:\lw\zl\air.2015.nc'
'enable print D:\lw\pm\rhumsd.gmf'
'set grads off'
'set grid off'
'set csmooth on'
'set parea 2.9 9.9 1.5 8'
'set lon 107.3'
'set lat 37.8'
'set zlog on'
'set xlabs 0908|0920|1008|1020|1108|1120|1208|1220'
lev1=1000
lev2=300
t1=880
t2=896
while (t1 <= t2)
while (lev2 <= lev1)
  'set lev 'lev1''
  'define p=lev'
  'define e=(6.112*exp((17.67*(air.4-273.16))/(air.4-29.65)))'
  'define qs=(0.622*e/(p-0.378*e))'
  'define q=qs*rhum.3/100'
  'define qdiv=hdivg(q*uwnd.1/9.8,q*vwnd.2/9.8)'
  'd qdiv'
  lev1=lev1-1
endwhile
t1=t1+1
endwhile
'print'
'disable print'
'reinit'
然而。。。错误太多把一开始的命令都顶掉了。。。

而且其实我不太确定我这个思路对不对。。。也想过用'collect 1 gr2stn( , )但是 貌似只能按站点和经纬度储存，强行按时间高度储存不行。。。
想了两天 最后还是来请教大神了