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