Type | Intent | Optional | Attributes | Name | ||
---|---|---|---|---|---|---|
real, | intent(out) | :: | output | |||
real, | intent(in) | :: | ustar | |||
real, | intent(in) | :: | blth |
subroutine momentumthickness(output,ustar,blth)
use modglobal, only : pi,Uinf
! use modinletdata, only : ubulk
! use modmpi, only : myid
implicit none
real, intent(in) :: ustar ! friction velocity
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 :: lam ! = Uinf/ustar
lam = Uinf / ustar
output = ((1. + C)/(kappa*lam) - (1./(( kappa**2)*(lam**2)))*(2. + 2.*C*(1.852/pi +1.) + (3./2.)*(C**2)))* blth
end subroutine momentumthickness