pmcum3p.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"
#include "mvsiz_p.inc"
#include "parit_c.inc"
#include "vectorize.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	pmcum3p (stir, fsky, fskyv, iadp, 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, nft)

Function/Subroutine Documentation

pmcum3p()

subroutine pmcum3p	(		stir,
			fsky,
			fskyv,
		integer, dimension(2,*)	iadp,
			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,
		integer, intent(in)	nft )

Definition at line 28 of file pmcum3p.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.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      "parit_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: NFT
      INTEGER NC1(MVSIZ),NC2(MVSIZ),NC3(MVSIZ)
      my_real
     .   fskyv(lsky,8),fsky(8,lsky),stir(*),
     .   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)
      INTEGER IADP(2,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II, N,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-------------------------
      IF (ivector == 1) THEN
#include "vectorize.inc"
        DO i=1,nel
          ii = i + nft
          n = iadp(1,ii)
          fskyv(n,4)=-m11(i)
          fskyv(n,5)=-m21(i)
          fskyv(n,6)=-m31(i)
          mx1(i) = m11(i)
          my1(i) = m21(i)
          mz1(i) = m31(i)
          fskyv(n,8)=stir(i)
          n = iadp(2,ii)
          fskyv(n,4)=-m12(i)
          fskyv(n,5)=-m22(i)
          fskyv(n,6)=-m32(i)
          mx2(i) = m12(i)
          my2(i) = m22(i)
          mz2(i) = m32(i)
          fskyv(n,8)=stir(i)
        ENDDO
      ELSE
        DO i=1,nel
          ii = i + nft
          n = iadp(1,ii)
          fsky(4,n)=-m11(i)
          fsky(5,n)=-m21(i)
          fsky(6,n)=-m31(i)
          mx1(i) = m11(i)
          my1(i) = m21(i)
          mz1(i) = m31(i)
          fsky(8,n)=stir(i)
          n = iadp(2,ii)
          fsky(4,n)=-m12(i)
          fsky(5,n)=-m22(i)
          fsky(6,n)=-m32(i)
          mx2(i) = m12(i)
          my2(i) = m22(i)
          mz2(i) = m32(i)
          fsky(8,n)=stir(i)
        ENDDO
      ENDIF
C-----------------------------------------------
      RETURN