1! PR middle-end/77719 2! { dg-do compile } 3! { dg-options "-O3 -ffast-math" } 4 5SUBROUTINE urep_egr(erep,derep,surr) 6 INTEGER, PARAMETER :: dp=8 7 REAL(dp), INTENT(inout) :: erep, derep(3) 8 REAL(dp), INTENT(in) :: surr(2) 9 REAL(dp) :: de_z, rz 10 INTEGER :: isp,spdim,jsp,nsp 11 IF (n_urpoly > 0) THEN 12 IF (r < spxr(1,1)) THEN 13 ispg: DO isp = 1,spdim ! condition ca) 14 IF (isp /= spdim) THEN 15 nsp = 5 ! condition cb 16 DO jsp = 0,nsp 17 IF( jsp <= 3 ) THEN 18 ELSE 19 erep = erep + surr(jsp-3)*rz**(jsp) 20 ENDIF 21 END DO 22 END IF 23 END DO ispg 24 END IF 25 END IF 26END SUBROUTINE urep_egr 27