DOUBLE PRECISION FUNCTION DMAG(Z)
      COMPLEX*16 Z
      DOUBLE PRECISION DIMAG,DREAL
      DMAG = DABS(DREAL(Z)) + DABS(DIMAG(Z))
      RETURN
      END