Q矢量散度---修改时间或者高度出的图都是一样的T_T

Diary

850  15日20时

850 16日20时

在论坛上找到的程序  可是图真的都一样

Reinit

clear

sdfopen e:\nc\g\hgt.2014.nc

define hgtc=hgt.1

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 850

set lon 90 120

set lat 30 45

set t 912

set gxout contour

set ccolor 2

set cthick 10

set clab forced

set clskip 1

d div

set mpdset cnriver

set map  3  1 8

draw map

printim f:\Q850hpa081620.png white x1200 y800

请大家帮忙~