r4sumg3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	r4sumg3 (jft, jlt, al, kx, r11, r21, r31, ke11, ke12, ke22)

Function/Subroutine Documentation

r4sumg3()

subroutine r4sumg3	(	integer	jft,
		integer	jlt,
			al,
			kx,
			r11,
			r21,
			r31,
			ke11,
			ke12,
			ke22 )

Definition at line 30 of file r4sumg3.F.

C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
#include      "implicit_f.inc"
C-----------------------------------------------
C   D U M M Y   A R G U M E N T S
C-----------------------------------------------
      INTEGER JFT,JLT
      my_real 
     .    al(*),kx(*),r11(*),r21(*),r31(*)
       my_real
     .    ke11(3,3,*),ke22(3,3,*),ke12(3,3,*)
C-----------------------------------------------
C   L O C A L   V A R I A B L E S
C-----------------------------------------------
      INTEGER I, J, EP
C---+----1----+----2----+----3----+----4----+----5----+----6----+----7----+----8
C---------------------------------------
C   TRANS LOCAL-->GLOBAL 
C---------------------------------------
       DO ep=jft,jlt
        ke11(1,1,ep)=kx(ep)*r11(ep)*r11(ep)
        ke11(1,2,ep)=kx(ep)*r11(ep)*r21(ep)
        ke11(1,3,ep)=kx(ep)*r11(ep)*r31(ep)
        ke11(2,2,ep)=kx(ep)*r21(ep)*r21(ep)
        ke11(2,3,ep)=kx(ep)*r21(ep)*r31(ep)
        ke11(3,3,ep)=kx(ep)*r31(ep)*r31(ep)
       ENDDO 
       DO i=1,3 
        DO j=i,3 
         DO ep=jft,jlt
          ke22(i,j,ep)=ke11(i,j,ep)
          ke12(i,j,ep)=-ke11(i,j,ep)
         ENDDO
        ENDDO
       ENDDO
       DO i=1,3 
        DO j=i+1,3 
         DO ep=jft,jlt
          ke12(j,i,ep)=-ke11(i,j,ep)
         ENDDO
        ENDDO
       ENDDO
C
      RETURN