upenric_last.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	upenric_last (ixc, nft, jft, jlt, iadc_crk, iel_crk, elcutc, nxlay, crkedge)

Function/Subroutine Documentation

upenric_last()

subroutine upenric_last	(	integer, dimension(nixc,*)	ixc,
		integer	nft,
		integer	jft,
		integer	jlt,
		integer, dimension(4,*)	iadc_crk,
		integer, dimension(*)	iel_crk,
		integer, dimension(2,*)	elcutc,
		integer	nxlay,
		type (xfem_edge_), dimension(*)	crkedge )

Definition at line 31 of file upenric_last.F.

C-----------------------------------------------
      USE crackxfem_mod
      use element_mod , only : nixc
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 IXC(NIXC,*),IADC_CRK(4,*),IEL_CRK(*),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,ILAY,IXEL,ICUT,LAYCUT,ELCRK,IADC1,IADC2,IADC3,IADC4
C=======================================================================
        DO i=jft,jlt
        ns   = i+nft                                                     
        icut = elcutc(1,ns)                                              
          IF (icut == 0) cycle
        elcutc(2,ns) = icut                                              
        elcrk = iel_crk(ns)                                              
        iadc1 = iadc_crk(1,elcrk)                                        
        iadc2 = iadc_crk(2,elcrk)                                        
        iadc3 = iadc_crk(3,elcrk)                                        
        iadc4 = iadc_crk(4,elcrk)                                        
          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)  
            crklvset(ilev)%ENR0(2,iadc4) = crklvset(ilev)%ENR0(1,iadc4)  
            ENDDO
          ENDDO
        ENDDO
C-------------
      RETURN