s4coork.F File Reference

#include "implicit_f.inc"
#include "scr18_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s4coork (x, ixs, x1, x2, x3, x4, y1, y2, y3, y4, z1, z2, z3, z4, offg, off, sav, nc1, nc2, nc3, nc4, ngl, mxt, ngeo, k11, k12, k13, k14, k22, k23, k24, k33, k34, k44, nel, ismstr)

Function/Subroutine Documentation

s4coork()

subroutine s4coork	(		x,
		integer, dimension(nixs,*)	ixs,
			x1,
			x2,
			x3,
			x4,
			y1,
			y2,
			y3,
			y4,
			z1,
			z2,
			z3,
			z4,
			offg,
			off,
		double precision, dimension(nel,9)	sav,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
		integer, dimension(*)	ngl,
		integer, dimension(*)	mxt,
		integer, dimension(*)	ngeo,
			k11,
			k12,
			k13,
			k14,
			k22,
			k23,
			k24,
			k33,
			k34,
			k44,
		integer	nel,
		integer, intent(in)	ismstr )

Definition at line 30 of file s4coork.F.

      use element_mod , only : nixs
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      "scr18_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: ISMSTR
      INTEGER NEL
      my_real
     .   x(3,*),
     .   x1(*), x2(*), x3(*), x4(*),  
     .   y1(*), y2(*), y3(*), y4(*),  
     .   z1(*), z2(*), z3(*), z4(*), 
     .   offg(*), off(*)
      double precision
     .  sav(nel,9)
      my_real
     .   k11(9,*)   ,k12(9,*)   ,k13(9,*)  ,k14(9,*)    ,k22(9,*)  ,
     .   k23(9,*)   ,k24(9,*)   ,k33(9,*)   ,k34(9,*)   ,k44(9,*)     
      INTEGER NC1(*), NC2(*), NC3(*), NC4(*),MXT(*), NGL(*),NGEO(*)
      INTEGER IXS(NIXS,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J,MXT_1
C     REAL
C-----------------------------------------------
C
      mxt_1 = ixs(1,1)
      DO i=1,nel
      ngeo(i)=ixs(10,i)
      ngl(i)=ixs(11,i)
      mxt(i)=mxt_1
      nc1(i)=ixs(2,i)
      nc2(i)=ixs(4,i)
      nc3(i)=ixs(7,i)
      nc4(i)=ixs(6,i)
      ENDDO
C----------------------------
C     NODAL COORDINATES     |
C----------------------------
      DO i=1,nel
        x1(i)=x(1,nc1(i))
        y1(i)=x(2,nc1(i))
        z1(i)=x(3,nc1(i))
        x2(i)=x(1,nc2(i))
        y2(i)=x(2,nc2(i))
        z2(i)=x(3,nc2(i))
        x3(i)=x(1,nc3(i))
        y3(i)=x(2,nc3(i))
        z3(i)=x(3,nc3(i))
        x4(i)=x(1,nc4(i))
        y4(i)=x(2,nc4(i))
        z4(i)=x(3,nc4(i))
        off(i) = min(one,abs(offg(i)))
      ENDDO
C-----------
      IF(ismstr==1.OR.(ismstr==2.AND.idtmin(1)==3))THEN
       DO i=1,nel
       IF(abs(offg(i))>one)THEN
        x1(i)=sav(i,1)
        y1(i)=sav(i,2)
        z1(i)=sav(i,3)
        x2(i)=sav(i,4)
        y2(i)=sav(i,5)
        z2(i)=sav(i,6)
        x3(i)=sav(i,7)
        y3(i)=sav(i,8)
        z3(i)=sav(i,9)
        x4(i)=zero
        y4(i)=zero
        z4(i)=zero
        off(i) = abs(offg(i))-one
       ELSE
        off(i) = offg(i)
       ENDIF
       ENDDO
C
      ENDIF
C-----------
      DO j=1,9      
      DO i=1,nel
        k11(j,i)=zero
        k12(j,i)=zero
        k13(j,i)=zero
        k14(j,i)=zero
        k22(j,i)=zero
        k23(j,i)=zero
        k24(j,i)=zero
        k33(j,i)=zero
        k34(j,i)=zero
        k44(j,i)=zero
      ENDDO
      ENDDO
C-----------
      RETURN