double precision FUNCTION ARCOS(X)
      implicit none
      double precision x
      include 'COMMON.GEO'
      IF (DABS(X).LT.1.0D0) GOTO 1
      ARCOS=PIPOL*(1.0d0-DSIGN(1.0D0,X))
      RETURN
    1 ARCOS=DACOS(X)
      RETURN
      END