pk_cdi_rl.f File Reference

Go to the source code of this file.

Functions/Subroutines
subroutine	pk_cdi_rl (lqd, rqd, reynolds, beta, cdi_rl)

Function/Subroutine Documentation

pk_cdi_rl()

subroutine pk_cdi_rl	(	real*8	lqd,
		real*8	rqd,
		real*8	reynolds,
		real*8	beta,
		real*8	cdi_rl
	)

!
       implicit none
!
      real*8 lqd,rqd,reynolds,beta,cdi_rl,rqd_cor,lrqd,cdi_r,glrqd
!
      rqd_cor=rqd
!
      if (rqd_cor.gt.lqd) then
         rqd_cor=lqd
      endif
!
      lrqd=lqd-rqd_cor
!      
      call pk_cdi_r(rqd_cor,reynolds,beta,cdi_r)
!      
      glrqd=(1d0+1.298d0*exp(-1.593d0*lrqd**2.33d0))
     &     *(0.435d0+0.021d0*lrqd)/(2.298d0*0.435d0)
!
      cdi_rl=1.d0-glrqd*(1.d0-cdi_r)
!
      return
!