水汽通量及水汽通量散度程序

wuxiao

'reinit'
'sdfopen e:\atmodate\pressure\air.2009.nc'
'sdfopen e:\atmodate\pressure\rhum.2009.nc'
'sdfopen e:\atmodate\pressure\uwnd.2009.nc'
'sdfopen e:\atmodate\pressure\vwnd.2009.nc'
'set mpdset cnworld'
'set lon 70 150'
'set lat 15 70'
'set lev 925'
t=302
while(t<=303)
'set t 't''
'set grads off'
'define p=lev'
'define e=exp(17.67*(air-273.15)/(air-29.65))*(rhum.2/100)*6.112'
'define q=622*e/(p-0.378*e)'
'define q1=uwnd.3*q/9.8'
'define q2=vwnd.4*q/9.8'
'set gxout shaded'
'set cint 5'
'd mag(q1,q2)'
'set gxout barb'
'set ccolor 2'
'd uwnd.3*2.5;vwnd.4*2.5'
'run cbarn.gs'
t=t+1
endwhile
'printim C:\Users\Administrator\Desktop\chutu\d.png white x1000 y800'

水汽通量散度
'reinit'
'sdfopen e:/atmodate/pressure/uwnd.2009.nc'
'sdfopen e:/atmodate/pressure/vwnd.2009.nc'
'sdfopen e:/atmodate/pressure/rhum.2009.nc'
'sdfopen e:/atmodate/pressure/air.2009.nc'
'set mpdset cnworld'
t=304
'set lev 925'
'set t 't
'set lon 70 150'
'set lat 15 70'
'set gxout shaded'
'define p=lev'
'define e=(6.112*exp((17.67*(air.4-273.16))/(air.4-29.65)))'
'define qs=(0.622*e/(p-0.378*e))'
'define q=qs*rhum.3/100'
'define qdiv=hdivg(q*uwnd.1/9.8,q*vwnd.2/9.8)'
'd qdiv'
'run cbarn.gs'
'printim e:/gs/qdiv.png white x1000 y800'

linda2520

有没有垂直整层的水汽通量散度剖面图呢，随时间变化的

qlguai

楼主你好，你的求水汽通量的gs文件，水汽通量是一个矢量对吧，你只要gxout barb 就好了啊，后面的这一对程序是什么意思啊？为什么还要输出阴影图？
'set gxout shaded'
'set cint 5'
'd mag(q1,q2)'

叶の幽灵

qlguai 发表于 2015-8-30 21:21
楼主你好，你的求水汽通量的gs文件，水汽通量是一个矢量对吧，你只要gxout barb 就好了啊，后面的这一对程 ...

为了直观显示值的大小吧

fyy123

请问rhum是什么？

fyy123

fyy123 发表于 2015-11-2 20:39
请问rhum是什么？

是grib2资料中的rhprs吗？

asu

为什么不是shum？

寻梦

怎么在程序运行的时候一直显示有错呢

小冰梨

{:eb502:}

karenlk

{:5_235:}