c3fcum3.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "com01_c.inc"
#include "com08_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	c3fcum3 (jft, jlt, f, f11, f12, f13, f21, f22, f23, f31, f32, f33, fzero, e1x, e1y, e1z, e2x, e2y, e2z, e3x, e3y, e3z)

Function/Subroutine Documentation

c3fcum3()

subroutine c3fcum3	(	integer	jft,
		integer	jlt,
			f,
			f11,
			f12,
			f13,
			f21,
			f22,
			f23,
			f31,
			f32,
			f33,
			fzero,
			e1x,
			e1y,
			e1z,
			e2x,
			e2y,
			e2z,
			e3x,
			e3y,
			e3z )

Definition at line 29 of file c3fcum3.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   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT, JLT
C     REAL
      my_real f(3,*),fzero(3,4,*),
     .     e1x(mvsiz), e1y(mvsiz), e1z(mvsiz),
     .     e2x(mvsiz), e2y(mvsiz), e2z(mvsiz),
     .     e3x(mvsiz), e3y(mvsiz), e3z(mvsiz),
     .     f11(mvsiz), f12(mvsiz), f13(mvsiz),
     .     f21(mvsiz), f22(mvsiz), f23(mvsiz),
     .     f31(mvsiz), f32(mvsiz), f33(mvsiz)
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
#include      "com01_c.inc"
#include      "com08_c.inc"
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
      my_real g11, g12, g13, g21, g22, g23, g31, g32, g33,
     .     temp1,temp2,temp3
C-----------------------------------------------
C               POUR ISIGI = 5 OU -5
C     CHANGEMENT DE REPERE DE FZERO POUR TT=DT1
C          REPERE GLOBAL ---> REPERE LOCAL
C-----------------------------------------------
      IF(iabs(isigi)==5.AND.ncycle==1.AND.irun==1) THEN
        DO i=jft,jlt
          DO j=1,3
            temp1=fzero(1,j,i)
            temp2=fzero(2,j,i)
            temp3=fzero(3,j,i)
            fzero(1,j,i)=temp1*e1x(i)+temp2*e1y(i)+temp3*e1z(i)
            fzero(2,j,i)=temp1*e2x(i)+temp2*e2y(i)+temp3*e2z(i)
            fzero(3,j,i)=temp1*e3x(i)+temp2*e3y(i)+temp3*e3z(i)
          ENDDO
        ENDDO
      ENDIF
C-----------------------------------------------
C               POUR ISIGI = 5 OU -5
C     PRISE EN COMPTE DES FORCES D EQUILIBRAGE
C      DANS L EVALUATION DES EFFORTS INTERNES
C-----------------------------------------------
      IF((iabs(isigi)==5.AND.ncycle>=1.AND.irun==1).OR.
     .   (iabs(isigi)==5.AND.irun>1)) THEN
        DO i=jft,jlt
          g11 = f11(i)+fzero(1,1,i)
          g21 = f21(i)+fzero(2,1,i)
          g31 = f31(i)+fzero(3,1,i)
          f11(i) = e1x(i)*g11 + e2x(i)*g21 + e3x(i)*g31
          f21(i) = e1y(i)*g11 + e2y(i)*g21 + e3y(i)*g31
          f31(i) = e1z(i)*g11 + e2z(i)*g21 + e3z(i)*g31
          g12 = f12(i)+fzero(1,2,i)
          g22 = f22(i)+fzero(2,2,i)
          g32 = f32(i)+fzero(3,2,i)
          f12(i) = e1x(i)*g12 + e2x(i)*g22 + e3x(i)*g32
          f22(i) = e1y(i)*g12 + e2y(i)*g22 + e3y(i)*g32
          f32(i) = e1z(i)*g12 + e2z(i)*g22 + e3z(i)*g32
          g13 = f13(i)+fzero(1,3,i)
          g23 = f23(i)+fzero(2,3,i)
          g33 = f33(i)+fzero(3,3,i)
          f13(i) = e1x(i)*g13 + e2x(i)*g23 + e3x(i)*g33
          f23(i) = e1y(i)*g13 + e2y(i)*g23 + e3y(i)*g33
          f33(i) = e1z(i)*g13 + e2z(i)*g23 + e3z(i)*g33
        ENDDO
      ENDIF
      IF(iabs(isigi)/=5.OR.tt==0) THEN
        DO i=jft,jlt
            g11 = f11(i)
            g21 = f21(i)
            g31 = f31(i)
            f11(i) = e1x(i)*g11 + e2x(i)*g21 + e3x(i)*g31
            f21(i) = e1y(i)*g11 + e2y(i)*g21 + e3y(i)*g31
            f31(i) = e1z(i)*g11 + e2z(i)*g21 + e3z(i)*g31
            g12 = f12(i)
            g22 = f22(i)
            g32 = f32(i)
            f12(i) = e1x(i)*g12 + e2x(i)*g22 + e3x(i)*g32
            f22(i) = e1y(i)*g12 + e2y(i)*g22 + e3y(i)*g32
            f32(i) = e1z(i)*g12 + e2z(i)*g22 + e3z(i)*g32
            g13 = f13(i)
            g23 = f23(i)
            g33 = f33(i)
            f13(i) = e1x(i)*g13 + e2x(i)*g23 + e3x(i)*g33
            f23(i) = e1y(i)*g13 + e2y(i)*g23 + e3y(i)*g33
            f33(i) = e1z(i)*g13 + e2z(i)*g23 + e3z(i)*g33
        ENDDO
      ENDIF
C
      RETURN