想弄一个水汽通量的时间剖面，结果居然是常数？

瑾珂

'reinit'
'sdfopen L:/DZQDL/nc/erainterim/u_wind.nc'
'sdfopen L:/DZQDL/nc/erainterim/v_wind.nc'
'sdfopen L:/DZQDL/nc/erainterim/specific_humidity2.nc'
'sdfopen L:/DZQDL/nc/erainterim/temp.nc'
'sdfopen L:/DZQDL/nc/erainterim/z.nc'
'set grid off'
'set grads off'
'set csmooth on'



'set t 21 28'
'set lon 96 108'
'set lat 20 30'


*------------计算比湿  单位：kg/kg--------------*
'define es=(6.11*exp((17.2693882*(t.4-273.16))/(t.4-35.86)))'
'define qs=0.62197*es/(z.5-0.378*es)'
*-------计算水汽通量(单层)  单位：(g*s)/kg=g/(s*hpa*cm)------*
'define qu=u.1*q.3*1e3/9.8'
'define qv=v.2*q.3*1e3/9.8'
*------------计算水汽通量散度------*
'define qdiv=hdivg(qu,qv)*1e5'
'set grads off'
'set parea 1 9.0 2 8'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set ylint 2'
'set xlint 2'
'set t 21 28'
*'set lon 102.04'
*'set lat 25.09'
*'set lon 102.10'
*'set lat 24.40'
'set lon 102.4'
'set lat 24.1'
'set gxout contour'
'd mag(qu,qv)'

'set gxout contour'
'set cthick 6'
'set ccolor 1'
'd qdiv'

'C:\grads\Contents\Cygwin\Versions\2.0.1.oga.1\i686\cbar_matlab 5 1 0 5'
*'run axis.gs -type b -position i -label on -interval 3 -sinterval 2 -lsize 0.2  -voffset -0.2'
*'run axis.gs -type l -position i -label on -interval 100 -lsize 0.2'
'printim l:/DZQDL/2014datafiles/divergence_time/div_t2_sanchahe.wmf white'
'disable gxprint'
*

某某晨

我也遇到这个问题了，请问楼主解决了吗