1 |       subroutine d1fds2(k2,a,ecg,phicg,phi2,der1)
   2 | 
   3 | c     Calcule d[f(a)]/ds2*1/f(a)
   4 | c
   5 | c     Entrées : k2 : coefficient d'absorption moyen dans la maille 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               phi2 : phi dans la maille 2
  10 | c
  11 | c      Sortie : der1 : dérivée première de f(a) par rapport à s1, prise en s1=s2=0, divisée
  12 | c                      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 k2,a,ecg,phicg,phi2,der1,pi
  24 |       Double precision t1,t2
  25 | 
  26 |       pi=4.D+0*datan(1.D+0)
  27 | 
  28 |       t1=k2*(phi2-phicg)/(2.*ecg)
  29 |       t2=1./phicg-(a-1.)**2/a
  30 | 
  31 |       der1=t1*t2
  32 | 
  33 |       return
  34 | 
  35 |       END