整层水汽通量散度GS文件，及阴影图GS文件

zxc

首先得感谢楼主的分享！有些问题的请教楼主和各位大侠：
在计算水汽通量时楼主用了一时间个循环，我不太明白，如下：

#该gs绘出水汽通量是阴影图
'reinit'
'open f:\tt.ctl'
'enable print d:\1243567788sqtl.gmf'
#经纬度设置
'set lat 20 35 '
'set lon 100 120'
'set lev 700'
#先设定时间循环
t=1
while(t<=9)
'c'
'set t 't''
#设置图片角标底图
'set grads off'
'set grid off'
'set mpdset cnworld'
#'set map  15 5 1'
#'draw map'
#设置划线格式
#'set cthick 9'
#'set xlopts 1 4 0.12'
#'set ylopts 1 4 0.12'
#'set cint 0.5'
#★重点：编程水气通量的算法
'define prs=lev'
#饱和水汽压的计算，单位：毫巴
'define es=6.112*exp(17.67*(Tk-273.15)/(Tk-29.65))'
#饱和比湿的计算  单位：克/克
'define qs=0.622*es/(lev-0.278*es)'
#rh是比湿，利用比湿求水汽
'define q=rh*qs/100'
#a为|V|
'define a=sqrt(u*u+v*v)'
#地球引力系数
'define g=9.8'
#以阴影图显示，命令一直到run cbarn.gs
'set gxout shaded'
'set cmin 10'
#10以上才画阴影图
'd a*q*1000/9.8'
#本来是q*|V|/g，但这里q的单位为：克/克，而水汽通量中的q一般取克、千克，所以乘以1000
#显示一下尺度标,画阴影图用
'run cbarn.gs'
#set gxout contour是重新设定图形输出格式为二维等值线图，叠加在阴影上，所以一般先画阴影

，再画等值线，不然阴影覆盖了等值线
'set gxout contour'
'set ccolor 0'
'set cmin 6'
只画值为6以上的等值线
'd a*q*1000/9.8'

'print'
pull dummy
t=t+1
endwhile
'disable print'

这个时间循环是指哪段时间，如果我想计算某一时次的水汽通量，如28日06时，怎样循环？

zxc

首先得感谢楼主的分享！有些问题的请教楼主：
在计算水汽通量时楼主用了一时间个循环，我不太明白，如下：

#该gs绘出水汽通量是阴影图
'reinit'
'open f:\tt.ctl'
'enable print d:\1243567788sqtl.gmf'
#经纬度设置
'set lat 20 35 '
'set lon 100 120'
'set lev 700'
#先设定时间循环
t=1
while(t<=9)
'c'
'set t 't''
#设置图片角标底图
'set grads off'
'set grid off'
'set mpdset cnworld'
#'set map  15 5 1'
#'draw map'
#设置划线格式
#'set cthick 9'
#'set xlopts 1 4 0.12'
#'set ylopts 1 4 0.12'
#'set cint 0.5'
#★重点：编程水气通量的算法
'define prs=lev'
#饱和水汽压的计算，单位：毫巴
'define es=6.112*exp(17.67*(Tk-273.15)/(Tk-29.65))'
#饱和比湿的计算  单位：克/克
'define qs=0.622*es/(lev-0.278*es)'
#rh是比湿，利用比湿求水汽
'define q=rh*qs/100'
#a为|V|
'define a=sqrt(u*u+v*v)'
#地球引力系数
'define g=9.8'
#以阴影图显示，命令一直到run cbarn.gs
'set gxout shaded'
'set cmin 10'
#10以上才画阴影图
'd a*q*1000/9.8'
#本来是q*|V|/g，但这里q的单位为：克/克，而水汽通量中的q一般取克、千克，所以乘以1000
#显示一下尺度标,画阴影图用
'run cbarn.gs'
#set gxout contour是重新设定图形输出格式为二维等值线图，叠加在阴影上，所以一般先画阴影

，再画等值线，不然阴影覆盖了等值线
'set gxout contour'
'set ccolor 0'
'set cmin 6'
只画值为6以上的等值线
'd a*q*1000/9.8'

'print'
pull dummy
t=t+1
endwhile
'disable print'

这个时间循环是指哪段时间，如果我想计算某一时次的水汽通量，如28日06时，怎样循环？

lulu

{:eb511:}{:eb511:}

luoluo969

到处都是热心的高手哦

huanggg

学习下，好东西

道狼52203

xuexi xuexi

huanggg

学习下学习下

766168584

多谢分享这个，找了半天了，呵呵

lulu

谢谢楼主分享。我下载来看了一下，楼主做的是单层的水汽通量散度，那整层大气的该如何求呢？我在用rhum.nc文件等求整层的时候，总会出现问题，布置楼主有没有经验呢？

albert007zh

楼主加油，谢谢分享