pmcum3.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	pmcum3 (m, stir, stifr, mx1, mx2, my1, my2, mz1, mz2, nc1, nc2, nc3, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z, m11, m12, m21, m22, m31, m32, nel)

Function/Subroutine Documentation

pmcum3()

subroutine pmcum3	(		m,
			stir,
			stifr,
			mx1,
			mx2,
			my1,
			my2,
			mz1,
			mz2,
		integer, dimension(mvsiz)	nc1,
		integer, dimension(mvsiz)	nc2,
		integer, dimension(mvsiz)	nc3,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z,
			m11,
			m12,
			m21,
			m22,
			m31,
			m32,
		integer, intent(in)	nel )

Definition at line 28 of file pmcum3.F.

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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER NC1(MVSIZ),NC2(MVSIZ),NC3(MVSIZ)
      my_real
     .   m(3,*),stir(*),stifr(*),
     .   mx1(mvsiz),my1(mvsiz),mz1(mvsiz),
     .   mx2(mvsiz),my2(mvsiz),mz2(mvsiz),
     .   e1x(mvsiz),e1y(mvsiz),e1z(mvsiz),
     .   e2x(mvsiz),e2y(mvsiz),e2z(mvsiz),
     .   e3x(mvsiz),e3y(mvsiz),e3z(mvsiz),
     .   m11(mvsiz),m12(mvsiz),m21(mvsiz),
     .   m22(mvsiz),m31(mvsiz),m32(mvsiz)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
      my_real
     .   g11(mvsiz),g21(mvsiz),g31(mvsiz),
     .   g12(mvsiz),g22(mvsiz),g32(mvsiz)
C-----------------------------------------------
      DO i=1,nel
        g11(i)=m11(i)
        g21(i)=m21(i)
        g31(i)=m31(i)
        g12(i)=m12(i)
        g22(i)=m22(i)
        g32(i)=m32(i)
      ENDDO
C
      DO i=1,nel
        m11(i)=e1x(i)*g11(i)+e2x(i)*g21(i)+e3x(i)*g31(i)
        m21(i)=e1y(i)*g11(i)+e2y(i)*g21(i)+e3y(i)*g31(i)
        m31(i)=e1z(i)*g11(i)+e2z(i)*g21(i)+e3z(i)*g31(i)
        m12(i)=e1x(i)*g12(i)+e2x(i)*g22(i)+e3x(i)*g32(i)
        m22(i)=e1y(i)*g12(i)+e2y(i)*g22(i)+e3y(i)*g32(i)
        m32(i)=e1z(i)*g12(i)+e2z(i)*g22(i)+e3z(i)*g32(i)
      ENDDO
C-------------------------
C     ASSEMBLAGE
C-------------------------
      DO i=1,nel
        m(1,nc1(i))=m(1,nc1(i))-m11(i)
        m(2,nc1(i))=m(2,nc1(i))-m21(i)
        m(3,nc1(i))=m(3,nc1(i))-m31(i)
        mx1(i) = m11(i)
        my1(i) = m21(i)
        mz1(i) = m31(i)
        stifr(nc1(i))=stifr(nc1(i))+stir(i)
C
        m(1,nc2(i))=m(1,nc2(i))-m12(i)
        m(2,nc2(i))=m(2,nc2(i))-m22(i)
        m(3,nc2(i))=m(3,nc2(i))-m32(i)
        mx2(i) = m12(i)
        my2(i) = m22(i)
        mz2(i) = m32(i)
        stifr(nc2(i))=stifr(nc2(i))+stir(i)
      ENDDO
C-----------------------------------------------
      RETURN