spmd_anim_ply_init.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	spmd_anim_ply_init ()

Function/Subroutine Documentation

spmd_anim_ply_init()

subroutine spmd_anim_ply_init

Definition at line 32 of file spmd_anim_ply_init.F.

C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
      USE plyxfem_mod
      USE spmd_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#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      "com01_c.inc"
#include      "param_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
#ifdef MPI
      INTEGER, DIMENSION(:),ALLOCATABLE :: PLYNODTAG,SENDBUF,RECBUF
      INTEGER MAXSENDSZ,RECSZ,PLYSZ
      INTEGER I,J,P,ND, ITAG, IDCMAXP(NSPMD-1), REQ(NSPMD-1),NODCOUNT
      INTEGER MSGOFF,MSGOFF2,STAT(MPI_STATUS_SIZE,NSPMD-1), IERR
C-----------------------------------------------
      DATA msgoff/7054/
C-----------------------------------------------
C
C Allocate Global PLY Structure
C
      IF (ispmd==0)THEN
        ALLOCATE (plyiadnod(nplymax,nspmd))
        plyiadnod = 0
      ELSE
         ALLOCATE (plyiadnod(1,1))
      ENDIF
      plyiadnod = 0
 
 
C Need to rebuild a global PLYNOD Structure on processor 0 to have
C a conversion array from Global Node ID to Internal PLYNOD ID
 
      DO i=1,nplymax
        IF (ispmd==0)THEN
C PROC 0, RECEIVE NODES NUM for each ply
           DO p=2,nspmd
              itag=msgoff
              CALL mpi_recv(recsz,1, mpi_integer, it_spmd(p),
     .                      itag, spmd_comm_world, stat, ierr)
              plyiadnod(i,p) = recsz
           ENDDO
        ELSE
           itag=msgoff
           CALL mpi_send(plynod(i)%PLYNUMNODS,1,mpi_integer,it_spmd(1),
     .                   itag,spmd_comm_world,ierr)
        ENDIF
      ENDDO
 
C
#endif
      RETURN