- 积分
- 3274
- 贡献
-
- 精华
- 在线时间
- 小时
- 注册时间
- 2012-12-1
- 最后登录
- 1970-1-1
|
数值模式
使用语言: |
ncl |
我看过提问的智慧: |
看过 |
自己思考时长(天): |
2 |
问题概况: |
在计算水汽通量时用到插值函数wrf_user_intrp3d,函数垂直坐标选择气压,不知道所用变量的气压值怎样输出 |
问题截图: |
- |
系统平台: |
linux |
登录后查看更多精彩内容~
您需要 登录 才可以下载或查看,没有帐号?立即注册
x
wrf_user_intrp3d(var3d,vert, 等),其中vert要和var3d的垂直层数对应,然而如何得到变量的这些层的值呢?
求助各位大神!
以下是读取和插值时的不正确ncl代码,我主要问的内容就是如何得到pressure_levels的值
do it=0,ntimes-1
tk = wrf_user_getvar(a,"tk",it) ; T in K
w = wrf_user_getvar(a,"wa",it) ; w averaged to mass points
p = wrf_user_getvar(a,"pressure",it) ;
rh = wrf_user_getvar(a,"rh",it) ; relative humidity;
z = wrf_user_getvar(a,"height",it)
pressure_levels = (/500./)
nlevels = dimsizes(pressure_levels) ; number of pressure levels
tk_plane = wrf_user_intrp3d(tk,p,"h",pressure_levels,0.,False)
w_plane = wrf_user_intrp3d( w,p,"h",pressure_levels,0.,False)
p_plane = wrf_user_intrp3d( p,p,"h",pressure_levels,0.,False)
z_plane = wrf_user_intrp3d( z,p,"h",pressure_levels,0.,False)
rh_plane = wrf_user_intrp3d(rh,p,"h",pressure_levels,0.,False)
e = p_plane*rh_plane/(0.622+rh_plane)
mm=e*w_plane*0.622/(287*tk_plane)
opts=True
opts@MainTile="lalala"
contour = wrf_contour(a,wks,mm(it,:,:),opts)
plot=wrf_map_overlays(a,wks,(/contour/),pltres,mpres)
end do
|
|