momentumthicknessmo Subroutine

public subroutine momentumthicknessmo(output, ustar, blth, lmo)

Uses

  • proc~~momentumthicknessmo~~UsesGraph proc~momentumthicknessmo momentumthicknessmo module~modglobal modglobal proc~momentumthicknessmo->module~modglobal

Arguments

Type IntentOptional Attributes Name
real, intent(out) :: output
real, intent(in) :: ustar
real, intent(in) :: blth
real, intent(in) :: lmo

Source Code

  subroutine momentumthicknessmo(output,ustar,blth,lmo)

    use modglobal, only : pi,Uinf
!    use modinletdata, only : ubulk
 !   use modmpi, only    : myid
    implicit none

    real, intent(in)                   :: ustar      ! friction velocity
    real, intent(in)                   :: lmo        ! Obukhov length
    real, intent(in)                   :: blth       ! boundary layer thickness
    real, intent(out)                  :: output     ! momentum thickness
    real :: B     = 5.0       ! Wake parameter
    real :: C     = 0.5       ! Coles parameter
    real :: kappa = 0.41      ! Von k�r�n constant
    real :: cmo   = 0.702     ! constant in MO theory (0.135*5.2)
    real :: lam               ! = Uinf/ustar

    lam = Uinf / ustar
    output = (1. + C + 0.5*cmo*blth/lmo)/(kappa*lam) - (1./(( kappa**2)*(lam**2)))*(2. + 2.*C*(1.852/pi +1.) + (3./2.)*(C**2) + (blth-0.25)*2.*cmo/lmo + (1. + 4./pi)*blth*C*cmo/lmo + (1./6.)*((cmo/lmo)**2)*(blth**2) )* blth

  end subroutine momentumthicknessmo