1 | include(dom.inc)
   2 | 
   3 |       SUBROUTINE gicd(s2,s3,x,cd)
   4 | 
   5 |         IMPLICIT NONE
   6 | 
   7 | !       IN
   8 |         DOM_REAL :: s2, s3, x
   9 | 
  10 | !       OUT
  11 |         DOM_REAL :: cd
  12 | 
  13 | !       LOCAL
  14 |         DOM_REAL :: g1, g2, g3, g4, g5, g6
  15 | 
  16 |         g1 = x / s3
  17 |         g2 = SQRT(s2/g1)
  18 |         g3 = -g2 * (1. - g1)
  19 |         g4 = -g2 * (1. + g1)
  20 | 
  21 |         CALL stcd(g3,g5)
  22 |         CALL stcd(g4,g6)
  23 | 
  24 |         cd = g5 + exp(2.*s2) * g6
  25 | 
  26 |       END SUBROUTINE gicd