upenritg_last.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	upenritg_last (ixtg, nft, jft, jlt, iad_crktg, iel_crktg, elcutc, nxlay, crkedge)

Function/Subroutine Documentation

upenritg_last()

subroutine upenritg_last	(	integer, dimension(nixtg,*)	ixtg,
		integer	nft,
		integer	jft,
		integer	jlt,
		integer, dimension(3,*)	iad_crktg,
		integer, dimension(*)	iel_crktg,
		integer, dimension(2,*)	elcutc,
		integer	nxlay,
		type (xfem_edge_), dimension(*)	crkedge )

Definition at line 31 of file upenritg_last.F.

C-----------------------------------------------
      USE crackxfem_mod
      use element_mod , only : 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"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER NFT,JFT,JLT,NXLAY
      INTEGER IXTG(NIXTG,*),IEL_CRKTG(*),IAD_CRKTG(3,*),ELCUTC(2,*)
      TYPE (XFEM_EDGE_)   , DIMENSION(*) :: CRKEDGE
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,II,NS,ILEV,IXEL,ILAY,ICUT,LAYCUT,IADC1,IADC2,IADC3,
     .    ELCRK,ELCRKTG
C=======================================================================
        DO i=jft,jlt
        ns   = i+nft
        icut = elcutc(1,ns)                                               
          IF (icut == 0) cycle
        elcutc(2,ns) = icut                                               
        elcrktg = iel_crktg(ns)                                              
        elcrk   = elcrktg + ecrkxfec                                            
        iadc1   = iad_crktg(1,elcrktg)                                      
        iadc2   = iad_crktg(2,elcrktg)                                      
        iadc3   = iad_crktg(3,elcrktg)                                      
          DO ilay = 1,nxlay
          laycut = crkedge(ilay)%LAYCUT(elcrk)                              
          IF (laycut == 0) cycle                                            
            ii = nxel*(ilay-1)    
          DO ixel=1,nxel                                                     
            ilev = ii + ixel                                                 
            crklvset(ilev)%ENR0(2,iadc1) = crklvset(ilev)%ENR0(1,iadc1)  
            crklvset(ilev)%ENR0(2,iadc2) = crklvset(ilev)%ENR0(1,iadc2)  
            crklvset(ilev)%ENR0(2,iadc3) = crklvset(ilev)%ENR0(1,iadc3)  
            ENDDO                 
            ENDDO
          ENDDO
C-----------
      RETURN