1 | include(dom.inc)
   2 | 
   3 |        SUBROUTINE EMISSIV_SNB(nfcelt,celldata,Lb,Lo,epsil,Tf,WVNB_SI,   &
   4 |      &                        ncells,nfacemax)
   5 | 
   6 |        IMPLICIT NONE
   7 | 
   8 |        DOM_REAL , PARAMETER :: pi = 3.141592654
   9 |        DOM_INT  :: ncells,nfacemax,ielt,m
  10 |        DOM_INT, DIMENSION (ncells) :: nfcelt
  11 | 
  12 |        DOM_REAL, DIMENSION(ncells) :: Lb
  13 |        DOM_REAL, DIMENSION(ncells,8) :: celldata
  14 |        DOM_REAL, DIMENSION(ncells,nfacemax) :: Lo, epsil, Tf
  15 |        DOM_REAL  :: blae
  16 |        DOM_REAL  :: WVNB_SI
  17 | 
  18 | !       print*, " Setting black luminance in domain and bndry"
  19 |        DO ielt=1,ncells
  20 |       
  21 | !          print*, "cell: ", ielt, "/",ncells
  22 |           DO m=1,nfcelt(ielt)
  23 |             
  24 |              !------------------------------!
  25 |              ! Luminances noires aux parois !
  26 |              !------------------------------!
  27 | 
  28 | !             print*, " face: ", m,"/",nfcelt(ielt)
  29 |              Lo(ielt,m)=epsil(ielt,m)/pi*blae(WVNB_SI,Tf(ielt,m))
  30 |             
  31 |           ENDDO
  32 |       
  33 |           !----------------------------------!
  34 |           ! Luminances noires dans le milieu !
  35 |           !----------------------------------!
  36 | 
  37 |           Lb(ielt)=blae(WVNB_SI,celldata(ielt,1))/pi
  38 | !          print*, "Lb     =", Lb(ielt)
  39 | !          print*, "T      =", celldata(ielt,1)
  40 | !          print*, "blae   =", blae(WVNB_SI,celldata(ielt,1))
  41 | !          print*, "WVNB_SI=", WVNB_SIB
  42 | 
  43 |        ENDDO
  44 | 
  45 |        END SUBROUTINE EMISSIV_SNB