surface_deactivation.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	surface_deactivation (ity, nrtm, glocal_surface_id, mseglo, mvoisin)

Function/Subroutine Documentation

surface_deactivation()

subroutine surface_deactivation	(	integer, intent(in)	ity,
		integer, intent(in)	nrtm,
		integer, intent(in)	glocal_surface_id,
		integer, dimension(nrtm), intent(in)	mseglo,
		integer, dimension(4,nrtm), intent(inout)	mvoisin )

Definition at line 28 of file surface_deactivation.F.

!$COMMENT
!       SURFACE_DEACTIVATION description
!           deactivation of neighbour surface
!       SURFACE_DEACTIVATION organization
!           - for a given K deactivated surface :
!              loop over the surfaces to deactivate the neighbour K
!$ENDCOMMENT
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
 
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
        INTEGER, INTENT(in) :: ITY ! kind of interface
        INTEGER, INTENT(in) :: NRTM ! number of surface
        INTEGER, INTENT(in) :: GLOCAL_SURFACE_ID ! global id of the surface
        INTEGER, DIMENSION(NRTM), INTENT(in) :: MSEGLO ! global id of the k surface 
        INTEGER, DIMENSION(4,NRTM), INTENT(inout) :: MVOISIN ! neighbouring array
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
        INTEGER :: J
C-----------------------------------------------
        ! ----------------------
        DO j=1,nrtm
            ! --------------
            ! deactivation of neighbouring 
            IF(mvoisin(1,j)==glocal_surface_id) mvoisin(1,j) = 0
            IF(mvoisin(2,j)==glocal_surface_id) mvoisin(2,j) = 0
            IF(mvoisin(3,j)==glocal_surface_id) mvoisin(3,j) = 0
            IF(mvoisin(4,j)==glocal_surface_id) mvoisin(4,j) = 0
            ! --------------
                
            ! --------------
            ! deactivation of the surface SURFACE_ID
            IF(ity==24.AND.mseglo(j)==glocal_surface_id) THEN
                mvoisin(1,j) = 0
                mvoisin(2,j) = 0
                mvoisin(3,j) = 0
                mvoisin(4,j) = 0
            ENDIF
            ! --------------
        ENDDO
        ! ----------------------
        RETURN