1!
2!     CalculiX - A 3-dimensional finite element program
3!     Copyright (C) 1998-2021 Guido Dhondt
4!
5!     This program is free software; you can redistribute it and/or
6!     modify it under the terms of the GNU General Public License as
7!     published by the Free Software Foundation(version 2);
8!
9!
10!     This program is distributed in the hope that it will be useful,
11!     but WITHOUT ANY WARRANTY; without even the implied warranty of
12!     MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13!     GNU General Public License for more details.
14!
15!     You should have received a copy of the GNU General Public License
16!     along with this program; if not, write to the Free Software
17!     Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
18!
19!cd incompressible for long orifices (eq.6)
20!
21!     author: Yannick Muller
22!
23       subroutine pk_cdi_rl(lqd,rqd,reynolds,beta,cdi_rl)
24!
25       implicit none
26!
27      real*8 lqd,rqd,reynolds,beta,cdi_rl,rqd_cor,lrqd,cdi_r,glrqd
28!
29      rqd_cor=rqd
30!
31      if (rqd_cor.gt.lqd) then
32         rqd_cor=lqd
33      endif
34!
35      lrqd=lqd-rqd_cor
36!
37      call pk_cdi_r(rqd_cor,reynolds,beta,cdi_r)
38!
39      glrqd=(1d0+1.298d0*exp(-1.593d0*lrqd**2.33d0))
40     &     *(0.435d0+0.021d0*lrqd)/(2.298d0*0.435d0)
41!
42      cdi_rl=1.d0-glrqd*(1.d0-cdi_r)
43!
44      return
45!
46      end
47