blthickness Subroutine

public subroutine blthickness(output, ustar)

Uses

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

Arguments

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

Source Code

  subroutine blthickness(output,ustar)

    use modglobal, only : numol,Uinf
    implicit none

       real, intent(in)                   :: ustar      ! friction velocity
       real, intent(out)                  :: output     !< BL thickness based on law of the wake
!       real, dimension(kb:ke)             :: mthick
!       real                               :: ucrit
!       real                               :: umax
!       integer :: k
       real :: B     = 5.0       ! Wake parameter
       real :: C     = 0.5       ! Coles parameter
       real :: kappa = 0.41      ! Von k�r�n constant
       real :: lam               ! = Uinf/ustar

       lam = Uinf / ustar
       output = (lam*numol/Uinf) * exp( kappa * (lam - B) - 2.*C)      ! See App. Lund et al.

   end subroutine blthickness