1 | include(dom.inc)
   2 | 
   3 |       SUBROUTINE stcd(g3,g5)
   4 | 
   5 |         IMPLICIT NONE
   6 | 
   7 | !       IN
   8 |         DOM_REAL g3
   9 | 
  10 | !       OUT
  11 |         DOM_REAL g5
  12 | 
  13 | !       LOCAL
  14 |         DOM_REAL s2, s1, erfc
  15 | 
  16 |         s1 = -g3/SQRT(2.)
  17 |         s2 = s1
  18 |         s2 = erfc(s2)
  19 |         g5 = s2/2.
  20 | 
  21 |       END SUBROUTINE stcd