q4cumu2p.F File Reference

#include "implicit_f.inc"
#include "parit_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	q4cumu2p (fsky, fskyv, iadq, fy1, fz1, fy2, fz2, fy3, fz3, fy4, fz4, stim, nel, nft)

Function/Subroutine Documentation

q4cumu2p()

subroutine q4cumu2p	(		fsky,
			fskyv,
		integer, dimension(4,*)	iadq,
			fy1,
			fz1,
			fy2,
			fz2,
			fy3,
			fz3,
			fy4,
			fz4,
			stim,
		integer, intent(in)	nel,
		integer, intent(in)	nft )

Definition at line 28 of file q4cumu2p.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      "parit_c.inc"
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 IADQ(4,*)
C     REAL
      my_real
     .    fsky(8,lsky), fskyv(lsky,8),
     .    fy1(*), fz1(*), fy2(*), fz2(*),
     .    fy3(*), fz3(*), fy4(*), fz4(*),
     .    stim(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I, II, K
C-----------------------------------------------
C   S o u r c e  L i n e s
C-----------------------------------------------
C
      DO i=1,nel
        stim(i)=fourth*stim(i)
      END DO
C
      IF(ivector==1) THEN
        DO i=1,nel
          ii = i + nft
C
          k = iadq(1,ii)
          fskyv(k,2) = fy1(i)
          fskyv(k,3) = fz1(i)
          fskyv(k,7) = stim(i)
C
          k = iadq(2,ii)
          fskyv(k,2) = fy2(i)
          fskyv(k,3) = fz2(i)
          fskyv(k,7) = stim(i)
C
          k = iadq(3,ii)
          fskyv(k,2) = fy3(i)
          fskyv(k,3) = fz3(i)
          fskyv(k,7) = stim(i)
C
          k = iadq(4,ii)
          fskyv(k,2) = fy4(i)
          fskyv(k,3) = fz4(i)
          fskyv(k,7) = stim(i)
        ENDDO
      ELSE
        DO i=1,nel
          ii = i + nft
C
          k = iadq(1,ii)
          fsky(2,k) = fy1(i)
          fsky(3,k) = fz1(i)
          fsky(7,k) = stim(i)
C
          k = iadq(2,ii)
          fsky(2,k) = fy2(i)
          fsky(3,k) = fz2(i)
          fsky(7,k) = stim(i)
C
          k = iadq(3,ii)
          fsky(2,k) = fy3(i)
          fsky(3,k) = fz3(i)
          fsky(7,k) = stim(i)
C
          k = iadq(4,ii)
          fsky(2,k) = fy4(i)
          fsky(3,k) = fz4(i)
          fsky(7,k) = stim(i)
        ENDDO
      ENDIF
C
      RETURN