1 |       subroutine d2fds1ds2(k1,k2,a,ecg,phicg,phi1,phi2,der2)
   2 | 
   3 | c     Calcule d2[f(a)]/ds1.ds2*1/f(a)
   4 | c
   5 | c     Entrées : k1,k2 : coefficients d'absorption moyens dans les mailles 1 et 2
   6 | c               a : coefficient d'absorption adimensionnel
   7 | c               ecg : épaisseur optique moyenne (Par Curtis-Godson) entre les 2 points
   8 | c               phicg : phi moyen, par Curtis-Godson, entre les deux points
   9 | c               phi1, phi2 : phi de la maille 1, de la maille 2
  10 | c
  11 | c      Sortie : der2 : dérivée seconde de f(a) par rapport à s1 et s2, prise en s1=s2=0,
  12 | c                      divisée par la valeur de f(a)
  13 | c
  14 | c
  15 | c    Pt d'émission   |                              |                 Pt d'absorption
  16 | c        +-----------|------------------------------|--------------------+
  17 | c             1      |              0               |          2
  18 | c          k1,s1,phi1              k0,s0,phi0                  k2,s2,phi2
  19 | c
  20 | 
  21 |       IMPLICIT NONE
  22 | 
  23 |       Double precision k1,k2,a,ecg,phicg,phi1,phi2,der2,pi
  24 |       Double precision B0,B1,B2,dB1
  25 | 
  26 |       pi=4.D+0*datan(1.D+0)
  27 | 
  28 |       B0=k2*(phi2-phicg)/(2.*ecg)
  29 |       B1=k1*(phi1-phicg)/(2.*ecg)
  30 |       B2=1./phicg-(a-1.)**2/a
  31 |       dB1=k1*k2*(2.*phicg-phi1-phi2)/(2.*ecg**2)
  32 | 
  33 |       der2=B0*B1*B2**2-2.*B0*B1/phicg**2+B2*dB1
  34 | 
  35 |       return
  36 | 
  37 |       END