[Wrf-users] compute altitudes in WRF

[Alba Badia (ENV)]

Hi everyone,

I am starting with the WRF-Chem. I just want to do a simple plot: Verticals profiles of O3.

To be sure, is this the way to compute the altitudes (height above sea level) in WRF?

for (k in 1: num_vertical_layers){
altitude[,,k,]=((PH[,,k,]+PH[,,k+1,] )/2 + (PHB[,,k,]+PHB[,,k+1,] )/2 )/G
}

where num_vertical_layers is the "bottom_top", G is the gravity, PHB is the Base State Geopotential and PH is the Perturbation Geopotential.

Thank you in advance,

Alba
-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://mailman.ucar.edu/pipermail/wrf-users/attachments/20150115/00ad824e/attachment.html