s4therm.F File Reference

#include "implicit_f.inc"
#include "mvsiz_p.inc"
#include "param_c.inc"

Go to the source code of this file.

Functions/Subroutines
subroutine	s4therm (pm, imat, vol, nc1, nc2, nc3, nc4, px1, px2, px3, px4, py1, py2, py3, py4, pz1, pz2, pz3, pz4, dt1, tempnc, tel, heat, fphi, offg, off, nel, theaccfact)

Function/Subroutine Documentation

s4therm()

subroutine s4therm	(		pm,
		integer, intent(in)	imat,
			vol,
		integer, dimension(*)	nc1,
		integer, dimension(*)	nc2,
		integer, dimension(*)	nc3,
		integer, dimension(*)	nc4,
			px1,
			px2,
			px3,
			px4,
			py1,
			py2,
			py3,
			py4,
			pz1,
			pz2,
			pz3,
			pz4,
			dt1,
			tempnc,
			tel,
			heat,
			fphi,
			offg,
			off,
		integer, intent(in)	nel,
		intent(in)	theaccfact )

Definition at line 28 of file s4therm.F.

C-----------------------------------------------
C   I m p l i c i t   T y p e s
C-----------------------------------------------
#include      "implicit_f.inc"
C-----------------------------------------------
C   G l o b a l   P a r a m e t e r s
C-----------------------------------------------
#include      "mvsiz_p.inc"
#include      "param_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) :: IMAT
      INTEGER NC1(*),NC2(*),NC3(*),NC4(*)
      my_real, INTENT(IN) :: theaccfact
      my_real
     .    vol(*), px1(*), px2(*), px3(*), px4(*), 
     .    py1(*), py2(*), py3(*), py4(*),
     .    pz1(*), pz2(*), pz3(*), pz4(*),
     .    tempnc(*), fphi(mvsiz,4), pm(npropm,*),heat(*),
     .    dt1,tel(*),off(*),offg(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J
C     REAL
      my_real
     .  ca, cb, kc, phix, phiy, phiz, a
C===================================================================
 
        ca = pm(75,imat)
        cb = pm(76,imat)
        DO i=1,nel
          IF(off(i)==zero.OR.offg(i)<=zero) cycle
 
C 
C - flux
C      
          phix =  tempnc(nc1(i))*px1(i) + tempnc(nc2(i))*px2(i)  +
     .            tempnc(nc3(i))*px3(i) + tempnc(nc4(i))*px4(i)      
          
          phiy =  tempnc(nc1(i))*py1(i) + tempnc(nc2(i))*py2(i)  +
     .            tempnc(nc3(i))*py3(i) + tempnc(nc4(i))*py4(i)  
      
          phiz =  tempnc(nc1(i))*pz1(i) + tempnc(nc2(i))*pz2(i)  +
     .            tempnc(nc3(i))*pz3(i) + tempnc(nc4(i))*pz4(i)        
C      
          kc = (ca + cb*tel(i))*dt1*vol(i)*theaccfact
          phix = kc*phix
          phiy = kc*phiy
          phiz = kc*phiz
C
C force thermique nodale
C
          a = fourth* heat(i)
          fphi(i,1) = a - (phix*px1(i) + phiy*py1(i) + pz1(i)*phiz)
          fphi(i,2) = a - (phix*px2(i) + phiy*py2(i) + pz2(i)*phiz)
          fphi(i,3) = a - (phix*px3(i) + phiy*py3(i) + pz3(i)*phiz)
          fphi(i,4) = a - (phix*px4(i) + phiy*py4(i) + pz4(i)*phiz)
       ENDDO
C 
      RETURN