qcumu2p.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	qcumu2p (f11, f12, f21, f22, ax1, ax2, t11, t12, t13, t14, t21, t22, t23, t24, fsky, fskyv, iadq, sti, nel, nft)

Function/Subroutine Documentation

qcumu2p()

subroutine qcumu2p	(		f11,
			f12,
			f21,
			f22,
			ax1,
			ax2,
			t11,
			t12,
			t13,
			t14,
			t21,
			t22,
			t23,
			t24,
			fsky,
			fskyv,
		integer, dimension(4,*)	iadq,
			sti,
		integer, intent(in)	nel,
		integer, intent(in)	nft )

Definition at line 28 of file qcumu2p.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
     .   fskyv(lsky,8),fsky(8,lsky),sti(*)
C     REAL
      my_real
     .   f11(*), f12(*), f21(*), f22(*), ax1(*), ax2(*),
     .   t11(*), t12(*), t13(*), t14(*), t21(*), t22(*), t23(*), t24(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER  I, II, K
C-------------------------
C     ASSEMBLE
C-------------------------
C
      IF(ivector==1) THEN
        DO i=1,nel
C
          ii=i+nft
          k = iadq(1,ii)
          fskyv(k,2)=-f11(i)-ax1(i)-t11(i)
          fskyv(k,3)=-f21(i)+ax2(i)-t21(i)
          fskyv(k,7)=sti(i)
C
          k = iadq(2,ii)
          fskyv(k,2)=-f12(i)-ax1(i)-t12(i)
          fskyv(k,3)=-f22(i)+ax2(i)-t22(i)
          fskyv(k,7)=sti(i)
C
          k = iadq(3,ii)
          fskyv(k,2)=+f11(i)-ax1(i)-t13(i)
          fskyv(k,3)=+f21(i)+ax2(i)-t23(i)
          fskyv(k,7)=sti(i)
C
          k = iadq(4,ii)
          fskyv(k,2)=+f12(i)-ax1(i)-t14(i)
          fskyv(k,3)=+f22(i)+ax2(i)-t24(i)
          fskyv(k,7)=sti(i)
C
        ENDDO
C
      ELSE
        DO i=1,nel
C
          ii=i+nft
          k = iadq(1,ii)
          fsky(2,k)=-f11(i)-ax1(i)-t11(i)
          fsky(3,k)=-f21(i)+ax2(i)-t21(i)
          fsky(7,k) = sti(i)
C
          k = iadq(2,ii)
          fsky(2,k)=-f12(i)-ax1(i)-t12(i)
          fsky(3,k)=-f22(i)+ax2(i)-t22(i)
          fsky(7,k) = sti(i)
C
          k = iadq(3,ii)
          fsky(2,k)=+f11(i)-ax1(i)-t13(i)
          fsky(3,k)=+f21(i)+ax2(i)-t23(i)
          fsky(7,k) = sti(i)
C
          k = iadq(4,ii)
          fsky(2,k)=+f12(i)-ax1(i)-t14(i)
          fsky(3,k)=+f22(i)+ax2(i)-t24(i)
          fsky(7,k) = sti(i)
C
        ENDDO
C
      ENDIF
C
      RETURN