1 | include(dom.inc)
   2 | 
   3 |       SUBROUTINE cdss(ph,cb,x,cd)
   4 | 
   5 |       IMPLICIT NONE
   6 | 
   7 | !     IN
   8 |       DOM_REAL   :: ph, cb, x
   9 | 
  10 | !     OUT
  11 |       DOM_REAL   :: cd
  12 | 
  13 | !     LOCAL
  14 |       DOM_REAL   :: r, k1,derfc
  15 | 
  16 |       r  = sqrt((ph*cb)/(2.*x))
  17 |       k1 = x / cb
  18 | 
  19 |       cd = 0.5 * derfc(r * (1 - k1)) +                                  &
  20 |      &     0.5 * derfc(r * (1 + k1)) * exp(2.*ph)
  21 | 
  22 |       END SUBROUTINE cdss