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      subroutine free_disc_pumping(node1,node2,nodem,nelem,lakon,kon,
20     &        ipkon,nactdog,identity,ielprop,prop,iflag,v,xflow,f,
21     &        nodef,idirf,df,cp,r,physcon,dvi,numf,set,shcon,
22     &        nshcon,rhcon,nrhcon,ntmat_,co,vold,mi,ttime,time,
23     &        iaxial,iplausi)
24!
25!     Free-Disc-Pump-Flow
26!
27      implicit none
28!
29      logical identity
30      character*8 lakon(*)
31      character*81 set(*)
32!
33      integer nelem,nactdog(0:3,*),node1,node2,nodem,numf,
34     &     ielprop(*),nodef(*),idirf(*),index,iflag,iaxial,
35     &     ipkon(*),kon(*),mi(*),nrhcon(*),ntmat_,nshcon(*),iplausi
36!
37      real*8 prop(*),v(0:mi(2),*),xflow,f,df(*),ttime,time,
38     &     physcon(*),co(3,*),vold(0:mi(2),*),r,cp,dvi,
39     &     shcon(0:3,ntmat_,*),rhcon(0:1,ntmat_,*)
40!
41!
42!
43      if (iflag.eq.0) then
44         identity=.true.
45!
46         if(nactdog(2,node1).ne.0)then
47            identity=.false.
48         elseif(nactdog(2,node2).ne.0)then
49            identity=.false.
50         elseif(nactdog(1,nodem).ne.0)then
51            identity=.false.
52         endif
53!
54      elseif ((iflag.eq.1).or.(iflag.eq.2).or.(iflag.eq.3))then
55!
56!    User defined flow element
57!
58      endif
59!
60      return
61      end
62
63
64