enthalpythickness Subroutine

public subroutine enthalpythickness(output, tinput, uinput)

Uses

  • proc~~enthalpythickness~~UsesGraph proc~enthalpythickness enthalpythickness module~modglobal modglobal proc~enthalpythickness->module~modglobal module~modinletdata modinletdata proc~enthalpythickness->module~modinletdata module~modmpi modmpi proc~enthalpythickness->module~modmpi module~modsurfdata modsurfdata proc~enthalpythickness->module~modsurfdata mpi mpi module~modmpi->mpi

Arguments

Type IntentOptional Attributes Name
real, intent(out) :: output
real, intent(in), dimension(kb:ke) :: tinput
real, intent(in), dimension(kb:ke) :: uinput

Called by

proc~~enthalpythickness~~CalledByGraph proc~enthalpythickness enthalpythickness proc~inletgen inletgen proc~inletgen->proc~enthalpythickness

Source Code

  subroutine enthalpythickness(output,tinput,uinput)

    use modglobal, only : jb,kb,ke,dzf !,Uinf
    use modinletdata, only : ubulk
    use modsurfdata, only : thls
    use modmpi, only    : myid
    implicit none

       real, dimension(kb:ke), intent(in) :: tinput  !< input temperature
       real, dimension(kb:ke), intent(in) :: uinput  !< input velocity
       real, intent(out)                  :: output  !< momentum thickness
       real, dimension(kb:ke)             :: ethick
       real thlsdummy
       integer :: k

       thlsdummy = thls
       if (tinput(ke) == thls) then
         thlsdummy = thls -0.000001
       end if
       do k=kb,ke
!         ethick(k) = (uinput(k)/uinput(ke)) * ((tinput(k) - tinput(ke)) /(thls - tinput(ke)) )*dzf(k)
         ethick(k) = (uinput(k)/uinput(ke)) * ((tinput(k) - tinput(ke)) /(thlsdummy - tinput(ke)) )*dzf(k)

       end do
       output   = sum(ethick)  ! enthalpy thickness
       if (output==0.) then
         output= 0.000001
       end if

  end subroutine enthalpythickness