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