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