function sindg (x)
  c june 1977 edition.   w. fullerton, c3, los alamos scientific lab.
        data raddeg / .017453292519943296e0 /
  c
        sindg = sin (raddeg*x)
  c
        if (amod(x,90.).ne.0.) return
        n = abs(x)/90.0 + 0.5
        n = mod (n, 2)
        if (n.eq.0) sindg = 0.
        if (n.eq.1) sindg = sign (1.0, sindg)
  c
        return
        end