ccoor3_crk.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com08_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	ccoor3_crk (jft, jlt, nft, ixc, ngl, pid, mat, ilev, iel_crk, iadc_crk, offg, off, sigy, thke, thk0, dt1c, vl1, vl2, vl3, vl4, vrl1, vrl2, vrl3, vrl4, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4)

Function/Subroutine Documentation

ccoor3_crk()

subroutine ccoor3_crk	(	integer	jft,
		integer	jlt,
		integer	nft,
		integer, dimension(nixc,*)	ixc,
		integer, dimension(mvsiz)	ngl,
		integer, dimension(mvsiz)	pid,
		integer, dimension(mvsiz)	mat,
		integer	ilev,
		integer, dimension(*)	iel_crk,
		integer, dimension(4,*)	iadc_crk,
			offg,
			off,
			sigy,
			thke,
			thk0,
			dt1c,
			vl1,
			vl2,
			vl3,
			vl4,
			vrl1,
			vrl2,
			vrl3,
			vrl4,
			x1,
			x2,
			x3,
			x4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4 )

Definition at line 31 of file ccoor3_crk.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   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      "com08_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT,JLT,ILEV
      INTEGER NGL(MVSIZ),IXC(NIXC,*),MAT(MVSIZ),PID(MVSIZ),
     .   IEL_CRK(*),IADC_CRK(4,*),NFT
C     REAL
      my_real
     .   offg(*),off(*),sigy(*),
     .   vl1(mvsiz,3),vl2(mvsiz,3),vl3(mvsiz,3),vl4(mvsiz,3),
     .   vrl1(mvsiz,3),vrl2(mvsiz,3),vrl3(mvsiz,3),vrl4(mvsiz,3),dt1c(*),thke(*),
     .   x1(mvsiz),x2(mvsiz),x3(mvsiz),x4(mvsiz),y1(mvsiz),
     .   y2(mvsiz),y3(mvsiz),y4(mvsiz),z1(mvsiz),z2(mvsiz),
     .   z3(mvsiz),z4(mvsiz),thk0(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,ELCRK,IADC1,IADC2,IADC3,IADC4
C     REAL
      my_real
     .   off_l
C-----------------------------------------------
C get elem from activ level: (superposed element)
C-----------------------------------------------
        DO i=jft,jlt
          elcrk = iel_crk(i+nft)
          iadc1 = iadc_crk(1,elcrk)
          iadc2 = iadc_crk(2,elcrk)
          iadc3 = iadc_crk(3,elcrk)
          iadc4 = iadc_crk(4,elcrk)
C--------------
C  COORDINATES:
C--------------
C  node 1:
          x1(i) = crkavx(ilev)%X(1,iadc1)
          y1(i) = crkavx(ilev)%X(2,iadc1)
          z1(i) = crkavx(ilev)%X(3,iadc1)
C  node 2:
          x2(i) = crkavx(ilev)%X(1,iadc2)
          y2(i) = crkavx(ilev)%X(2,iadc2)
          z2(i) = crkavx(ilev)%X(3,iadc2)
C  node 3:
          x3(i) = crkavx(ilev)%X(1,iadc3)
          y3(i) = crkavx(ilev)%X(2,iadc3)
          z3(i) = crkavx(ilev)%X(3,iadc3)
C  node 4:
          x4(i) = crkavx(ilev)%X(1,iadc4)
          y4(i) = crkavx(ilev)%X(2,iadc4)
          z4(i) = crkavx(ilev)%X(3,iadc4)
C--------------
C  VELOCITIES :
C--------------
C  node 1:
          vl1(i,1)  = crkavx(ilev)%V(1,iadc1)
          vl1(i,2)  = crkavx(ilev)%V(2,iadc1)
          vl1(i,3)  = crkavx(ilev)%V(3,iadc1)
          vrl1(i,1) = crkavx(ilev)%VR(1,iadc1)
          vrl1(i,2) = crkavx(ilev)%VR(2,iadc1)
          vrl1(i,3) = crkavx(ilev)%VR(3,iadc1)
C  node 2:
          vl2(i,1)  = crkavx(ilev)%V(1,iadc2)
          vl2(i,2)  = crkavx(ilev)%V(2,iadc2)
          vl2(i,3)  = crkavx(ilev)%V(3,iadc2)
          vrl2(i,1) = crkavx(ilev)%VR(1,iadc2)
          vrl2(i,2) = crkavx(ilev)%VR(2,iadc2)
          vrl2(i,3) = crkavx(ilev)%VR(3,iadc2)
C  node 3:
          vl3(i,1)  = crkavx(ilev)%V(1,iadc3)
          vl3(i,2)  = crkavx(ilev)%V(2,iadc3)
          vl3(i,3)  = crkavx(ilev)%V(3,iadc3)
          vrl3(i,1) = crkavx(ilev)%VR(1,iadc3)
          vrl3(i,2) = crkavx(ilev)%VR(2,iadc3)
          vrl3(i,3) = crkavx(ilev)%VR(3,iadc3)
C  node 4:
          vl4(i,1)  = crkavx(ilev)%V(1,iadc4)
          vl4(i,2)  = crkavx(ilev)%V(2,iadc4)
          vl4(i,3)  = crkavx(ilev)%V(3,iadc4)
          vrl4(i,1) = crkavx(ilev)%VR(1,iadc4)
          vrl4(i,2) = crkavx(ilev)%VR(2,iadc4)
          vrl4(i,3) = crkavx(ilev)%VR(3,iadc4)
        END DO
C-----------------------------------------------
      DO i=jft,jlt
        mat(i)=ixc(1,i)
        ngl(i)=ixc(7,i)
        pid(i)=ixc(6,i)
        thk0(i) = thke(i)
      END DO
C----------------------------------
      off_l = zero
      DO i=jft,jlt
        dt1c(i) = dt1
        off(i) = min(one,abs(offg(i)))
        off_l  = min(off_l,offg(i))
        sigy(i) = 1.e30
      ENDDO
C----------------------------------
      IF(off_l<zero)THEN
        DO i=jft,jlt
         IF(offg(i)<zero)THEN
           vl1(i,1)=zero
           vl1(i,2)=zero
           vl1(i,3)=zero
           vl2(i,1)=zero
           vl2(i,2)=zero
           vl2(i,3)=zero
           vl3(i,1)=zero
           vl3(i,2)=zero
           vl3(i,3)=zero
           vl4(i,1)=zero
           vl4(i,2)=zero
           vl4(i,3)=zero
           vrl1(i,1)=zero
           vrl1(i,2)=zero
           vrl1(i,3)=zero
           vrl2(i,1)=zero
           vrl2(i,2)=zero
           vrl2(i,3)=zero
           vrl3(i,1)=zero
           vrl3(i,2)=zero
           vrl3(i,3)=zero
           vrl4(i,1)=zero
           vrl4(i,2)=zero
           vrl4(i,3)=zero
         ENDIF
        ENDDO
      ENDIF
C----------------------------------
      RETURN