求助！求助！在论坛里找到了一个画整层水汽通量的gs，不知道自己画出来的是什么变量

林国英

'reinit'
'sdfopen E:\uwind\uwnd.2015.nc'*纬向风
'sdfopen E:\vwind\vwnd.2015.nc'*经向风
'sdfopen E:\shum\shum.2015.nc'*相对湿度
'enable print e:\sy1\ma.gmf'
i=693
while(i<=788)
'set lon 50 170'
'set lat 0 80'
'set z 8'
'define dqx=cdiff(shum.3,x)'
'define dqy=cdiff(shum.3,y)'
'define dx=cdiff(lon,x)*3.1416/180'
'define dy=cdiff(lat,y)*3.1416/180'
'define sq=-(1.0/9.8))*((uwnd*dqx)/(cos(lat*3.1416/180)*dx)+vwnd.2*dqy/dy)/6.37E6'
'define uwind=uwnd*shum.3'
'define vwind=vwnd.2*shum.3'
'c'
i=i+1
endwhile
'define pj=ave(sq,t=693,t=788,1)'
'set lon 60 160'
'set lat 10 60'
'set z 9'
'set grid off'
'set grads off'
'set mpdset cnriver'
'set gxout shaded'
'd pj*1000000000'
'cbarn'
'define u=ave(uwind,t=693,t=788,1)'
'define v=ave(vwind,t=693,t=788,1)'
'set gxout vectour'
'd u;v'
'draw title Moisture Advection'
'print'
'disable print'
;
*****************************************************************************************************
嫌我自己改的啰嗦的，这里呈上原程序
'reinit'
'sdfopen E:\data\reanalysisdata\u.nc'
'sdfopen E:\data\reanalysisdata\v.nc'
'sdfopen E:\data\reanalysisdata\sh.nc'
'set gxout shaded'
'set lon 105 145'
'set lat 10 40'
'set t 23'
'set z 9'
'define dqx=cdiff(q.3,x)'
'define dqy=cdiff(q.3,y)'
'define dx=cdiff(lon,x)*3.1416/180'
'define dy=cdiff(lat,y)*3.1416/180'
'd -(1.0/9.8))*((u.1*dqx)/(cos(lat*3.1416/180)*dx)+v.2*dqy/dy)/6.37E6'
'cbarn 0.5 1'
'set gxout vector'
*'set gxout stream'
*'define u=u.1*q.3'
*'define v=v.2*q.3'
d skip(u.1*q.3,4);skip(v.2*q.3,4)'
*'d u;v'
'draw title Moisture Advection'
'printim E:\picture\MoistureAdvection23,9.png white'
;

爱地理不爱物理

第八高度层上的该时段内的平均水汽通量和水汽平流

林国英

爱地理不爱物理 发表于 2016-5-11 09:06
第八高度层上的该时段内的平均水汽通量和水汽平流

那我想要整层水汽的话是set z 1 8吗？

爱地理不爱物理

林国英 发表于 2016-5-11 20:57
那我想要整层水汽的话是set z 1 8吗？

看到你问这种问题我觉得我三言两语也解释不清，就算解释了你估计也听不懂，建议你好好看看整层水汽通量的定义以及计算公式，百度或者看专业书或者课件，或者问同学问老师，很多资源可以利用。呵呵呵～～～