1 | c.......................................................................
   2 | c
   3 | c     calcul de la derivee de l'emittance monochromatique du corps noir 
   4 | c     par rapport a la temperature 
   5 | c     a un nombre d'onde et une temperature donnee
   6 | c
   7 | c     en entree : * le nombre d'onde blanu (m-1)
   8 | c                 * la temperature blat (K)
   9 | c
  10 | c     en sortie : * la derivee de l' emittance dblae
  11 | c
  12 | c.......................................................................
  13 |       subroutine black_dt(blanu,blat,dblae)
  14 | c.......................................................................
  15 |       implicit double precision (a-h,o-z)
  16 |       include 'cecile.inc'  !pour pi
  17 | c.......................................................................
  18 |       double precision sigma,c0,h,cbol,rind,c,c1,c2
  19 | !      double precision pi
  20 | c.......................................................................
  21 |       sigma=5.6693D-8
  22 | !!      pi=datan(1.D+0)*4.D+0
  23 |       c0=2.9979D+8
  24 |       h=6.6262D-34
  25 |       cbol=1.3806D-23
  26 |       rind=1.D+0
  27 |       c=c0/rind
  28 |       c1=h*(c**2)
  29 |       c2=h*c/cbol
  30 | c.......................................................................
  31 |       if (blat .eq. 0.d+0) then
  32 |          dblae = 0.d+0
  33 |          else
  34 | 
  35 |          dblae=2.D+0*pi*c1*c2*blanu**4*(dexp(c2*blanu/blat))/
  36 |      &         (blat*(dexp(c2*blanu/blat)-1.D+0))**2
  37 |       endif
  38 | 
  39 |       if ((c2*blanu/blat).gt. 700.d+0) then
  40 |          !cas rare toujour provoquer par faible temperature
  41 |          !exp(tre grand nombre) = infinity
  42 |          dblae = 0.d+0
  43 |       endif
  44 | 
  45 | !      print*, blanu, (c2*blanu/blat),dexp(c2*blanu/blat),dblae
  46 | !      if ((in.eq.0).and.(iin.eq.39)) then
  47 | !      print *, in, iin, dblae, blat, blanu
  48 | !      endif
  49 | 
  50 | c     ca fait zero en zero
  51 | c     les erreurs due a une temperature proche de zero ne sont pas evitees
  52 | c     le probleme ne s'est jamais pose: amaury le 1999 IV 14 
  53 | c.......................................................................
  54 |       return
  55 |       end
  56 | c.......................................................................
  57 | 
  58 | 
  59 | 
  60 |