spmd_fvb_aoff.F File Reference

#include "implicit_f.inc"
#include "spmd.inc"
#include "task_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	spmd_fvb_aoff (fvel2fa)

Function/Subroutine Documentation

spmd_fvb_aoff()

subroutine spmd_fvb_aoff

(

integer, dimension(*)

fvel2fa

)

Definition at line 33 of file spmd_fvb_aoff.F.

C-----------------------------------------------
C   M o d u l e s 
C-----------------------------------------------
      USE fvbag_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
        USE spmd_comm_world_mod, ONLY : spmd_comm_world
#include      "implicit_f.inc"
C-----------------------------------------------------------------
C   M e s s a g e   P a s s i n g
C-----------------------------------------------
#include "spmd.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER FVEL2FA(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
#ifdef MPI
      INTEGER ELOFF, I, PMAIN, J, K, KK, N, NN, N1, N2, N3, NND,
     .        ITAG, NNTR, IERR, STAT(MPI_STATUS_SIZE), JJ, IAD,
     .        LEN, MSGOFF,MSGOFF2
      INTEGER, DIMENSION(:), ALLOCATABLE :: OFFTR, IBUF, ITAGT
C-----------------------------------------------
      DATA msgoff/7041/
      DATA msgoff2/7042/
C-----------------------------------------------
C
      eloff=0
      DO i=1,nfvbag
         pmain=fvspmd(i)%PMAIN
         IF (ispmd==0) THEN
            IF (ispmd==pmain-1) THEN
               ALLOCATE(offtr(fvdata(i)%NNTR))
               DO j=1,fvdata(i)%NNTR
                  offtr(j)=0
               ENDDO
               DO j=1,fvdata(i)%NPOLH
                  DO k=fvdata(i)%IFVPADR(j),fvdata(i)%IFVPADR(j+1)-1
                     kk=fvdata(i)%IFVPOLH(k)
                     DO n=fvdata(i)%IFVTADR(kk),
     .                         fvdata(i)%IFVTADR(kk+1)-1
                        nn=fvdata(i)%IFVPOLY(n)
                        IF (nn>0) THEN
                           n1=fvdata(i)%IFVTRI_ANIM(1,nn)
                           n2=fvdata(i)%IFVTRI_ANIM(2,nn)
                           n3=fvdata(i)%IFVTRI_ANIM(3,nn)
                           nnd=1
                           IF (n2/=n1) nnd=nnd+1
                           IF (n3/=n2.AND.n3/=n1) nnd=nnd+1
C
                           nn=fvel2fa(eloff+nn)-eloff
                           IF (nnd==3) offtr(nn)=1
                        ENDIF
                     ENDDO
                  ENDDO
               ENDDO
               eloff=eloff+fvdata(i)%NNTR
               CALL write_c_c(offtr,fvdata(i)%NNTR)
               DEALLOCATE(offtr)
            ELSE
               itag=msgoff
               CALL mpi_recv(nntr, 1, mpi_integer, it_spmd(pmain),
     .                       itag, spmd_comm_world, stat, ierr)
C
               len=2*nntr
               ALLOCATE(ibuf(len), offtr(nntr))
               itag=msgoff2
               CALL mpi_recv(ibuf, len, mpi_integer, it_spmd(pmain),
     .                       itag, spmd_comm_world, stat, ierr)
C
               DO j=1,nntr
                  offtr(j)=0
               ENDDO
               jj=0
               DO j=1,nntr
                  nn=ibuf(jj+1)
                  IF (nn /=0 ) THEN
                    nn=fvel2fa(eloff+nn)-eloff
                    offtr(nn)=ibuf(jj+2)
                  ENDIF
                  jj=jj+2
               ENDDO
               eloff=eloff+nntr
               CALL write_c_c(offtr,nntr)
               DEALLOCATE(ibuf, offtr)
            ENDIF
         ELSE
            IF (ispmd==pmain-1) THEN
               itag=msgoff
               CALL mpi_send(fvdata(i)%NNTR, 1, mpi_integer,
     .                       it_spmd(1), itag, spmd_comm_world,
     .                       ierr)
C
               len=2*fvdata(i)%NNTR
               ALLOCATE(ibuf(len), itagt(fvdata(i)%NNTR))
               DO j=1,len
                  ibuf(j)=0
               ENDDO
               DO j=1,fvdata(i)%NNTR
                  itagt(j)=0
               ENDDO
               iad=0
               DO j=1,fvdata(i)%NPOLH
                  DO k=fvdata(i)%IFVPADR(j),fvdata(i)%IFVPADR(j+1)-1
                     kk=fvdata(i)%IFVPOLH(k)
                     DO n=fvdata(i)%IFVTADR(kk),
     .                         fvdata(i)%IFVTADR(kk+1)-1
                        nn=fvdata(i)%IFVPOLY(n)
                        IF (nn>0.AND.itagt(nn)==0) THEN
                           n1=fvdata(i)%IFVTRI_ANIM(1,nn)
                           n2=fvdata(i)%IFVTRI_ANIM(2,nn)
                           n3=fvdata(i)%IFVTRI_ANIM(3,nn)
                           nnd=1
                           IF (n2/=n1) nnd=nnd+1
                           IF (n3/=n2.AND.n3/=n1) nnd=nnd+1
C
                           ibuf(iad+1)=nn
                           IF (nnd==3) ibuf(iad+2)=1
                           iad=iad+2
                           itagt(nn)=1
                        ENDIF
                     ENDDO
                  ENDDO
               ENDDO
               itag=msgoff2
               CALL mpi_send(ibuf, len, mpi_integer, it_spmd(1),
     .                       itag, spmd_comm_world,  ierr)
C
               DEALLOCATE(ibuf, itagt)              
            ENDIF
         ENDIF
      ENDDO
C
#endif
      RETURN