m26th.F File Reference

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

Go to the source code of this file.

Functions/Subroutines
subroutine	m26th (mat, rho, t, xk, pm, sesame, z, nel, nft)

Function/Subroutine Documentation

m26th()

subroutine m26th	(	integer, dimension(*)	mat,
			rho,
			t,
			xk,
			pm,
			sesame,
			z,
		integer, intent(in)	nel,
		integer, intent(in)	nft )

Definition at line 30 of file m26th.F.

C----------------------------------------
C     CALCUL DE LA CONDUCTIVITE THERMIQUE
C----------------------------------------
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      "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) :: NFT
      INTEGER MAT(*)
      my_real
     .   rho(*), t(*), xk(*), pm(npropm,*), sesame(*), z(*)
C-----------------------------------------------
C   L o c a l   V a r i a b l e s
C-----------------------------------------------
      INTEGER I, J, MX, NR, NT, IDR, IDT, IDQ
      my_real
     .   zz, deltat, xlamb, stefan, rl, tl, opac, dxdr, rossel, xkr
      real*8 ne, na, atom
C-----------------------------------------------
      DATA na /6.0225e+23/
C-----------------------------------------------
      !----------------------------
      !     CONDUCTION ELECTRONIQUE
      !----------------------------
      DO i=1,nel
        j=i+nft
        mx = mat(i)
        zz = max(em10,z(i))
        deltat= threep44 * zep26 * log(zz) / zz
        IF(deltat>zero)THEN
          deltat= one / ( one +  deltat)
        ELSE
          deltat= zero
        ENDIF
        atom  = pm(37,mx)
        ne    = rho(i)*na*zz/atom
        xlamb = pm(36,mx)*t(i)**three_half/sqrt(ne)
        xlamb = max(one,xlamb)
        xlamb = max(em10, log(xlamb))
        xk(j) = xk(j) + zep4*deltat*pm(35,mx)*t(i)**twop5 / (zz*xlamb)
      ENDDO
 
      !----------------------------
      !     RADIATION
      !----------------------------
      DO i=1,nel
        j=i+nft
        mx = mat(i)
        stefan = pm(51,mx)
        IF(stefan>zero.AND.t(i)>ep04)THEN
          nr     = nint(pm(48,mx))
          nt     = nint(pm(49,mx))
          idr    = nint(pm(50,mx))
          idt = idr + nr
          idq = idt + nt
          rl = log10(rho(i))
          tl = log10(t(i))
          CALL mintp_rt(sesame(idr),nr,
     +        sesame(idt),nt,sesame(idq),rl,tl,opac,dxdr)
          opac = ten**opac
          rossel = one / ( rho(i) * opac)
          xkr = sixteen * stefan * t(i)**3 * rossel * third
          xk(j) = xk(j) + xkr
        ENDIF
        xk(j) = min(xk(j),pm(52,mx))
      ENDDO
C
      RETURN