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