global_to_local_point.f

Go to the documentation of this file.

!                   **************************************
                    INTEGER FUNCTION global_to_local_point
!                   **************************************
!
     &(ipoin,mesh)
!
!***********************************************************************
! BIEF   V6P2                                   21/08/2010
!***********************************************************************
!
!brief    In parallel, returns the local value of a point given in
!+    global value. 0 is returned if point not in subdomain. In scalar
!+    mode returns IPOIN. This function replaces array MESH%KNOGL%I
!
!warning   Naive programming with a mere loop, will be optimised
!+    in further versions. Only linear points are considered.
!
!history  J-M HERVOUET (LNHE)
!+        30/07/2012
!+        V6P2
!+
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!| IPOIN          |-->| THE POINT NUMBER IN ORIGINAL SCALAR MESH
!| MESH           |-->| MESH STRUCTURE
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief, ex_global_to_local_point => global_to_local_point
!
      USE declarations_special
      IMPLICIT NONE
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER         , INTENT(IN) :: IPOIN
      TYPE(bief_mesh) , INTENT(IN) :: MESH
!
!+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+-+
!
      INTEGER I
!
!-----------------------------------------------------------------------
!
      IF(ncsize.GT.1) THEN
        global_to_local_point=0
        DO i=1,mesh%NPOIN
          IF(mesh%KNOLG%I(i).EQ.ipoin) THEN
            global_to_local_point=i
          ENDIF
        ENDDO
      ELSE
        global_to_local_point=ipoin
      ENDIF
!
!-----------------------------------------------------------------------
!
      RETURN
      END