Mapping vertical levels to physical altitude

Topics related to running an idealized case
Post Reply
mfryba
Posts: 2
Joined: Thu Dec 31, 2020 6:09 pm

Mapping vertical levels to physical altitude

Post by mfryba » Thu Dec 31, 2020 6:20 pm

Hello, I'm not a weather person per se (physicist) and am trying to do some basic WRF runs to generate "realistic" radar returns for a separate simulation. I did find the vapor_wrf Python code to compute dBZ. I ran the basic em_quarter_ss case to get something to start with. I understand the X and Y grids are Cartesian and DX and DY give the spacing. The Z coordinate is a bit tougher since I'm not well versed in pressure level jargon. To generate simulated 3-D radar returns I need to map those Z coordinates to real altitude. It appears some algorithm using DZBOT, DZSTRETCH_? and perhaps some other variables is required. Any tips? Thanks!

Ming Chen
Posts: 1455
Joined: Mon Apr 23, 2018 9:42 pm

Re: Mapping vertical levels to physical altitude

Post by Ming Chen » Tue Jan 05, 2021 5:31 pm

Hi,
DZBOT, DZSTRETCH etc. are designed for setting vertical eta levels. For ideal cases, vertical levels are calculate based on ZTOP, which corresponds to real attitude.
In wrfout, you can find variables PH and PHB, (PH + PHB)/9.8 will give you the real altitude of each model level. You can do vertical interpolation to put the variables in any levels you want for dBZ calculation.
WRF Help Desk

mfryba
Posts: 2
Joined: Thu Dec 31, 2020 6:09 pm

Re: Mapping vertical levels to physical altitude

Post by mfryba » Tue Jan 05, 2021 6:53 pm

Ming Chen wrote:
Tue Jan 05, 2021 5:31 pm
Hi,
DZBOT, DZSTRETCH etc. are designed for setting vertical eta levels. For ideal cases, vertical levels are calculate based on ZTOP, which corresponds to real attitude.
In wrfout, you can find variables PH and PHB, (PH + PHB)/9.8 will give you the real altitude of each model level. You can do vertical interpolation to put the variables in any levels you want for dBZ calculation.
Thanks! I started reviewing the Technical Note that has some details; appreciate the quick summary because it is clear that once I potentially get away from an ideal case the mapping of eta levels to altitude is non-trivial, and even the specification of eta levels and how that maps to nominal dry pressure is a bit hairy. Curious as to units; I see PH and PHB have units m2 s-2; what is the units of that 9.8 factor and how does that relate to physical constants? Oh, I just found on page 5.142 of the Users Guide that same equation, and the 9.81 is clearly just acceleration due to gravity g, which has units of m s-2, so that makes sense that the result is in meters. :D

Ming Chen
Posts: 1455
Joined: Mon Apr 23, 2018 9:42 pm

Re: Mapping vertical levels to physical altitude

Post by Ming Chen » Wed Jan 20, 2021 11:02 pm

PH and PHB are geopotential, which is a physical variable defined in geoscience. You can easily find its meaning in the website.
WRF Help Desk

Post Reply

Return to “Idealized Cases”