cdk6fcum3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	cdk6fcum3 (jft, jlt, r11, r12, r13, r21, r22, r23, r31, r32, r33, f11, f12, f13, f21, f22, f23, f31, f32, f33, f14, f15, f16, f24, f25, f26, f34, f35, f36)

Function/Subroutine Documentation

cdk6fcum3()

subroutine cdk6fcum3	(	integer	jft,
		integer	jlt,
			r11,
			r12,
			r13,
			r21,
			r22,
			r23,
			r31,
			r32,
			r33,
			f11,
			f12,
			f13,
			f21,
			f22,
			f23,
			f31,
			f32,
			f33,
			f14,
			f15,
			f16,
			f24,
			f25,
			f26,
			f34,
			f35,
			f36 )

Definition at line 28 of file cdk6fcum3.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#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
     .   r11(*),r12(*),r13(*),r21(*),r22(*),r23(*),
     .   r31(*),r32(*),r33(*),
     .   f11(*), f12(*), f13(*),f21(*), f22(*), f23(*),
     .   f31(*), f32(*), f33(*),
     .   f14(*), f15(*), f16(*),f24(*), f25(*), f26(*),
     .   f34(*), f35(*), f36(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
      my_real lx,ly,lz
C-----------------------------------------------
C   L o c a l Coordinate-> Global 
C-----------------------------------------------
C
      DO i=jft,jlt
       lx =r11(i)*f11(i)+ r12(i)*f21(i)+r13(i)*f31(i)
       ly =r21(i)*f11(i)+ r22(i)*f21(i)+r23(i)*f31(i)
       lz =r31(i)*f11(i)+ r32(i)*f21(i)+r33(i)*f31(i)
       f11(i) =lx
       f21(i) =ly
       f31(i) =lz
       lx =r11(i)*f12(i)+ r12(i)*f22(i)+r13(i)*f32(i)
       ly =r21(i)*f12(i)+ r22(i)*f22(i)+r23(i)*f32(i)
       lz =r31(i)*f12(i)+ r32(i)*f22(i)+r33(i)*f32(i)
       f12(i) =lx
       f22(i) =ly
       f32(i) =lz
       lx =r11(i)*f13(i)+ r12(i)*f23(i)+r13(i)*f33(i)
       ly =r21(i)*f13(i)+ r22(i)*f23(i)+r23(i)*f33(i)
       lz =r31(i)*f13(i)+ r32(i)*f23(i)+r33(i)*f33(i)
       f13(i) =lx
       f23(i) =ly
       f33(i) =lz
      ENDDO
C
      DO i=jft,jlt
       lx =r11(i)*f14(i)+ r12(i)*f24(i)+r13(i)*f34(i)
       ly =r21(i)*f14(i)+ r22(i)*f24(i)+r23(i)*f34(i)
       lz =r31(i)*f14(i)+ r32(i)*f24(i)+r33(i)*f34(i)
       f14(i) =lx
       f24(i) =ly
       f34(i) =lz
       lx =r11(i)*f15(i)+ r12(i)*f25(i)+r13(i)*f35(i)
       ly =r21(i)*f15(i)+ r22(i)*f25(i)+r23(i)*f35(i)
       lz =r31(i)*f15(i)+ r32(i)*f25(i)+r33(i)*f35(i)
       f15(i) =lx
       f25(i) =ly
       f35(i) =lz
       lx =r11(i)*f16(i)+ r12(i)*f26(i)+r13(i)*f36(i)
       ly =r21(i)*f16(i)+ r22(i)*f26(i)+r23(i)*f36(i)
       lz =r31(i)*f16(i)+ r32(i)*f26(i)+r33(i)*f36(i)
       f16(i) =lx
       f26(i) =ly
       f36(i) =lz
      ENDDO
C
      RETURN