grelem_sav.F File Reference

#include "implicit_f.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	grelem_sav (jft, jlt, gresav, igrth, grth, off, ei, ek, xm, ym, zm, xmas, xcg, ycg, zcg, xxm, yym, zzm, ixx, iyy, izz, ixy, iyz, izx, rei, rek, flag)

Function/Subroutine Documentation

grelem_sav()

subroutine grelem_sav	(	integer	jft,
		integer	jlt,
			gresav,
		integer, dimension(*)	igrth,
		integer, dimension(*)	grth,
			off,
			ei,
			ek,
			xm,
			ym,
			zm,
			xmas,
			xcg,
			ycg,
			zcg,
			xxm,
			yym,
			zzm,
			ixx,
			iyy,
			izz,
			ixy,
			iyz,
			izx,
			rei,
			rek,
		integer	flag )

Definition at line 48 of file grelem_sav.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-----------------------------------------------
#include      "param_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER JFT,JLT,FLAG
      INTEGER IGRTH(*),GRTH(*)
      my_real
     .   gresav(npsav,*),off(*),
     .   xmas(*), ei(*),rei(*),rek(*),
     .   ek(*), xm(*), ym(*), zm(*),
     .   xxm(*), yym(*), zzm(*), 
     .   xcg(*), ycg(*), zcg(*),
     .   ixx(*), iyy(*), izz(*),
     .   ixy(*), iyz(*), izx(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,J
C-----------------------------------------------
      DO i=jft,jlt
        IF (igrth(i) /= igrth(i+1)) THEN
          IF(off(i) /= zero)THEN
            DO j = igrth(i),igrth(i+1)-1
              gresav(1,grth(j)) = gresav(1,grth(j)) + ei(i)
              gresav(2,grth(j)) = gresav(2,grth(j)) + ek(i)
              gresav(3,grth(j)) = gresav(3,grth(j)) + xm(i)
              gresav(4,grth(j)) = gresav(4,grth(j)) + ym(i)
              gresav(5,grth(j)) = gresav(5,grth(j)) + zm(i)
              gresav(6,grth(j)) = gresav(6,grth(j)) + xmas(i)
            ENDDO
          ENDIF
        ENDIF
      ENDDO
      IF (flag /= 0)THEN
        DO i=jft,jlt
          IF (igrth(i) /= igrth(i+1)) THEN
            DO j = igrth(i),igrth(i+1)-1
              gresav(9,grth(j))  = gresav(9,grth(j))  + xcg(i)
              gresav(10,grth(j)) = gresav(10,grth(j)) + ycg(i)
              gresav(11,grth(j)) = gresav(11,grth(j)) + zcg(i)
              gresav(12,grth(j)) = gresav(12,grth(j)) + xxm(i)
              gresav(13,grth(j)) = gresav(13,grth(j)) + yym(i)
              gresav(14,grth(j)) = gresav(14,grth(j)) + zzm(i)
              gresav(15,grth(j)) = gresav(15,grth(j)) + ixx(i)
              gresav(16,grth(j)) = gresav(16,grth(j)) + iyy(i)
              gresav(17,grth(j)) = gresav(17,grth(j)) + izz(i)
              gresav(18,grth(j)) = gresav(18,grth(j)) + ixy(i)
              gresav(19,grth(j)) = gresav(19,grth(j)) + iyz(i)
              gresav(20,grth(j)) = gresav(20,grth(j)) + izx(i)
              gresav(21,grth(j)) = gresav(21,grth(j)) + rei(i)
              gresav(22,grth(j)) = gresav(22,grth(j)) + rek(i)
            ENDDO
          ENDIF
        ENDDO
      ENDIF
      RETURN