关于Grads中Vint函数及整层水汽通量的计算单位量级问题

六月薇雨

学习学习

Jinnee

厉害

1048651077

谢谢，很有用

Amandaaa

学习学习，困扰好久

瑾珂

到底有没有g

栤唥de嶶笑

瑾珂 发表于 2022-11-28 09:03
到底有没有g

整层的没有g，单层的有g

瑾珂

我计算整层的水汽通量总是出错，以下是我的gs，老师能不能帮忙看看是哪里出问题？

'reinit'
'open J:\0813-032\DJBY\20y_winter\2017_1\fnl_201701.ctl'


i=1
*while(i<=21)
'set t 'i''
'set gxout shaded'
'set mpdset cnworld'
'set map 1 1 9'
'set csmooth on'
'set grads off'
'set grid off'


'set lev 1000 1'
'define tp=TMPprs'
'define rh=RHprs'
'define prs=lev'
'define uw=UGRDprs'
'define vw=VGRDprs'

*------------计算比湿  单位：kg/kg--------------*
'define es=(6.11*exp((17.2693882*(tp-273.16))/(tp-35.86)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define qq=rh*qs/100'
*-------计算水汽通量(单层)  单位：(g*s)/kg=g/(s*hpa*cm)------*
'define qu=uw*qq*1000'
'define qv=vw*qq*1000'

'define pres=PRESsfc/100'
'define uqx=vint(pres,qu,300)'
'define vqx=vint(pres,qv,300)'

*------------计算水汽通量散度------*

'define aqdiv=hdivg(uqx,vqx)*100000'

'set grads off'
'set grid off'
'set parea 1 9.5 0.5 8'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set xlint 5'
'set ylint 2'
'set gxout shaded'
'set cint 2'
'd mag(uqx,vqx)'
'run J:\0813-032\2022222snow\gs\output.gs'
'set gxout barb'
'set ccolor 2'
'd uqx;vqx'
'set gxout contour'
'set ccolor 9'
'set cthick 6'
'd aqdiv'
*'cbarm'
*'D:\grads\OpenGrADS2.1.a2\opengrads\Contents\Cygwin\Versions\2.1.a2.oga.1\i686\cbarm 1 1 1'
'D:\grads\OpenGrADS2.1.a2\opengrads\Contents\Cygwin\Versions\2.1.a2.oga.1\i686\cbar_interp1'


'q dim'
it=sublin(result,5)
it=subwrd(it,6)
'draw title 'it'-jifen-waterflux'
'printim  J:\0813-032\DJBY\ca_saf\201701\waterflu\'it'-jifen-waflux.wmf white'
*c
*i=i+1
*endwhile
'disable gxprint'
;

栤唥de嶶笑

瑾珂 发表于 2022-11-28 10:09
我计算整层的水汽通量总是出错，以下是我的gs，老师能不能帮忙看看是哪里出问题？

'reinit'

aqu 不是变量和函数名，请检查脚本

瑾珂

栤唥de嶶笑 发表于 2022-11-28 11:41
aqu 不是变量和函数名，请检查脚本

老师，后面我把‘set lev 1000 1’删掉后能运行了，就是地图又不见了，也不知道我的这个单位换算对不对，感觉量级大
'reinit'
'open J:\0813-032\DJBY\20y_winter\2017_1\fnl_201701.ctl'


i=1
*while(i<=21)
'set t 'i''



*'set lev 1000 1'
'define tp=TMPprs'
'define rh=RHprs'
'define prs=lev'
'define uw=UGRDprs'
'define vw=VGRDprs'

*------------计算比湿  单位：kg/kg--------------*
'define es=(6.11*exp((17.2693882*(tp-273.16))/(tp-35.86)))'
'define qs=0.62197*es/(prs-0.378*es)'
'define qq=rh*qs/100'
*-------计算水汽通量(单层)  单位：(g*s)/kg=g/(s*hpa*cm)------*
'define qu=uw*qq*1000/9.8'
'define qv=vw*qq*1000/9.8'

'define pres=PRESsfc/100'
'define uqx=vint(pres,qu,300)'
'define vqx=vint(pres,qv,300)'

*------------计算水汽通量散度------*

'define aqdiv=hdivg(uqx,vqx)*100000'

'set lat 15 35'
'set lon 90 120'
'set gxout shaded'
'set mpdset cnworld'
'set map 1 1 9'
'set csmooth on'
'set grads off'
'set grid off'
'set parea 1 9.5 0.5 8'
'set xlopts 1 6 0.2'
'set ylopts 1 6 0.2'
'set xlint 5'
'set ylint 2'
'set gxout shaded'
'set cint 1'
'd mag(uqx,vqx)'
'run J:\0813-032\2022222snow\gs\output.gs'
'set gxout barb'
'set ccolor 2'
'd uqx;vqx'
'set gxout contour'
'set ccolor 9'
'set cthick 6'
'd aqdiv'
*'cbarm'
*'D:\grads\OpenGrADS2.1.a2\opengrads\Contents\Cygwin\Versions\2.1.a2.oga.1\i686\cbarm 1 1 1'
'D:\grads\OpenGrADS2.1.a2\opengrads\Contents\Cygwin\Versions\2.1.a2.oga.1\i686\cbar_interp1'


'q dim'
it=sublin(result,5)
it=subwrd(it,6)
'draw title 'it'-jifen-waterflux'
'printim  J:\0813-032\DJBY\ca_saf\201701\waterflu\'it'-jifen-waflux.wmf white'
*c
*i=i+1
*endwhile
'disable gxprint'
;

逝者如斯

，学习学习！