整层水汽通量散度计算及其单位

成信院

*------------------------------------------------------
'reinit'
*------------------------------------------------------
'sdfopen C:\Users\i\Desktop\lunwen\fengchang\uwnd.mon.mean.nc'
'sdfopen C:\Users\i\Desktop\lunwen\fengchang\vwnd.mon.mean.nc'
'sdfopen C:\Users\i\Desktop\lunwen\fengchang\shum.mon.mean.nc'
'sdfopen C:\Users\i\Desktop\lunwen\fengchang\hgt.mon.mean.nc'
'set mpdset cnworld'
*------------------------------------------------------
'clear'
'set grads off'
'set lon 70 130'
'set lat 0 40'

nz=1
while(nz<=8)
'set z'nz''
'define qu=uwnd.1*shum.3/9.8'
'd qu'
nz=nz+1
endwhile

nz=1
while(nz<=8)
'set z'nz'
'define qv=vwnd.2*shum.3/9.8'
'd qv'
nz=nz+1
endwhile

nz=1
while(nz<=8)
'set z'nz''
'define qu=uwnd.1*shum.3/9.8'
'define qv=vwnd.2*shum.3/9.8'
'd hdivg(qu,qv)'
nz=nz+1
endwhile

'set z 1'
'define quall=vint(hgt.4,qu,300)*9.8'
'define qvall=vint(hgt.4,qv,300)*9.8'
'd quall'
'd qvall'
'd hdivg(quall,qvall)'


'draw title (c)'
'enable print C:\Users\i\Desktop\lunwen\fengchang\quv.gmf'
'print'
*------------------------------------------------------
'disable print'
*------------------------------------------------------
这是我的gs，但是结果是错的，请楼主帮忙看看哪儿错了。还有我想做的事63,69,79,10年的3-4月的整层水汽通量散度，谢谢楼主

楠楠_

楼主真棒！！！

luoxincheng

很好 顶一个

yujie

好东西啊，谢谢楼主

如梦荼蘼

赞一个~~~抱走了~~~谢谢楼主~~~

不可一世

{:5_213:}

lt9171

为什么图上没有单位显示

在水一方，

{:5_213:}

kc121

网上有很多种方法，不知道哪种计算法才是正确的。

lideqin

很好  谢谢楼主分享！