比湿。水汽通量，散度

nicocle

哪位大神能帮我写写，比湿，水汽通量，水汽通量散度的FNL资料计算公式？写在GS文件的那种？？

东风急流

'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'
tt=302
while(tt<=303)
'set t 'tt''
'set grads off'
'set grid 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)'
*'d q1;q2'
'set gxout barb'
'set ccolor 2'
'd uwnd.3*2.5;vwnd.4*2.5'
'run cbarn.gs'
tt=tt+1
endwhile
'disable print'
'printim C:\Users\Administrator\Desktop\chutu\d.png white x1000 y800'

nicocle

我的是FNL中的GRIB2资料

zhuzhejl

用GRADS提取后FORTRAN写

nicocle

zhuzhejl 发表于 2015-5-10 08:02
用GRADS提取后FORTRAN写

不是可以用GS文件算吗？

东风急流

'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'

东风急流

上面两个不是fnl资料的，但可以照写，我也正在画，大家一起交流吧

qq88521113

看来你对grads的研究还不够深入，在学习的间隙应该适当的看一下相关资料，比如摇一把{:5_231:}{:5_231:}

Kazoo

站点降水grads做图，批量fnl数据生成ctl并画图，暴雨中尺度分析相关物理量gs
http://bbs.06climate.com/forum.p ... 4&fromuid=50724
(出处: 气象家园)

简_单998

学习了，谢谢