flow1.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "flowcom.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	flow1 (iflow, rflow, nbgauge, a)

Function/Subroutine Documentation

flow1()

subroutine flow1	(	integer, dimension(*)	iflow,
			rflow,
		integer	nbgauge,
			a )

Definition at line 30 of file flow1.F.

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