iqela3.F File Reference

#include "implicit_f.inc"
#include "com01_c.inc"
#include "com04_c.inc"
#include "scr08_a_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	iqela3 (v, w, irect, crst, nsv, irtl, msr, nsn)

Function/Subroutine Documentation

iqela3()

subroutine iqela3	(	dimension(3,numnod), intent(in)	v,
		dimension(3,numnod), intent(inout)	w,
		integer, dimension(4,*), intent(in)	irect,
		dimension(2,*), intent(in)	crst,
		integer, dimension(*), intent(in)	nsv,
		integer, dimension(*), intent(in)	irtl,
		integer, dimension(*), intent(in)	msr,
		integer, intent(in)	nsn )

Definition at line 31 of file iqela3.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      "com01_c.inc"
#include      "com04_c.inc"
#include      "scr08_a_c.inc"
C-----------------------------------------------
C   D u m m y   A r g u m e n t s
C-----------------------------------------------
      INTEGER,INTENT(IN)    :: IRECT(4,*), NSV(*), IRTL(*), MSR(*)
      my_real,INTENT(IN)    :: v(3,numnod), crst(2,*)
      INTEGER, INTENT(IN)   :: NSN
      my_real,INTENT(INOUT) :: w(3,numnod)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER :: NIR, II, N, L, JJ, NN
      my_real :: h(4), ss, tt, vmx, vmy, vmz
C-----------------------------------------------
C   S o u r c e   L i n e s
C-----------------------------------------------
      nir=2
      IF(n2d == 0)nir=4      
      DO  ii=1,nsn
        n        = nsv(ii)
        l        = irtl(ii)
        DO jj=1,nir
          nn     = irect(jj,l)
          ix(jj) = msr(nn)
        ENDDO
        ss       = crst(1,ii)
        tt       = crst(2,ii)
        CALL shapeh(h,ss,tt)
        vmx      = zero
        vmy      = zero
        vmz      = zero     
        DO jj=1,nir
        vmx      = vmx+v(1,ix(jj))*h(jj)
        vmy      = vmy+v(2,ix(jj))*h(jj)
        vmz      = vmz+v(3,ix(jj))*h(jj)
        ENDDO
        !GRID VELOCITY
        w(1,n)   = vmx
        w(2,n)   = vmy
        w(3,n)   = vmz
      ENDDO ! NEXT II
C-----------------------------------------------
      RETURN