fcumu3p.F File Reference

#include "implicit_f.inc"
#include "comlock.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	fcumu3p (offg, sti, fsky, iads, f11, f21, f31, f12, f22, f32, f13, f23, f33, f14, f24, f34, f15, f25, f35, f16, f26, f36, f17, f27, f37, f18, f28, f38, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, nel, nft)

Function/Subroutine Documentation

fcumu3p()

subroutine fcumu3p	(		offg,
			sti,
			fsky,
		integer, dimension(8,*)	iads,
			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, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
		integer, dimension(*)	nc5,
		integer, dimension(*)	nc6,
		integer, dimension(*)	nc7,
		integer, dimension(*)	nc8,
		integer, intent(in)	nel,
		integer, intent(in)	nft )

Definition at line 28 of file fcumu3p.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
#include      "comlock.inc"
C-----------------------------------------------
C   C o m m o n   B l o c k s
C-----------------------------------------------
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER, INTENT(IN) :: NEL
      INTEGER, INTENT(IN) :: NFT
      INTEGER NC1(*), NC2(*), NC3(*), NC4(*), NC5(*), NC6(*), NC7(*), NC8(*)
      my_real  offg(*),fsky(3,*),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(*)
C
      INTEGER IADS(8,*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, II, K, J
      my_real off_l
C-----------------------------------------------
      off_l = 0.
      DO i=1,nel
        off_l = min(off_l,offg(i))
      ENDDO
      IF(off_l < zero)THEN
        DO i=1,nel
         IF(offg(i) >= zero)cycle
           f11(i)=zero
           f21(i)=zero
           f31(i)=zero
           f12(i)=zero
           f22(i)=zero
           f32(i)=zero
           f13(i)=zero
           f23(i)=zero
           f33(i)=zero
           f14(i)=zero
           f24(i)=zero
           f34(i)=zero
           f15(i)=zero
           f25(i)=zero
           f35(i)=zero
           f16(i)=zero
           f26(i)=zero
           f36(i)=zero
           f17(i)=zero
           f27(i)=zero
           f37(i)=zero
           f18(i)=zero
           f28(i)=zero
           f38(i)=zero
           sti(i)=zero
        ENDDO
      ENDIF
C
C     because 2*Mnodal = 2*Melement/8
      DO i=1,nel
        sti(i)=fourth*sti(i)
      END DO
C
cc       RETURN
C
        DO i=1,nel
          ii=i + nft
          k = iads(1,ii)
c          F21(I)=F11(I)
c          F11(I) = ZERO
                      
          fsky(1,k)=f11(i)
          fsky(2,k)=f21(i)
          fsky(3,k)=f31(i) 
C
          k = iads(7,ii)
          fsky(1,k)=f17(i)
          fsky(2,k)=f27(i)
          fsky(3,k)=f37(i)       
C
          k = iads(8,ii)
          fsky(1,k)=f18(i)
          fsky(2,k)=f28(i)
          fsky(3,k)=f38(i)    
C
          k = iads(3,ii)
          fsky(1,k)=f13(i)
          fsky(2,k)=f23(i)
          fsky(3,k)=f33(i)    
C
          k = iads(5,ii)
          fsky(1,k)=f15(i)
          fsky(2,k)=f25(i)
          fsky(3,k)=f35(i)    
C
          k = iads(4,ii)
          fsky(1,k)=f14(i)
          fsky(2,k)=f24(i)
          fsky(3,k)=f34(i)    
C
          k = iads(6,ii)
          fsky(1,k)=f16(i)
          fsky(2,k)=f26(i)
          fsky(3,k)=f36(i) 
C          
          k = iads(2,ii)
          fsky(1,k) = f12(i)
          fsky(2,k)=f22(i)
          fsky(3,k)=f32(i)   
        ENDDO
        
      RETURN