s8ederig3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s8ederig3 (px1, px2, px3, px4, px5, px6, px7, px8, py1, py2, py3, py4, py5, py6, py7, py8, pz1, pz2, pz3, pz4, pz5, pz6, pz7, pz8, aji1, aji2, aji3, aji4, aji5, aji6, aji7, aji8, aji9, pr, ps, pt, nel)

Function/Subroutine Documentation

s8ederig3()

subroutine s8ederig3	(		px1,
			px2,
			px3,
			px4,
			px5,
			px6,
			px7,
			px8,
			py1,
			py2,
			py3,
			py4,
			py5,
			py6,
			py7,
			py8,
			pz1,
			pz2,
			pz3,
			pz4,
			pz5,
			pz6,
			pz7,
			pz8,
			aji1,
			aji2,
			aji3,
			aji4,
			aji5,
			aji6,
			aji7,
			aji8,
			aji9,
			pr,
			ps,
			pt,
		integer, intent(in)	nel )

Definition at line 29 of file s8ederig3.F.

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-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      my_real
     .   px1(*), px2(*), px3(*), px4(*),
     .   px5(*), px6(*), px7(*), px8(*),
     .   py1(*), py2(*), py3(*), py4(*),
     .   py5(*), py6(*), py7(*), py8(*),
     .   pz1(*), pz2(*), pz3(*), pz4(*), 
     .   pz5(*), pz6(*), pz7(*), pz8(*),
     .   aji1(*), aji2(*), aji3(*),
     .   aji4(*), aji5(*), aji6(*),
     .   aji7(*), aji8(*), aji9(*),
     .   pr(8),ps(8),pt(8)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J 
C     REAL
C                                                                     12
      my_real
     .   a1pr1,a1pr3,a1pr5,a1pr7,
     .   a4pr1,a4pr3,a4pr5,a4pr7,
     .   a7pr1,a7pr3,a7pr5,a7pr7,
     .   a2ps1,a2ps2,a2ps5,a2ps6,
     .   a5ps1,a5ps2,a5ps5,a5ps6,
     .   a8ps1,a8ps2,a8ps5,a8ps6,
     .   a3pt1,a3pt2,a3pt3,a3pt4,
     .   a6pt1,a6pt2,a6pt3,a6pt4,
     .   a9pt1,a9pt2,a9pt3,a9pt4
C-----------------------------------------------
cc-----pxI
      DO i=1,nel
        a1pr1=aji1(i)*pr(1)
        a1pr3=aji1(i)*pr(3)
        a1pr5=aji1(i)*pr(5)
        a1pr7=aji1(i)*pr(7)
        a2ps1=aji2(i)*ps(1)
        a2ps2=aji2(i)*ps(2)
        a2ps5=aji2(i)*ps(5)
        a2ps6=aji2(i)*ps(6)
        a3pt1=aji3(i)*pt(1)
        a3pt2=aji3(i)*pt(2)
        a3pt3=aji3(i)*pt(3)
        a3pt4=aji3(i)*pt(4)
       px1(i)= a1pr1+a2ps1+a3pt1
       px2(i)=-a1pr1+a2ps2+a3pt2
       px3(i)= a1pr3-a2ps2+a3pt3
       px4(i)=-a1pr3-a2ps1+a3pt4
       px5(i)= a1pr5+a2ps5-a3pt1
       px6(i)=-a1pr5+a2ps6-a3pt2
       px7(i)= a1pr7-a2ps6-a3pt3
       px8(i)=-a1pr7-a2ps5-a3pt4
      ENDDO
C
      DO i=1,nel
C
        a4pr1=aji4(i)*pr(1)
        a4pr3=aji4(i)*pr(3)
        a4pr5=aji4(i)*pr(5)
        a4pr7=aji4(i)*pr(7)
        a5ps1=aji5(i)*ps(1)
        a5ps2=aji5(i)*ps(2)
        a5ps5=aji5(i)*ps(5)
        a5ps6=aji5(i)*ps(6)
        a6pt1=aji6(i)*pt(1)
        a6pt2=aji6(i)*pt(2)
        a6pt3=aji6(i)*pt(3)
        a6pt4=aji6(i)*pt(4)
       py1(i)= a4pr1+a5ps1+a6pt1
       py2(i)=-a4pr1+a5ps2+a6pt2
       py3(i)= a4pr3-a5ps2+a6pt3
       py4(i)=-a4pr3-a5ps1+a6pt4
       py5(i)= a4pr5+a5ps5-a6pt1
       py6(i)=-a4pr5+a5ps6-a6pt2
       py7(i)= a4pr7-a5ps6-a6pt3
       py8(i)=-a4pr7-a5ps5-a6pt4
      ENDDO
C
      DO i=1,nel
C
        a7pr1=aji7(i)*pr(1)
        a7pr3=aji7(i)*pr(3)
        a7pr5=aji7(i)*pr(5)
        a7pr7=aji7(i)*pr(7)
        a8ps1=aji8(i)*ps(1)
        a8ps2=aji8(i)*ps(2)
        a8ps5=aji8(i)*ps(5)
        a8ps6=aji8(i)*ps(6)
        a9pt1=aji9(i)*pt(1)
        a9pt2=aji9(i)*pt(2)
        a9pt3=aji9(i)*pt(3)
        a9pt4=aji9(i)*pt(4)
       pz1(i)= a7pr1+a8ps1+a9pt1
       pz2(i)=-a7pr1+a8ps2+a9pt2
       pz3(i)= a7pr3-a8ps2+a9pt3
       pz4(i)=-a7pr3-a8ps1+a9pt4
       pz5(i)= a7pr5+a8ps5-a9pt1
       pz6(i)=-a7pr5+a8ps6-a9pt2
       pz7(i)= a7pr7-a8ps6-a9pt3
       pz8(i)=-a7pr7-a8ps5-a9pt4
      ENDDO
C
      RETURN