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. height and surface
     real(r8), intent(in) :: zldis ! reference height "minus" zero displacement heght [m]
     real(r8), intent(in) :: z0m   ! roughness length, momentum [m]
 
     real(r8), intent(out) :: um   ! wind speed including the stability effect [m/s]
     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