<HTML> <BODY BGCOLOR=#ccccdd LINK=#0000aa VLINK=#0000ff ALINK=#ff0000 ><BASE TARGET="bottom_target"><PRE>
<A NAME='DA_GET_HEIGHT'><A href='../../html_code/gen_be/da_get_height.inc.html#DA_GET_HEIGHT' TARGET='top_target'><IMG SRC="../../gif/bar_red.gif" border=0></A>
subroutine da_get_height( input_file, dim1, dim2, dim3, height),4
!---------------------------------------------------------------------------
! Purpose: Calculates T, RH from input WRF file.
!---------------------------------------------------------------------------
implicit none
character(len=200), intent(in) :: input_file ! NETCDF file nane.
integer, intent(in) :: dim1, dim2, dim3 ! Dimensions.
real, intent(out) :: height(1:dim1,1:dim2,1:dim3) ! Height.
character(len=10) :: var ! Variable to search for.
integer :: k ! Loop counter.
real :: gravity_inv ! 1/gravity.
real :: phb(1:dim1,1:dim2) ! Base state geopotential.
real :: ph(1:dim1,1:dim2) ! Perturbation geopotential.
real :: phfull(1:dim1,1:dim2,1:dim3+1) ! Geopotential.
if (trace_use) call da_trace_entry
("da_get_height")
gravity_inv = 1.0 / gravity
do k = 1, dim3+1
var = "PHB"
call da_get_field( input_file, var, 3, dim1, dim2, dim3+1, k, phb)
var = "PH"
call da_get_field( input_file, var, 3, dim1, dim2, dim3+1, k, ph)
phfull(:,:,k) = phb + ph ! Calculate full geopotential on full(w) model levels:
end do
do k = 1, dim3
height(:,:,k) = 0.5 *( phfull(:,:,k+1) + phfull(:,:,k)) * gravity_inv
end do
if (trace_use) call da_trace_exit("da_get_height")
end subroutine da_get_height