momentumthickness Subroutine

public subroutine momentumthickness(output, ustar, blth)

Uses

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

Arguments

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

Source Code

  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