grads、程序运行、整层水汽通量散度

1048651077

'reinit'
'sdfopen E:\grads\ERA5\zuvq1.nc'
'sdfopen E:\grads\ERA5\sp1.nc'
'set grid off'  
'set grads off'  

'set lon 60 150'  
'set lat 0 60'  

'set lev 1000 300'
'set t 1 61'
'define uq=u.1*q.1'
'define vq=v.1*q.1'
'define sqtlsd=hdivg(uq,vq)'
'define zsqtlsd=vint(sp.2(lev=0)/100,sqtlsd,300)'
'set lev 500'
'set t 61'*需确定高度和时间
'define zsd2=ave(zsqtlsd,t=1,t=61)*1e4'
define_colors
'set gxout shaded'
'set clevs -10 -8 -6 -4 -2 0 2 4 6 8 10'
'set ccols 47 46 45 44 43 42 41 24 25 26 27 28 29'
'd zsd2'
'run E:\grads\ERA5\cbarn.gs '
'printim E:\grads\ERA5\a.png white x1000 y800'

想问一下各位大神：
①这样计算整层水汽通量散度有没有错误？（能成功运行，并出图）但我不确定出来的图是不是“时间1-61的平均整层水汽通量散度图，还有就是我的量级是1e4，而文献基本是1e5，是哪出了问题？（u，v：m/s；q：kg/kg；sp：Pa）
②grads在运行时，是不是lev=1000→t=1→uq、vq、sqtlsd，接着位势高度不变，t=2→uq、vq、sqtlsd·····t=61；再接着lev=925→t=1→uq、vq、sqtlsd·····依次类推下去？那每一层，每个时次的uq、vq、sqtlsd是否会被独立记录，还是说会被覆盖？
③程序中的”define zsqtlsd···“是否放错位置？