anim_crk_init.F File Reference

#include "implicit_f.inc"
#include "com04_c.inc"
#include "com_xfem1.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	anim_xfe_init (ixc, ixtg, inod_crk, iel_crk, iadc_crk, iadtg_crk)

Function/Subroutine Documentation

anim_xfe_init()

subroutine anim_xfe_init	(	integer, dimension(nixc,*)	ixc,
		integer, dimension(nixtg,*)	ixtg,
		integer, dimension(*)	inod_crk,
		integer, dimension(*)	iel_crk,
		integer, dimension(4,*)	iadc_crk,
		integer, dimension(3,*)	iadtg_crk )

Definition at line 30 of file anim_crk_init.F.

C-----------------------------------------------
C   M o d u l e s
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   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com04_c.inc"
#include      "com_xfem1.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER IXC(NIXC,*),IXTG(NIXTG,*),INOD_CRK(*),
     .        IEL_CRK(*),IADC_CRK(4,*),IADTG_CRK(3,*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER N,NN,I,II,J,ILEV,NCOUNT,IADC(4)
C=======================================================================
      DO ilev=1,nlevmax
        ncount = 0
c
        DO i=1,numelc
          ii = iel_crk(i)        ! n  systeme xfem
          IF (ii > 0) THEN
            DO j=1,4
              n  = ixc(j+1,i)
              nn = inod_crk(n)   ! n  systeme xfem
              IF (nn > 0) THEN
                ncount  = ncount  + 1
                iadc(j) = iadc_crk(j,ii)   ! adresse sky nodal xfem
                crknod(ilev)%NOD2IAD(ncount) = iadc(j)
              ENDIF
            ENDDO
         ENDIF
        ENDDO
c
        DO i=1,numeltg
          ii = iel_crk(i+numelc)
          IF (ii > 0) THEN
            DO j=1,3
              n  = ixtg(j+1,i)
              nn = inod_crk(n)
              iadc(j) = iadtg_crk(j,ii)
              IF (nn > 0) THEN
                ncount  = ncount  + 1
                crknod(ilev)%NOD2IAD(ncount) = iadc(j)
              ENDIF
            ENDDO
c           ! need one more additional node to reconstruct triangle cut elem
            ncount  = ncount  + 1
            crknod(ilev)%NOD2IAD(ncount) = iadc(3)
          ENDIF
        ENDDO
      ENDDO      !  DO ILEV=1,NLEVMAX
c-----------
      RETURN