c_ncrkxfem.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	c_ncrkxfem (nodglob, inod_crkxfem, inod_l, numnod_l, numnodcrkxfe_l, index, proc, ixc, ixtg, cep_xfe, nodlocal, nodlevxf_l, nodlevxf, nodglobxfe, nod_xfe_l, crkshell)

Function/Subroutine Documentation

c_ncrkxfem()

subroutine c_ncrkxfem	(	integer, dimension(*)	nodglob,
		integer, dimension(*)	inod_crkxfem,
		integer, dimension(*)	inod_l,
		integer	numnod_l,
		integer	numnodcrkxfe_l,
		integer, dimension(*)	index,
		integer	proc,
		integer, dimension(nixc,*)	ixc,
		integer, dimension(nixtg,*)	ixtg,
		integer, dimension(*)	cep_xfe,
		integer, dimension(*)	nodlocal,
		integer, dimension(*)	nodlevxf_l,
		integer, dimension(*)	nodlevxf,
		integer, dimension(*)	nodglobxfe,
		integer	nod_xfe_l,
		type (xfem_shell_), dimension(nlevmax)	crkshell )

Definition at line 29 of file c_ncrkxfem.F.

C-----------------------------------------------
      USE xfem2def_mod
      use element_mod , only : nixc,nixtg
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      "com_xfem1.inc"
#include      "com04_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER  NODGLOB(*),INOD_CRKXFEM(*),INOD_L(*),
     .         NUMNOD_L,NUMNODCRKXFE_L,INDEX(*),PROC,
     .         IXC(NIXC,*),IXTG(NIXTG,*),CEP_XFE(*),
     .         NODLOCAL(*),NODLEVXF_L(*),NODLEVXF(*),
     .         NODGLOBXFE(*),NOD_XFE_L
      TYPE (XFEM_SHELL_) , DIMENSION(NLEVMAX)       :: CRKSHELL
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,NL_L,J,K,NOD,ELTYP,ELEM,
     .        INOD_CRK,NOD_XFE_G,NELCRK
      INTEGER, DIMENSION(:), ALLOCATABLE :: NODTAG
C=======================================================================
!     1d array
      ALLOCATE( nodtag(0:numnod_l+1) )
! ----------------------------------
c     Phantoma nodes each ply
C-----------------------------------------------
      nelcrk = 0
      DO k=1,nlevmax
        DO i=1,crkshell(k)%CRKNUMSHELL
          elem  = crkshell(k)%PHANTOML(i)
          eltyp = crkshell(k)%ELTYPE(i)
          IF (cep_xfe(i) == proc) THEN
            IF(eltyp == 4) THEN
               DO j=1,eltyp
                nod = ixc(j+1,elem)
                IF (nod > 0) THEN
                  IF (inod_crkxfem(nod) > 0) THEN        ! N noeud local xfem
                    nod_xfe_g = crkshell(k)%XNODEG(j,i)  ! ID GLOB NODE PHANTOME SUR PLY
                    nod_xfe_l = nod_xfe_l + 1            ! Local ID Node Phantome Sur Ply
                    nodglobxfe(nod_xfe_l) = nod_xfe_g    ! Id local -> Id global (phant)
                  ENDIF
                ENDIF
              ENDDO
            ELSEIF (eltyp == 3) THEN
              DO j=1,eltyp
               nod = ixtg(j+1,elem)
               IF (nod > 0) THEN
                 IF(inod_crkxfem(nod) > 0)THEN
                   nod_xfe_g = crkshell(k)%XNODEG(j,i)
                   nod_xfe_l = nod_xfe_l + 1
                   nodglobxfe(nod_xfe_l) = nod_xfe_g
                  ENDIF
                ENDIF
              ENDDO
C
c             add one more node (as sh4) for animation files  (3N -> 4N)
C
              nod_xfe_g = crkshell(k)%XNODEG(4,i)
              nod_xfe_l = nod_xfe_l + 1
              nodglobxfe(nod_xfe_l) = nod_xfe_g
            END IF
          ENDIF
        ENDDO
        nelcrk = nelcrk + crkshell(k)%CRKNUMSHELL     ! Nb elements total sur nlevmax
      ENDDO
      numnodxfe = nod_xfe_l                           ! Nb noeuds total sur nlevmax
C
      nodtag(1:numnod_l) = 0
      k = 1 ! the same as K=1,NLEVMAX
      DO i=1,crkshell(k)%CRKNUMSHELL
        eltyp = crkshell(k)%ELTYPE(i)
        elem  = crkshell(k)%PHANTOML(i)
        IF (cep_xfe(i) == proc) THEN
          IF (eltyp == 4) THEN
            DO j=1,eltyp
              nod = ixc(j+1,elem)
              nodtag(nodlocal(nod))=nod
            ENDDO
          ELSEIF (eltyp == 3) THEN
            DO j=1,eltyp
              nod = ixtg(j+1,elem)
              nodtag(nodlocal(nod))=nod
            ENDDO
          END IF
        END IF
      END DO
C---
      nl_l = 0
      DO i=1,numnod_l
        nod = nodtag(i)
        IF (nod > 0) THEN
          IF (inod_crkxfem(nod) > 0) THEN
            nl_l = nl_l + 1
            inod_l(i) = nl_l
            index(nl_l) = inod_crkxfem(nod)
            inod_crk = inod_crkxfem(nod)
            nodlevxf_l(nl_l) = nodlevxf(inod_crk) ! number of copies of a standard xfem node
          ENDIF
        ENDIF
      ENDDO
C---
      numnodcrkxfe_l = nl_l
C---
! ----------------------------------
!     1d array
      DEALLOCATE( nodtag )
! ----------------------------------
      RETURN