1 | include(dom.inc)
   2 | 
   3 |       DOM_REAL FUNCTION BLAE(nu,T)
   4 | 
   5 | !     ----------------------------------------------------------------!
   6 | !     les erreurs due a une temperature proche de zero ne sont pas    !
   7 | !     evitees. le probleme ne s'est jamais pose: amaury le 1999 II 22 !
   8 | !     ----------------------------------------------------------------!
   9 | 
  10 |         IMPLICIT NONE
  11 | 
  12 |         include 'dom_constants.h'
  13 | 
  14 | !       IN
  15 |         DOM_REAL, INTENT(IN) :: nu, T
  16 | 
  17 | !       LOCAL
  18 |         DOM_REAL:: rind, c, c1, c2, nu3
  19 | 
  20 | !---    Refraction index 'rind' to be entered as input !
  21 |         rind=1.
  22 | 
  23 |         c=c0/rind
  24 | 
  25 |         c1  = Planckcst * (c*c)
  26 |         c2  = Planckcst * c / Boltzmann
  27 |         nu3 = nu*nu*nu
  28 | 
  29 |         IF (T.eq.0.) THEN
  30 |           blae = 0.
  31 |         ELSE
  32 |           blae= 2.*pi*c1*nu3 / ( exp(c2*nu/T) - 1. )
  33 |         ENDIF
  34 | 
  35 |       RETURN
  36 |       END FUNCTION BLAE