s8cumu3.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s8cumu3 (offg, e, nc, stifn, sti, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nel)

Function/Subroutine Documentation

s8cumu3()

subroutine s8cumu3	(		offg,
			e,
		integer, dimension(8,*)	nc,
			stifn,
			sti,
			f11,
			f21,
			f31,
			f12,
			f22,
			f32,
			f13,
			f23,
			f33,
			f14,
			f24,
			f34,
			f15,
			f25,
			f35,
			f16,
			f26,
			f36,
			f17,
			f27,
			f37,
			f18,
			f28,
			f38,
		integer, intent(in)	nel )

Definition at line 28 of file s8cumu3.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, INTENT(IN) :: NEL
      INTEGER NC(8,*)
      my_real
     .   offg(*),e(3,*),stifn(*),sti(*),
     .   f11(*),f21(*),f31(*),f12(*),f22(*),f32(*),
     .   f13(*),f23(*),f33(*),f14(*),f24(*),f34(*),
     .   f15(*),f25(*),f35(*),f16(*),f26(*),f36(*),
     .   f17(*),f27(*),f37(*),f18(*),f28(*),f38(*)
CMasParINCLUDE 'scumu3.intmap.inc'
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I
C-----------------------------------------------
      my_real
     .   off_l
C-----------------------------------------------
      off_l = 0.
      DO i=1,nel
Cf small3b        IF(OFF(I)<1.)OFFG(I) = OFF(I)
        off_l = min(off_l,offg(i))
      ENDDO
      IF(off_l<0.)THEN
        DO i=1,nel
          IF(offg(i)<0.)THEN
            f11(i)=0.
            f21(i)=0.
            f31(i)=0.
            f12(i)=0.
            f22(i)=0.
            f32(i)=0.
            f13(i)=0.
            f23(i)=0.
            f33(i)=0.
            f14(i)=0.
            f24(i)=0.
            f34(i)=0.
            f15(i)=0.
            f25(i)=0.
            f35(i)=0.
            f16(i)=0.
            f26(i)=0.
            f36(i)=0.
            f17(i)=0.
            f27(i)=0.
            f37(i)=0.
            f18(i)=0.
            f28(i)=0.
            f38(i)=0.
            sti(i)=0.
          ENDIF
        ENDDO
      ENDIF
C
C     because 2*Mnodal = 2*Melement/8
      DO i=1,nel
        sti(i)=fourth*sti(i)
      END DO
C
      DO 100 i=1,nel
C
      e(1,nc(1,i))=e(1,nc(1,i))+f11(i)
      e(2,nc(1,i))=e(2,nc(1,i))+f21(i)
      e(3,nc(1,i))=e(3,nc(1,i))+f31(i)
      stifn(nc(1,i))=stifn(nc(1,i))+sti(i)
C
      e(1,nc(7,i))=e(1,nc(7,i))+f17(i)
      e(2,nc(7,i))=e(2,nc(7,i))+f27(i)
      e(3,nc(7,i))=e(3,nc(7,i))+f37(i)
      stifn(nc(7,i))=stifn(nc(7,i))+sti(i)
C
      e(1,nc(2,i))=e(1,nc(2,i))+f12(i)
      e(2,nc(2,i))=e(2,nc(2,i))+f22(i)
      e(3,nc(2,i))=e(3,nc(2,i))+f32(i)
      stifn(nc(2,i))=stifn(nc(2,i))+sti(i)
C
      e(1,nc(8,i))=e(1,nc(8,i))+f18(i)
      e(2,nc(8,i))=e(2,nc(8,i))+f28(i)
      e(3,nc(8,i))=e(3,nc(8,i))+f38(i)
      stifn(nc(8,i))=stifn(nc(8,i))+sti(i)
C
      e(1,nc(3,i))=e(1,nc(3,i))+f13(i)
      e(2,nc(3,i))=e(2,nc(3,i))+f23(i)
      e(3,nc(3,i))=e(3,nc(3,i))+f33(i)
      stifn(nc(3,i))=stifn(nc(3,i))+sti(i)
C
      e(1,nc(5,i))=e(1,nc(5,i))+f15(i)
      e(2,nc(5,i))=e(2,nc(5,i))+f25(i)
      e(3,nc(5,i))=e(3,nc(5,i))+f35(i)
      stifn(nc(5,i))=stifn(nc(5,i))+sti(i)
C
      e(1,nc(4,i))=e(1,nc(4,i))+f14(i)
      e(2,nc(4,i))=e(2,nc(4,i))+f24(i)
      e(3,nc(4,i))=e(3,nc(4,i))+f34(i)
      stifn(nc(4,i))=stifn(nc(4,i))+sti(i)
C
      e(1,nc(6,i))=e(1,nc(6,i))+f16(i)
      e(2,nc(6,i))=e(2,nc(6,i))+f26(i)
      e(3,nc(6,i))=e(3,nc(6,i))+f36(i)
      stifn(nc(6,i))=stifn(nc(6,i))+sti(i)
C
  100 CONTINUE
      RETURN