parcom2.f

Go to the documentation of this file.

!                   ******************
                    SUBROUTINE parcom2
!                   ******************
!
     &( x1 , x2 , x3 , npoin , nplan , icom , ian , mesh )
!
!***********************************************************************
! BIEF   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    COMPLEMENTS A VECTOR AT THE INTERFACES BETWEEN
!+                SUB-DOMAINS.
!+
!+            X CAN BE A BLOCK OF VECTORS. IN THIS CASE, ALL THE
!+                VECTORS IN THE BLOCK ARE TREATED.
!
!note     IMPORTANT : FROM RELEASE 5.9 ON, IDENTICAL VALUES ARE
!+                     ENSURED AT INTERFACE POINTS SO THAT DIFFERENT
!+                     PROCESSORS WILL ALWAYS MAKE THE SAME DECISION
!+                     IN TESTS ON REAL NUMBERS.
!
!warning  IF THE VECTORS HAVE A SECOND DIMENSION, IT IS
!+            IGNORED FOR THE TIME BEING
!
!history  J-M HERVOUET (LNHE)
!+        13/08/08
!+        V5P9
!+   AFTER REINHARD HINKELMANN (HANNOVER UNI.)
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        13/07/2010
!+        V6P0
!+   Translation of French comments within the FORTRAN sources into
!+   English comments
!
!history  N.DURAND (HRW), S.E.BOURBAN (HRW)
!+        21/08/2010
!+        V6P0
!+   Creation of DOXYGEN tags for automated documentation and
!+   cross-referencing of the FORTRAN sources
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IAN            |-->| NUMBER OF VECTORS TO BE CONDIDERED (1, 2 OR 3)
!| ICOM           |-->| COMMUNICATION MODE
!|                |   | = 1 : VALUE WITH MAXIMUM ABSOLUTE VALUE TAKEN
!|                |   | = 2 : CONTRIBUTIONS ADDED
!|                |   | = 3 : MAXIMUM CONTRIBUTION RETAINED
!|                |   | = 4 : MINIMUM CONTRIBUTION RETAINED
!| MESH           |-->| MESH STRUCTURE
!| NPLAN          |-->| NUMBER OF PLANES IN THE 3D MESH OF PRISMS
!| NPOIN          |-->| NUMBER OF POINTS
!| X1             |<->| VECTOR TO BE COMPLETED
!| X2             |<->| VECTOR TO BE COMPLETED
!| X3             |<->| VECTOR TO BE COMPLETED
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_parcom2 => parcom2
!
      USE declarations_special
      IMPLICIT NONE
!
!-----------------------------------------------------------------------
!
      INTEGER, INTENT(IN) :: ICOM,NPOIN,NPLAN,IAN
!
!     STRUCTURES: VECTORS OR BLOCKS
!
      TYPE(bief_mesh) , INTENT(INOUT) :: MESH
      DOUBLE PRECISION, INTENT(INOUT) :: X1(npoin,nplan)
      DOUBLE PRECISION, INTENT(INOUT) :: X2(npoin,nplan)
      DOUBLE PRECISION, INTENT(INOUT) :: X3(npoin,nplan)
!
!-----------------------------------------------------------------------
!
      CALL paraco(x1,x2,x3,npoin,icom,ian,nplan,
     &            mesh%NB_NEIGHB,mesh%NB_NEIGHB_PT%I,mesh%LIST_SEND%I,
     &            mesh%NH_COM%I,mesh%NH_COM%DIM1,mesh%BUF_SEND%R,
     &            mesh%BUF_RECV%R,mesh%BUF_SEND%DIM1)
!
!     RELEASE 5.9 : ENSURES SAME VALUES AT INTERFACE POINTS
!                   SHARED BY SEVERAL SUB-DOMAINS
!
      IF(icom.EQ.2.AND.ncsize.GT.2) THEN
!
      CALL paraco(x1,x2,x3,npoin,1,ian,nplan,
     &            mesh%NB_NEIGHB,mesh%NB_NEIGHB_PT%I,mesh%LIST_SEND%I,
     &            mesh%NH_COM%I,mesh%NH_COM%DIM1,mesh%BUF_SEND%R,
     &            mesh%BUF_RECV%R,mesh%BUF_SEND%DIM1)
!
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END