爱气象,爱气象家园! 

气象家园

 找回密码
 立即注册

QQ登录

只需一步,快速开始

搜索
查看: 3753|回复: 0

请问wrf的结果中如何输出bottom_up_stag的值

[复制链接]
发表于 2015-10-10 21:19:19 | 显示全部楼层 |阅读模式
数值模式
使用语言: 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

密码修改失败请联系微信:mofangbao
您需要登录后才可以回帖 登录 | 立即注册

本版积分规则

Copyright ©2011-2014 bbs.06climate.com All Rights Reserved.  Powered by Discuz! (京ICP-10201084)

本站信息均由会员发表,不代表气象家园立场,禁止在本站发表与国家法律相抵触言论

快速回复 返回顶部 返回列表