i25free_bound.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	i25free_bound (nrtm, mvoisin, irect, stifm, nrtm_free, free_irect_id)

Function/Subroutine Documentation

i25free_bound()

subroutine i25free_bound	(	integer, intent(in)	nrtm,
		integer, dimension(4,nrtm), intent(in)	mvoisin,
		integer, dimension(4,nrtm), intent(in)	irect,
		dimension(nrtm), intent(in)	stifm,
		integer, intent(inout)	nrtm_free,
		integer, dimension(nrtm), intent(inout)	free_irect_id )

Definition at line 31 of file i25free_bound.F.

C Get the list of free IRECT (i.e. at least one edge is at the boundary)
C-----------------------------------------------
C   M o d u l e s
C-----------------------------------------------
#ifdef WITH_ASSERT
      USE debug_mod
#endif
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NRTM
      INTEGER, INTENT(IN) :: MVOISIN(4,NRTM),IRECT(4,NRTM) 
      my_real, INTENT(IN) :: stifm(nrtm)
      INTEGER, INTENT(INOUT) :: NRTM_FREE ! out: number of free irect
      INTEGER, INTENT(INOUT) :: FREE_IRECT_ID(NRTM)! [out] max size = NRTM, list of free irect
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: IRM,IEDG
      LOGICAL :: FREE
 
      nrtm_free= 0
      DO irm=1,nrtm
        IF(stifm(irm)<=zero)cycle
        free = .false.
        DO iedg=1,4
          IF( mvoisin(iedg,irm)==0 )THEN
            IF(.NOT.(irect(3,irm)==irect(4,irm).AND.iedg==3))THEN
              free = .true.
            ENDIF
           ENDIF
        ENDDO
        IF(free) THEN
           nrtm_free = nrtm_free + 1
           free_irect_id(nrtm_free) = irm
        ENDIF
      ENDDO