activ_xfem.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com_xfem1.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	activ_xfem (iparg, nft, jft, jlt, nxlay, ng, elcutc, iel_crk, ity, crkedge)

Function/Subroutine Documentation

activ_xfem()

subroutine activ_xfem	(	integer, dimension(nparg,*)	iparg,
		integer	nft,
		integer	jft,
		integer	jlt,
		integer	nxlay,
		integer	ng,
		integer, dimension(2,*)	elcutc,
		integer, dimension(*)	iel_crk,
		integer	ity,
		type (xfem_edge_), dimension(*)	crkedge )

Definition at line 30 of file activ_xfem.F.

C-----------------------------------------------
      USE crackxfem_mod
C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com_xfem1.inc"
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NFT,JFT,JLT,NXLAY,IPARG(NPARG,*),NG,
     .        IEL_CRK(*),ELCUTC(2,*),ITY
      TYPE (XFEM_EDGE_)   , DIMENSION(*) :: CRKEDGE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,ELCRK,IELCUT,ILAYCUT,ILAY,IGON
      INTEGER ISHEON(MVSIZ)
C=======================================================================
      DO i=jft,jlt
        ielcut = elcutc(1,i+nft)
        isheon(i) = 0
        IF (ielcut /= 0) THEN
          elcrk = iel_crk(i+nft)
          IF (ity == 7) elcrk = elcrk + ecrkxfec
          DO ilay = 1,nxlay
            ilaycut = crkedge(ilay)%LAYCUT(elcrk)
            IF (abs(ilaycut) == 1) isheon(i) = 1
          ENDDO
        ENDIF
      ENDDO
C
      igon = 0
      DO i=jft,jlt
        IF (isheon(i) > 0) THEN
          igon = 1
        ENDIF
      ENDDO
C----------------------------------------
C     TEST FOR NEW GROUP ACTIVATION
C----------------------------------------
      IF (igon == 1) iparg(70,ng) = igon
C-----------
      RETURN