Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(out) | :: | output | |||
real, | intent(in) | :: | ustar |
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