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 |       IF(ph.gt.0) THEN 
  20 |         cd = 0.5 * derfc(r * (1 - k1)) +                                &
  21 |      &       0.5 * derfc(r * (1 + k1)) * exp(2.*ph)
  22 |       ELSE
  23 |         cd = 1.0
  24 |       ENDIF
  25 | 
  26 |       END SUBROUTINE cdss