s10coork.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	s10coork (x, ixs, ixs10, xx, yy, zz, offg, off, sav, nc, ngl, mxt, ngeo, k11, k12, k13, k14, k15, k16, k17, k18, k19, k10, k22, k23, k24, k25, k26, k27, k28, k29, k20, k33, k34, k35, k36, k37, k38, k39, k30, k44, k45, k46, k47, k48, k49, k40, k55, k56, k57, k58, k59, k50, k66, k67, k68, k69, k60, k77, k78, k79, k70, k88, k89, k80, k99, k90, k00, nel, ismstr)

Function/Subroutine Documentation

s10coork()

subroutine s10coork	(		x,
		integer, dimension(nixs,*)	ixs,
		integer, dimension(6,*)	ixs10,
			xx,
			yy,
			zz,
			offg,
			off,
		double precision, dimension(nel,30)	sav,
		integer, dimension(mvsiz,10)	nc,
		integer, dimension(*)	ngl,
		integer, dimension(*)	mxt,
		integer, dimension(*)	ngeo,
			k11,
			k12,
			k13,
			k14,
			k15,
			k16,
			k17,
			k18,
			k19,
			k10,
			k22,
			k23,
			k24,
			k25,
			k26,
			k27,
			k28,
			k29,
			k20,
			k33,
			k34,
			k35,
			k36,
			k37,
			k38,
			k39,
			k30,
			k44,
			k45,
			k46,
			k47,
			k48,
			k49,
			k40,
			k55,
			k56,
			k57,
			k58,
			k59,
			k50,
			k66,
			k67,
			k68,
			k69,
			k60,
			k77,
			k78,
			k79,
			k70,
			k88,
			k89,
			k80,
			k99,
			k90,
			k00,
		integer	nel,
		integer, intent(in)	ismstr )

Definition at line 30 of file s10coork.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   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      "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
C     REAL
      my_real
     .   x(3,*),
     .  xx(mvsiz,10), yy(mvsiz,10), zz(mvsiz,10),
     .   offg(*), off(*)
      DOUBLE PRECISION 
     .   SAV(NEL,30)
      my_real
     .   k11(9,*)   ,k12(9,*)   ,k13(9,*)    ,k14(9,*)    ,k15(9,*)   ,
     .   k16(9,*)   ,k17(9,*)   ,k18(9,*)    ,k19(9,*)    ,k10(9,*)   ,
     .   k22(9,*)   ,k23(9,*)   ,k24(9,*)    ,k25(9,*)    ,k26(9,*)   ,
     .   k27(9,*)   ,k28(9,*)   ,k29(9,*)    ,k20(9,*)    ,k33(9,*)   ,
     .   k34(9,*)   ,k35(9,*)   ,k36(9,*)    ,k37(9,*)    ,k38(9,*)   ,
     .   k39(9,*)   ,k30(9,*)   ,k44(9,*)   ,k45(9,*)    ,k46(9,*)    ,
     .   k47(9,*)   ,k48(9,*)   ,k49(9,*)   ,k40(9,*)    ,k55(9,*)   ,
     .   k56(9,*)   ,k57(9,*)   ,k58(9,*)   ,k59(9,*)    ,k50(9,*)   ,
     .   k66(9,*)   ,k67(9,*)   ,k68(9,*)   ,k69(9,*)    ,k60(9,*)   ,
     .   k77(9,*)   ,k78(9,*)   ,k79(9,*)   ,k70(9,*)    ,k88(9,*)   ,
     .   k89(9,*)   ,k80(9,*)   ,k99(9,*)   ,k90(9,*)   ,k00(9,*)   
      INTEGER NC(MVSIZ,10),MXT(*), NGL(*),NGEO(*)
      INTEGER IXS(NIXS,*), IXS10(6,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, IPERM1(10),IPERM2(10),N,N1,N2,NN,IUN,J,MXT_1
C     REAL
      DATA iperm1/0,0,0,0,1,2,3,1,2,3/
      DATA iperm2/0,0,0,0,2,3,1,4,4,4/
C-----------------------------------------------
      iun=1
C
      mxt_1 = ixs(1,1)
      DO i=1,nel
        ngeo(i)=ixs(10,i)
        ngl(i)=ixs(11,i)
        mxt(i)=mxt_1
        nc(i,1)=ixs(2,i)
        nc(i,2)=ixs(4,i)
        nc(i,3)=ixs(7,i)
        nc(i,4)=ixs(6,i)
        nc(i,5) =ixs10(1,i)
        nc(i,6) =ixs10(2,i)
        nc(i,7) =ixs10(3,i)
        nc(i,8) =ixs10(4,i)
        nc(i,9) =ixs10(5,i)
        nc(i,10)=ixs10(6,i)
      ENDDO
C----------------------------
C     NODAL COORDINATES     |
C----------------------------
      DO n=1,10
       DO i=1,nel
        nn = max(iun,nc(i,n))
        xx(i,n)=x(1,nn)
        yy(i,n)=x(2,nn)
        zz(i,n)=x(3,nn)
        off(i) = min(one,abs(offg(i)))
       ENDDO
      ENDDO
C-----------
      IF(ismstr==1.OR.(ismstr==2.AND.idtmin(1)==3))THEN
      DO n=1,10
       DO i=1,nel
       IF(abs(offg(i))>one)THEN
        xx(i,n)=sav(i,n)
        yy(i,n)=sav(i,n+10)
        zz(i,n)=sav(i,n+20)
        off(i) = abs(offg(i))-one
       ELSE
        off(i) = offg(i)
       ENDIF
       ENDDO
      ENDDO
C
      ENDIF
C
      DO n=5,10
        n1=iperm1(n)
        n2=iperm2(n)
        DO i=1,nel
          IF(nc(i,n)==0)THEN
            xx(i,n) = half*(xx(i,n1)+xx(i,n2))
            yy(i,n) = half*(yy(i,n1)+yy(i,n2))
            zz(i,n) = half*(zz(i,n1)+zz(i,n2))
          ENDIF
        ENDDO
      ENDDO
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
        k15(j,i)=zero
        k16(j,i)=zero
        k17(j,i)=zero
        k18(j,i)=zero
        k19(j,i)=zero
        k10(j,i)=zero
        k22(j,i)=zero
        k23(j,i)=zero
        k24(j,i)=zero
        k25(j,i)=zero
        k26(j,i)=zero
        k27(j,i)=zero
        k28(j,i)=zero
        k29(j,i)=zero
        k20(j,i)=zero
        k33(j,i)=zero
        k34(j,i)=zero
        k35(j,i)=zero
        k36(j,i)=zero
        k37(j,i)=zero
        k38(j,i)=zero
        k39(j,i)=zero
        k30(j,i)=zero
        k44(j,i)=zero
        k45(j,i)=zero
        k46(j,i)=zero
        k47(j,i)=zero
        k48(j,i)=zero
        k49(j,i)=zero
        k40(j,i)=zero
        k55(j,i)=zero
        k56(j,i)=zero
        k57(j,i)=zero
        k58(j,i)=zero
        k59(j,i)=zero
        k50(j,i)=zero
        k66(j,i)=zero
        k67(j,i)=zero
        k68(j,i)=zero
        k69(j,i)=zero
        k60(j,i)=zero
        k77(j,i)=zero
        k78(j,i)=zero
        k79(j,i)=zero
        k70(j,i)=zero
        k88(j,i)=zero
        k89(j,i)=zero
        k80(j,i)=zero
        k99(j,i)=zero
        k90(j,i)=zero
        k00(j,i)=zero
      ENDDO
      ENDDO
C-----------
      RETURN