1! { dg-additional-options "-Ofast" }
2MODULE module_ra_gfdleta
3      INTEGER, PARAMETER              :: NBLY=15
4      REAL   , SAVE :: EM1(28,180),EM1WDE(28,180),TABLE1(28,180),     &
5                           TABLE2(28,180),TABLE3(28,180),EM3(28,180), &
6                           SOURCE(28,NBLY), DSRCE(28,NBLY)
7CONTAINS
8      SUBROUTINE TABLE
9 INTEGER, PARAMETER :: NBLX=47
10 INTEGER , PARAMETER:: NBLW = 163
11      REAL ::  &
12               SUM(28,180),PERTSM(28,180),SUM3(28,180),       &
13               SUMWDE(28,180),SRCWD(28,NBLX),SRC1NB(28,NBLW), &
14               DBDTNB(28,NBLW)
15      REAL ::  &
16               ZMASS(181),ZROOT(181),SC(28),DSC(28),XTEMV(28), &
17               TFOUR(28),FORTCU(28),X(28),X1(28),X2(180),SRCS(28), &
18               R1T(28),R2(28),S2(28),T3(28),R1WD(28)
19      REAL ::  EXPO(180),FAC(180)
20      I = 0
21      DO 417 J=121,180
22      FAC(J)=ZMASS(J)*(ONE-(ONE+X2(J))*EXPO(J))/(X2(J)*X2(J))
23417   CONTINUE
24      DO 421 J=121,180
25      SUM3(I,J)=SUM3(I,J)+DBDTNB(I,N)*FAC(J)
26421   CONTINUE
27      IF (CENT.GT.160. .AND. CENT.LT.560.) THEN
28         DO 420 J=1,180
29         DO 420 I=1,28
30         SUMWDE(I,J)=SUMWDE(I,J)+SRC1NB(I,N)*EXPO(J)
31420      CONTINUE
32      ENDIF
33      DO 433 J=121,180
34      EM3(I,J)=SUM3(I,J)/FORTCU(I)
35433   CONTINUE
36      DO 501 I=1,28
37      EM1WDE(I,J)=SUMWDE(I,J)/TFOUR(I)
38501   CONTINUE
39      END SUBROUTINE TABLE
40      END MODULE module_RA_GFDLETA
41