1 | include(dom.inc)
   2 | 
   3 |       DOM_REAL FUNCTION BLAE(blanu,blat)
   4 | 
   5 |       IMPLICIT NONE
   6 | 
   7 |       DOM_REAL, INTENT(IN) :: blanu, blat
   8 |       DOM_REAL:: sigma,pi,c0,hp,cbol,rind,c,c1,c2
   9 | 
  10 |       sigma=5.667E-08
  11 |       pi=3.141592654
  12 |       c0=2.9979E8
  13 |       hp=6.626E-34
  14 |       cbol=1.3805E-23
  15 |       rind=1.
  16 |       c=c0/rind
  17 |       c1=hp*(c**2)
  18 |       c2=hp*c/cbol
  19 | 
  20 |       IF (blat .eq. 0.) THEN
  21 |          blae = 0.
  22 |       ELSE
  23 |          blae=2.*pi*c1*blanu**3/(exp(c2*blanu/blat)-1.)
  24 |       ENDIF
  25 | 
  26 | 
  27 |       RETURN
  28 |       END FUNCTION BLAE