flow1.F

Go to the documentation of this file.

Copyright>        OpenRadioss
Copyright>        Copyright (C) 1986-2025 Altair Engineering Inc.
Copyright>
Copyright>        This program is free software: you can redistribute it and/or modify
Copyright>        it under the terms of the GNU Affero General Public License as published by
Copyright>        the Free Software Foundation, either version 3 of the License, or
Copyright>        (at your option) any later version.
Copyright>
Copyright>        This program is distributed in the hope that it will be useful,
Copyright>        but WITHOUT ANY WARRANTY; without even the implied warranty of
Copyright>        MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
Copyright>        GNU Affero General Public License for more details.
Copyright>
Copyright>        You should have received a copy of the GNU Affero General Public License
Copyright>        along with this program.  If not, see <https://www.gnu.org/licenses/>.
Copyright>
Copyright>
Copyright>        Commercial Alternative: Altair Radioss Software
Copyright>
Copyright>        As an alternative to this open-source version, Altair also offers Altair Radioss
Copyright>        software under a commercial license.  Contact Altair to discuss further if the
Copyright>        commercial version may interest you: https://www.altair.com/radioss/.
!||====================================================================
!||    flow1    ../engine/source/fluid/flow1.F
!||--- called by ------------------------------------------------------
!||    resol    ../engine/source/engine/resol.F
!||--- calls      -----------------------------------------------------
!||    daaacc   ../engine/source/fluid/daaacc.F
!||====================================================================
      SUBROUTINE flow1(IFLOW, RFLOW, NBGAUGE, A)
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "flowcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IFLOW(*), NBGAUGE
      my_real rflow(*), a(3,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  N, NEL, NNO, II1, II2, II3, II4, II5, II6
      INTEGER  IR1, IR2, IR3, IR4, IR5, IR6, IR7, IR8
      INTEGER  IADI, IADR, ITYP, JFORM, NDIM, FREESURF
C
      iadi=0
      iadr=0
      ii4 = 0
      ii6 = 0
      DO n=1,nflow
         ityp=iflow(iadi+2)
         IF (ityp==3) THEN
            jform    = iflow(iadi+4)
            nno      = iflow(iadi+5)
            nel      = iflow(iadi+6)
            freesurf = iflow(iadi+25)
            ii1=iadi+1
            ii2=ii1+niflow
            ii3=ii2+nno
            IF(jform==1) THEN
               ii4=ii3+3*nel
               ii5=ii4+nno
               ii6=ii5
               ndim=3
            ELSEIF(jform==2) THEN
               ii4=ii3+5*nel
               ii5=ii4+nno
               ii6=ii5+nbgauge
               ndim=5
            ENDIF
            ir1=iadr+1
            ir2=ir1+nrflow
            ir3=ir2+nel*3
            ir4=ir3+nel*freesurf
            ir5=ir4+nel
            ir6=ir5+nel*freesurf
            ir7=ir6+nel*freesurf
            ir8=ir7+nel*nel
            CALL daaacc(ndim, nno, nel, iflow(ii1), iflow(ii2), iflow(ii3), iflow(ii4), iflow(ii6), 
     .                  a,  rflow(ir2), rflow(ir8))
         ENDIF
         IF(n < nflow) THEN
           iadi=iadi+iflow(iadi+14)  
           iadr=iadr+iflow(iadi+15)
         ENDIF
      ENDDO
C
      RETURN
      END