s6ctherm.F File Reference

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

Go to the source code of this file.

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

Function/Subroutine Documentation

s6ctherm()

subroutine s6ctherm	(	dimension(npropm,nummat), intent(in)	pm,
		integer, intent(in)	imat,
		dimension(mvsiz), intent(in)	vol,
		integer, dimension(mvsiz), intent(in)	nc1,
		integer, dimension(mvsiz), intent(in)	nc2,
		integer, dimension(mvsiz), intent(in)	nc3,
		integer, dimension(mvsiz), intent(in)	nc4,
		integer, dimension(mvsiz), intent(in)	nc5,
		integer, dimension(mvsiz), intent(in)	nc6,
		dimension(mvsiz), intent(in)	px1,
		dimension(mvsiz), intent(in)	px2,
		dimension(mvsiz), intent(in)	px3,
		dimension(mvsiz), intent(in)	px4,
		dimension(mvsiz), intent(in)	py1,
		dimension(mvsiz), intent(in)	py2,
		dimension(mvsiz), intent(in)	py3,
		dimension(mvsiz), intent(in)	py4,
		dimension(mvsiz), intent(in)	pz1,
		dimension(mvsiz), intent(in)	pz2,
		dimension(mvsiz), intent(in)	pz3,
		dimension(mvsiz), intent(in)	pz4,
		intent(in)	dt1,
		dimension(numnod), intent(in)	tempnc,
		dimension(mvsiz), intent(in)	tel,
		dimension(mvsiz), intent(in)	heat,
		dimension(mvsiz,6), intent(inout)	fphi,
		dimension(mvsiz), intent(in)	offg,
		dimension(mvsiz), intent(in)	off,
		integer, intent(in)	nel,
		intent(in)	theaccfact )

Definition at line 28 of file s6ctherm.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"
#include      "com04_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, INTENT(IN) :: NC1(MVSIZ),NC2(MVSIZ),
     .   NC3(MVSIZ),NC4(MVSIZ),NC5(MVSIZ),NC6(MVSIZ)
      my_real, INTENT(IN) :: theaccfact
      my_real, INTENT(IN) :: dt1
      my_real, INTENT(IN) :: vol(mvsiz), px1(mvsiz), py1(mvsiz),pz1(mvsiz),
     .   px2(mvsiz), py2(mvsiz), pz2(mvsiz),
     .   px3(mvsiz), py3(mvsiz),pz3(mvsiz),
     .   px4(mvsiz), py4(mvsiz),pz4(mvsiz),
     .   tempnc(numnod), pm(npropm,nummat),
     .   heat(mvsiz), tel(mvsiz),off(mvsiz),offg(mvsiz)
      my_real, INTENT(INOUT) :: fphi(mvsiz,6)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I
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 .AND. offg(i)>zero) THEN 
c          
          ! Flux initialization 
          phix = zero
          phiy = zero
          phiz = zero
c
          ! Computing flux
          phix =  tempnc(nc1(i))*(px1(i)-px4(i)) + tempnc(nc2(i))*(px2(i)-px4(i)) +
     .            tempnc(nc3(i))*(px3(i)-px4(i)) + tempnc(nc4(i))*(px1(i)+px4(i)) +
     .            tempnc(nc5(i))*(px2(i)+px4(i)) + tempnc(nc6(i))*(px3(i)+px4(i))   
          
          phiy =  tempnc(nc1(i))*(py1(i)-py4(i)) + tempnc(nc2(i))*(py2(i)-py4(i)) +
     .            tempnc(nc3(i))*(py3(i)-py4(i)) + tempnc(nc4(i))*(py1(i)+py4(i)) +
     .            tempnc(nc5(i))*(py2(i)+py4(i)) + tempnc(nc6(i))*(py3(i)+py4(i)) 
      
          phiz =  tempnc(nc1(i))*(pz1(i)-pz4(i)) + tempnc(nc2(i))*(pz2(i)-pz4(i)) +
     .            tempnc(nc3(i))*(pz3(i)-pz4(i)) + tempnc(nc4(i))*(pz1(i)+pz4(i)) +
     .            tempnc(nc5(i))*(pz2(i)+pz4(i)) + tempnc(nc6(i))*(pz3(i)+pz4(i)) 
c
          kc = (ca + cb*tel(i))*dt1*vol(i)*theaccfact           
          phix = kc*phix
          phiy = kc*phiy
          phiz = kc*phiz 
c
          ! Computing nodal thermic forces
          a = heat(i)*one_over_6
          fphi(i,1) = fphi(i,1) +
     .        a - (phix*(px1(i)-px4(i)) + phiy*(py1(i)-py4(i)) + phiz*(pz1(i)-pz4(i)))
          fphi(i,2) = fphi(i,2) +
     .        a - (phix*(px2(i)-px4(i)) + phiy*(py2(i)-py4(i)) + phiz*(pz2(i)-pz4(i)))
          fphi(i,3) = fphi(i,3) +
     .        a - (phix*(px3(i)-px4(i)) + phiy*(py3(i)-py4(i)) + phiz*(pz3(i)-pz4(i)))
          fphi(i,4) = fphi(i,4) +
     .        a - (phix*(px1(i)+px4(i)) + phiy*(py1(i)+py4(i)) + phiz*(pz1(i)+pz4(i)))
          fphi(i,5) = fphi(i,5) +
     .        a - (phix*(px2(i)+px4(i)) + phiy*(py2(i)+py4(i)) + phiz*(pz2(i)+pz4(i)))
          fphi(i,6) = fphi(i,6) +
     .        a - (phix*(px3(i)+px4(i)) + phiy*(py3(i)+py4(i)) + phiz*(pz3(i)+pz4(i)))
        ENDIF
      ENDDO
C