s8ztempel.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s8ztempel (lft, llt, ni, ksi, eta, zeta, nc1, nc2, nc3, nc4, nc5, nc6, nc7, nc8, temp, tempel)

Function/Subroutine Documentation

s8ztempel()

subroutine s8ztempel	(	integer	lft,
		integer	llt,
			ni,
			ksi,
			eta,
			zeta,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
		integer, dimension(*)	nc5,
		integer, dimension(*)	nc6,
		integer, dimension(*)	nc7,
		integer, dimension(*)	nc8,
			temp,
			tempel )

Definition at line 31 of file s8ztempel.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 LFT,LLT,NC1(*), NC2(*), NC3(*), NC4(*), 
     .        NC5(*), NC6(*), NC7(*), NC8(*)
C     REAL
      my_real
     .   ksi,eta,zeta,ni(*),temp(*),tempel(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
C     REAL
      my_real
     .   rp,sp,tp,sm,rm,tm
C-----------------------------------------------
C
C isoparametric shape function
C
      rp = one + ksi            
      sp = one + eta            
      tp = one + zeta           
      rm = one - ksi            
      sm = one - eta            
      tm = one - zeta
      ni(1) = one_over_8*rm*sm*tm
      ni(2) = one_over_8*rp*sm*tm
      ni(3) = one_over_8*rp*sp*tm
      ni(4) = one_over_8*rm*sp*tm
      ni(5) = one_over_8*rm*sm*tp
      ni(6) = one_over_8*rp*sm*tp
      ni(7) = one_over_8*rp*sp*tp
      ni(8) = one_over_8*rm*sp*tp
      DO i=lft,llt
         tempel(i)=   ni(1)*temp(nc1(i)) +  ni(2)*temp(nc2(i)) +
     .                ni(3)*temp(nc3(i)) +  ni(4)*temp(nc4(i)) +
     .                ni(5)*temp(nc5(i)) +  ni(6)*temp(nc6(i)) +
     .                ni(7)*temp(nc7(i)) +  ni(8)*temp(nc8(i))  
      ENDDO
C-----------
      RETURN