请帮看下问题出在哪里？ 湿Q矢量的

Diary

reinit
sdfopen e:\nc\g\hgt.2014.nc
sdfopen e:\nc\g\uwnd.2014.nc
sdfopen e:\nc\g\vwnd.2014.nc

set lev 1000 10
set t 900 920
set lon 0 360
set lat -90 90

set grads off
set grid off
set mpdset cnworld

define hgtc=hgt.1
define u=uwnd.2
define v=vwnd.3

define f=2*7.292*sin(lat*3.14159/180)

define dudx=cdiff(u,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))
define dvdx=cdiff(v,x)/(111177*cdiff(lon,x)*cos(lat*3.1415/180))

define dudy= cdiff(u,y)/(111177*cdiff(lat,y))
define dvdy= cdiff(v,y)/(111177*cdiff(lat,y))

define dudp= u(z-1)-u(z+1)/(100*(lev(z-1)-lev(z+1)))
define dvdp= v(z-1)-v(z+1)/(100*(lev(z-1)-lev(z+1)))

define qx=f*(dudx*dvdp-dvdx*dudp)
define qy=f*(dudy*dvdp-dvdy*dudp)
define qdiv=hdivg(qx,qy)*1e6

set lon 90 120
set lat 30 45
set t 905
set lev 850

set gxout contour
set ccolor 1
set cint 2
set clskip 1
set clopts -1 -1 0.11

d qdiv

set mpdset cnriver
set map  3  1 8
draw map

printim f:\Q.png white x1200 y800

湿Q矢量
上面的程序是根据论坛上的未完成程序自己修改的
没有提示有错误

问题1.修改维度 不出图 或者出一点图
2.不知道dudp dvdp的表达方式对不对
3.实际公式后面还有两项 原程序没有加 不知道加不加上
图中最下面两个公式是要求的 实际程序中只有f*(*****)的部分

谢谢大家

Diary

自己顶顶顶顶顶顶顶顶顶

earryzrp

二楼给出了湿q矢量的表达式，你的gs文件中少了两项，一个是和位温有关的那项，还有非绝热加热项