cbatempel.F File Reference

#include "implicit_f.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	cbatempel (jft, jlt, ng, ixc, temp, tempel)

Function/Subroutine Documentation

cbatempel()

subroutine cbatempel	(	integer	jft,
		integer	jlt,
		integer	ng,
		integer, dimension(nixc,*)	ixc,
			temp,
			tempel )

Definition at line 30 of file cbatempel.F.

      use element_mod , only : nixc
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 JFT, JLT,NG,IXC(NIXC,*)
      my_real temp(*),tempel(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I,KSI ,ETA ,SF1,SF2 ,SF3 ,SF4
      my_real 
     .   vpg(2,4),pg1,pg
       parameter(pg=.577350269189626)
       parameter(pg1=-.577350269189626)
C--------------------------
C     INITIALISATION
C-------------------------
      DATA vpg/pg1,pg1,pg,pg1,pg,pg,pg1,pg/
C-----------------------------------------------
C --- Temperature in Gauss point NG
       ksi = vpg(1,ng)
       eta = vpg(2,ng)
 
       sf1 = (1-ksi)*(1-eta)
       sf2 = (1+ksi)*(1-eta)
       sf3 = (1+ksi)*(1+eta)
       sf4 = (1-ksi)*(1+eta)
 
        DO i=jft,jlt          
          tempel(i) = fourth * (sf1 *temp(ixc(2,i)) + sf2*temp(ixc(3,i)) + 
     .                          sf3 *temp(ixc(4,i)) + sf4*temp(ixc(5,i)) )
       ENDDO
C
      RETURN