masque_artemis.f

Go to the documentation of this file.

!                   *************************
                    SUBROUTINE masque_artemis
!                   *************************
!
!
!***********************************************************************
! ARTEMIS   V6P1                                   21/08/2010
!***********************************************************************
!
!brief    FILLS THE ARRAYS MASK1, MASK2, MASK3, MASK4, MASK5
!+
!+      MASK1: CORRESPONDS TO INCIDENT WAVES (KINC)
!+      MASK2: CORRESPONDS TO FREE EXIT (KSORT)
!+      MASK3: CORRESPONDS TO SOLID BOUNDARY (KLOG)
!+      MASK4: CORRESPONDS TO IMPOSED WAVES (KENT)
!+      MASK5: CORRESPONDS TO INCIDENT POTENTIAL (KPOT)
!
!history  D. AELBRECHT (LNH)
!+        06/07/1999
!+        V5P1
!+
!
!history  C. DENIS (SINETICS)
!+        18/03/2010
!+        V6P0
!+
!
!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
!
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~
!
      USE bief
      USE declarations_telemac
      USE declarations_artemis
!
      USE declarations_special
      IMPLICIT NONE
!
!
      INTEGER IK
!
!-----------------------------------------------------------------------
!
!
!     INITIALISES (SETS TO 0) ALL MASKING VECTORS
!
      CALL os('X=C     ', x=mask1, c=0.d0)
      CALL os('X=C     ', x=mask2, c=0.d0)
      CALL os('X=C     ', x=mask3, c=0.d0)
      CALL os('X=C     ', x=mask4, c=0.d0)
      CALL os('X=C     ', x=mask5, c=0.d0)
!
      DO ik=1,nptfr
!
        IF (lihbor%I(ik).EQ.klog) THEN
          mask3%R(ik) = 1.d0
        ELSEIF (lihbor%I(mesh%KP1BOR%I(ik)).NE.klog) THEN
          IF (lihbor%I(ik).EQ.kinc) THEN
            mask1%R(ik) = 1.d0
          ENDIF
          IF (lihbor%I(ik).EQ.ksort) THEN
            mask2%R(ik) = 1.d0
          ENDIF
          IF (lihbor%I(ik).EQ.kent) THEN
            mask4%R(ik) = 1.d0
          ENDIF
          IF (lihbor%I(ik).EQ.kpot) THEN
            mask5%R(ik) = 1.d0
          ENDIF
        ELSE
          mask3%R(ik) = 1.d0
        ENDIF
      ENDDO ! IK
!
      RETURN
      END