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