MoninObukIni

INTERFACE:

   subroutine MoninObukIni (ur, thv, dthv, zldis, z0m, um, obu)

DESCRIPTION:

   Initialization of the Monin-Obukhov length.
   The scheme is based on the work of Zeng et al. (1998):
   Intercomparison of bulk aerodynamic algorithms for the computation
   of sea surface fluxes using TOGA CORE and TAO data. J. Climate,
   Vol. 11, 2628-2644.

USES:

     use clm_varcon, only : grav

ARGUMENTS:

     implicit none
     real(r8), intent(in)  :: ur    ! wind speed at reference height [m/s]
     real(r8), intent(in)  :: thv   ! virtual potential temperature (kelvin)
     real(r8), intent(in)  :: dthv  ! diff of vir. poten. temp. between ref. h
     real(r8), intent(in)  :: zldis ! reference height "minus" zero displaceme
     real(r8), intent(in)  :: z0m   ! roughness length, momentum [m]
     real(r8), intent(out) :: um    ! wind speed including the stability effec
     real(r8), intent(out) :: obu   ! monin-obukhov length (m)

CALLED FROM:

   subroutine BareGroundFluxes in module BareGroundFluxesMod.F90
   subroutine BiogeophysicsLake in module BiogeophysicsLakeMod.F90
   subroutine CanopyFluxes in module CanopyFluxesMod.F90

REVISION HISTORY:

   15 September 1999: Yongjiu Dai; Initial code
   15 December 1999:  Paul Houser and Jon Radakovich; F90 Revision