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