WRF------------单层水汽通量

`;--------------------------------WRFOUT D01

a             = addfile(".nc","r")

u                    = wrf_user_getvar(a,"ua",-1)     
v                    = wrf_user_getvar(a,"va",-1)        
p                    = wrf_user_getvar(a,"pressure",-1)  ;
z                    = wrf_user_getvar(a,"height",-1)
qvapor               = wrf_user_getvar(a,"QVAPOR",-1)    ; Qv

qv       = qvapor*1000.
qv@units = "g/kg"
g    = 9.80665          ; m/s2
g@units = "m/s2"
u_500   = wrf_user_intrp3d(u,p,"h",500.,0.,False)
v_500   = wrf_user_intrp3d(v,p,"h",500.,0.,False)
qv_500  = wrf_user_intrp3d(qv,p,"h",500.,0.,False)
spd = (u_500*u_500 + v_500*v_500)^(0.5)            
spd@description = "Wind Speed"
spd@units = "m/s"

;单层水汽通量 flux = 比湿(g/kg)*风V(m/s)/g(m/s2)

flux_5       = qv_500*spd/g     ;g*(s-1*hPa-1*cm-1)

flux_500     = flux_5*10        ;10-1*g*(s-1*hPa-1*cm-1)

avg_flux_500 =dim_avg_n_Wrap(flux_500,0)

avg_u        =dim_avg_n(u_500,0)

avg_v        =dim_avg_n(v_500,0)

copy_VarCoords(qv_500(0,:,:),avg_flux_500

copy_VarCoords(u_500(0,:,:),avg_u)

copy_VarCoords(v_500(0,:,:),avg_v)

delete([/flux_5,flux_500,spd,qv_500,u_500,v_500/])

;--------------换经纬度
lat2d = a->XLAT(0,:,:)
lon2d = a->XLONG(0,:,:)
lat1d = lat2d(:,0)
lon1d = lon2d(0,:)

grid_flux_500 = rcm2rgrid_Wrap(lat2d,lon2d,avg_flux_500,lat1d,lon1d,0)
grid_u        = rcm2rgrid_Wrap(lat2d,lon2d,avg_u,lat1d,lon1d,0)
grid_v        = rcm2rgrid_Wrap(lat2d,lon2d,avg_v,lat1d,lon1d,0)

delete([/avg_flux_500,avg_u,avg_v/])`

WRF------------单层水汽通量

  • 8
    点赞
  • 10
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值