! { dg-do compile } ! { dg-additional-options "-Ofast -std=legacy" } MODULE module_ra_gfdleta INTEGER, PARAMETER :: NBLY=15 REAL , SAVE :: EM1(28,180),EM1WDE(28,180),TABLE1(28,180), & TABLE2(28,180),TABLE3(28,180),EM3(28,180), & SOURCE(28,NBLY), DSRCE(28,NBLY) CONTAINS SUBROUTINE TABLE INTEGER, PARAMETER :: NBLX=47 INTEGER , PARAMETER:: NBLW = 163 REAL :: & SUM(28,180),PERTSM(28,180),SUM3(28,180), & SUMWDE(28,180),SRCWD(28,NBLX),SRC1NB(28,NBLW), & DBDTNB(28,NBLW) REAL :: & ZMASS(181),ZROOT(181),SC(28),DSC(28),XTEMV(28), & TFOUR(28),FORTCU(28),X(28),X1(28),X2(180),SRCS(28), & R1T(28),R2(28),S2(28),T3(28),R1WD(28) REAL :: EXPO(180),FAC(180) I = 0 DO 417 J=121,180 FAC(J)=ZMASS(J)*(ONE-(ONE+X2(J))*EXPO(J))/(X2(J)*X2(J)) 417 CONTINUE DO 421 J=121,180 SUM3(I,J)=SUM3(I,J)+DBDTNB(I,N)*FAC(J) 421 CONTINUE IF (CENT.GT.160. .AND. CENT.LT.560.) THEN DO 420 J=1,180 DO 420 I=1,28 SUMWDE(I,J)=SUMWDE(I,J)+SRC1NB(I,N)*EXPO(J) 420 CONTINUE ENDIF DO 433 J=121,180 EM3(I,J)=SUM3(I,J)/FORTCU(I) 433 CONTINUE DO 501 I=1,28 EM1WDE(I,J)=SUMWDE(I,J)/TFOUR(I) 501 CONTINUE END SUBROUTINE TABLE END MODULE module_RA_GFDLETA