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