! { dg-do compile }
! { dg-additional-options "-Ofast" }
      SUBROUTINE GAUBON(NV,PTS,PP)
      IMPLICIT DOUBLE PRECISION(A-H,O-Z)
      PARAMETER (MXSP=250)
      DIMENSION PTS(3,10),PP(3)
      COMMON /PCMPLY/ XE(MXSP),YE(MXSP),ZE(MXSP)
      DATA PI/3.141592653589793D+00/
      DATA ZERO/0.0D+00/
      DO I = 1, NV
      PP(1) = PP(1) + (PTS(1,I)-XE(NS))
      PP(2) = PP(2) + (PTS(2,I)-YE(NS))
      PP(3) = PP(3) + (PTS(3,I)-ZE(NS))
      ENDDO
      END