grads画Q矢量与风场叠加出问题了，请教

寻路人

'reinit'
'sdfopen e:\hgt.2014.nc'
'sdfopen e:\uwnd.2014.nc'
'sdfopen e:\vwnd.2014.nc'
'set mpdset cnworld'
'set lon 50 140'
'set lat 5 70'
'set t 146'
'define hgtc=hgt.1'
'define u=uwnd.2'
'define v=vwnd.3'
'define ug= -9.8*cdiff(hgtc,y)/(111177*cdiff(lat,y))'
'define vg= 9.8*cdiff(hgtc,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dugdx=cdiff(ug,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dvgdx=cdiff(vg,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dtdx=cdiff(hgtc,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dugdy= cdiff(ug,y)/(111177*cdiff(lat,y))'
'define dvgdy= cdiff(vg,y)/(111177*cdiff(lat,y))'
'define dtdy= cdiff(hgtc,y)/(111177*cdiff(lat,y))'
'define qx=(-287.05/lev)*(dugdx*dtdx+dvgdx*dtdy)*100000000000000'
'define qy=(-287.05/lev)*(dugdy*dtdx+dvgdy*dtdy)*100000000000000'
'define div=hdivg(qx,qy)*100000'
'set mpdset cnworld'
'set grads off'
'set grid off'
'set lev 700'
'set lon 90 130'
'set lat 10 50'
'set t 146'
'd div'
'set gxout vector'
'set ccolor 1'
'd u;v'
'printim e:\Q700hpa0614.gif gif white x1200 y800'
;

youyou

你用下面的GS试试，我这里用了一下是能出图的，希望于你有用！
'reinit'
'sdfopen D:\2014shuju\hgt.2014.nc'
'sdfopen D:\2014shuju\uwnd.2014.nc'
'sdfopen D:\2014shuju\vwnd.2014.nc'

'define hgtc=hgt.1'
'define u=uwnd.2'
'define v=vwnd.3'
'define ug= -9.8*cdiff(hgtc,y)/(111177*cdiff(lat,y))'
'define vg= 9.8*cdiff(hgtc,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dugdx=cdiff(ug,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dvgdx=cdiff(vg,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dtdx=cdiff(hgtc,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))'
'define dugdy= cdiff(ug,y)/(111177*cdiff(lat,y))'
'define dvgdy= cdiff(vg,y)/(111177*cdiff(lat,y))'
'define dtdy= cdiff(hgtc,y)/(111177*cdiff(lat,y))'
'define qx=(-287.05/lev)*(dugdx*dtdx+dvgdx*dtdy)*100000000000000'
'define qy=(-287.05/lev)*(dugdy*dtdx+dvgdy*dtdy)*100000000000000'
'define div=hdivg(qx,qy)*100000'

'set mpdset cnworld'
'set map 15 4 15'

'set grads off'
'set grid off'

'set lev 700'
'set lon 70 140'
'set lat 15 55'
'set t 146'

'set gxout contour'
'set ccolor 2'
'set cthick 10'
'set clab forced'
'set clskip 1'
'd div'

'set gxout barb'
'set ccolor 1'
'd 2.5*u;2.5*v'

'set mpdset cnriver'
'set map  3  1 8'
'draw map'

'draw title 26MAY2014_700hPa_Q'
'printim f:\Q700hpa0614.png white x1200 y800'

寻路人

上面求Q矢量的脚本是引用本论坛里某位同志上传的，在此表示感谢