1 | include(dom.inc)
2 |
3 | SUBROUTINE EMISSIV(nfcelt, celldata, Lb, Lo, epsil, Tf, ncells, &
4 | & nfacemax)
5 |
6 | IMPLICIT NONE
7 |
8 | DOM_INT :: ielt,m,ncells,nfacemax
9 | DOM_INT, DIMENSION (ncells) :: nfcelt
10 |
11 | DOM_REAL,DIMENSION(ncells,8) :: celldata
12 | DOM_REAL,DIMENSION(ncells) :: Lb
13 | DOM_REAL, DIMENSION (ncells,nfacemax) :: Lo, epsil, Tf
14 | DOM_REAL :: planck
15 |
16 |
17 | DO ielt=1,ncells
18 |
19 | !------------------------------!
20 | ! Luminances noires aux parois !
21 | !------------------------------!
22 |
23 | DO m=1,nfcelt(ielt)
24 | Lo(ielt,m)=epsil(ielt,m)*planck(Tf(ielt,m))
25 | ENDDO
26 |
27 | !----------------------------------!
28 | ! Luminances noires dans le milieu !
29 | !----------------------------------!
30 |
31 | Lb(ielt)=planck(celldata(ielt,1))
32 |
33 | ENDDO
34 |
35 | END SUBROUTINE EMISSIV
emissiv.F could be called by: