flowdec.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	flowdec (iflow)

Function/Subroutine Documentation

flowdec()

subroutine flowdec

(

integer, dimension(*)

iflow

)

Definition at line 30 of file flowdec.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
 
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      "com04_c.inc"
#include      "flowcom.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IFLOW(*)
C-----------------------------------------------
C   E x t e r n a l   F u n c t i o n s
C-----------------------------------------------
      INTEGER  NLOCAL
      EXTERNAL nlocal   
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER IADI, I, ITYP, NINOUT, NNO, NEL, NNN, NBLOC, IADL, PP,
     .        JJ, J, NPROW, NPCOL, K, P, JFORM
      INTEGER PGRID(NSPMD,NSPMD)
C
      iadi=0
      DO i=1,nflow
         ityp =iflow(iadi+2)
         nno  =iflow(iadi+5)
         nel  =iflow(iadi+6)
         nbloc=iflow(iadi+12)
         nprow=iflow(iadi+18)
         npcol=iflow(iadi+19)
         pp=0
         DO j=1,nprow
            DO k=1,npcol
               pp=pp+1
               pgrid(j,k)=pp
            ENDDO
         ENDDO
 
         IF (ityp==1) THEN
C
            ninout=iflow(iadi+4)
            nnn=iflow(iadi+7)
            iadl=iadi+niflow+nno+3*nel+ninout*niioflow+nnn+nel
            pp=1
            jj=0
            DO j=1,nno
               jj=jj+1
               IF (jj==nbloc+1) THEN
                  pp=pp+1
                  IF (pp==nprow+1) pp=1
                  jj=1
               ENDIF
               iflow(iadl+j)=pp  !IBUFR
            ENDDO
            pp=1
            jj=0
            DO j=1,nno
               jj=jj+1
               IF (jj==nbloc+1) THEN
                  pp=pp+1
                  IF (pp==npcol+1) pp=1
                  jj=1
               ENDIF
               iflow(iadl+nno+j)=pp  !IBUFC
            ENDDO
C
            iadl=iadl+4*nno+2*nnn
            pp=1
            jj=0
            DO j=1,nel
               jj=jj+1
               IF (jj==nbloc+1) THEN
                  pp=pp+1
                  IF (pp==nprow+1) pp=1
                  jj=1
               ENDIF
               iflow(iadl+j)=pp   !IBUFELR
            ENDDO
            pp=1
            jj=0
            DO j=1,nel
               jj=jj+1
               IF (jj==nbloc+1) THEN
                  pp=pp+1
                  IF (pp==npcol+1) pp=1
                  jj=1
               ENDIF
               iflow(iadl+nel+j)=pp  !IBUFELC
            ENDDO
C CNP nombre de processeurs possedant chaque noeud
            iadl=iadi+niflow+nno+3*nel+ninout*niioflow+nnn+nel+3*nno
            DO j=1,nno
               iflow(iadl+j)=0
            ENDDO
            DO p=1,nspmd
               DO j=1,nno
                  jj=iflow(iadi+niflow+j)
                  IF (nlocal(jj,p)==1) iflow(iadl+j)=iflow(iadl+j)+1
               ENDDO
            ENDDO
C CNPI nombre de processeurs possedant chaque noeud auxiliaire
            iadl=iadl+nno+nnn
            DO j=1,nnn
               iflow(iadl+j)=0
            ENDDO
            DO p=1,nspmd
               DO j=1,nnn
                  jj=iflow(iadi+niflow+nno+3*nel+ninout*niioflow+j)           
                  IF (nlocal(jj,p)==1) iflow(iadl+j)=iflow(iadl+j)+1
               ENDDO
            ENDDO
C            
         ELSEIF(ityp == 3) THEN
            jform=iflow(iadi+4)
C CNP nombre de processeurs possedant chaque noeud
            iadl = 0 
            IF(jform==1) iadl=iadi+niflow+nno+3*nel+nno
            IF(jform==2) iadl=iadi+niflow+nno+5*nel+nno+nbgauge
            DO j=1,nno
               iflow(iadl+j)=0
            ENDDO
            DO p=1,nspmd
               DO j=1,nno
                  jj=iflow(iadi+niflow+j)
                  IF (nlocal(jj,p)==1) iflow(iadl+j)=iflow(iadl+j)+1
               ENDDO
            ENDDO
C IBUFELR ligne des elements dans la Process Grid
            iadl=iadl+nno+nel
            pp=1
            jj=0
            DO j=1,nel
               jj=jj+1
               IF (jj==nbloc+1) THEN
                  pp=pp+1
                  IF (pp==nprow+1) pp=1
                  jj=1
               ENDIF
               iflow(iadl+j)=pp   !IBUFELR
            ENDDO
C IBUFELC colonne des elements dans la Process Grid
            iadl=iadl+nel
            pp=1
            jj=0
            DO j=1,nel
               jj=jj+1
               IF (jj==nbloc+1) THEN
                  pp=pp+1
                  IF (pp==npcol+1) pp=1
                  jj=1
               ENDIF
               iflow(iadl+j)=pp  !IBUFELC
            ENDDO
         ENDIF
         iadi=iadi+iflow(iadi+14)
      ENDDO
C
      RETURN