spmd_fr_poff.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	spmd_fr_poff (fr_wall, fs, len)

Function/Subroutine Documentation

spmd_fr_poff()

subroutine spmd_fr_poff	(	integer, dimension(*)	fr_wall,
			fs,
		integer	len )

Definition at line 33 of file spmd_fr_poff.F.

C cumul force sur frontiere en parith/Off
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      "com01_c.inc"
#include      "task_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER LEN, FR_WALL(*)
      my_real
     .        fs(len)
C-----------------------------------------------
C   L o c a l  V a r i a b l e s
C-----------------------------------------------
#ifdef MPI
      INTEGER MSGOFF,MSGOFF2,MSGTYP,INFO,P,PMAIN, LOC_PROC, J
      DATA msgoff/110/
      DATA msgoff2/111/
      INTEGER STATUS(MPI_STATUS_SIZE),IERROR
      my_real
     .        ftmp(len)
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
      pmain = fr_wall(nspmd+2)
      loc_proc = ispmd+1
      IF(loc_proc==pmain) THEN
        DO p = 1, nspmd
          IF(loc_proc/=p.AND.fr_wall(p)/=0)THEN
            msgtyp = msgoff
            CALL mpi_recv(ftmp  ,len           ,real  ,it_spmd(p),
     .                    msgtyp,spmd_comm_world,status,ierror    )
C
            DO j = 1, len
              fs(j) = fs(j) + ftmp(j)
            END DO
          ENDIF
        ENDDO
C
        DO p = 1, nspmd
          IF(loc_proc/=p.AND.fr_wall(p)/=0)THEN
            msgtyp = msgoff2
            CALL mpi_send(fs    ,len           ,real  ,it_spmd(p),
     .                    msgtyp,spmd_comm_world,ierror)
          ENDIF
        ENDDO
C
      ELSE
        msgtyp = msgoff
        CALL mpi_send(fs    ,len           ,real  ,it_spmd(pmain),
     .                msgtyp,spmd_comm_world,ierror)
        msgtyp = msgoff2
        CALL mpi_recv(fs    ,len           ,real  ,it_spmd(pmain),
     .                msgtyp,spmd_comm_world,status,ierror    )
      ENDIF
C
#endif
      RETURN