1! { dg-do compile } 2 3 SUBROUTINE MLIST(MOLsp,PBCx,PBCy,PBCz, X0) 4 IMPLICIT NONE 5 INTEGER, PARAMETER :: NM=16384 6 INTEGER :: MOLsp, i 7 REAL :: PBCx, PBCy, PBCz, boxjmp, HALf=1./2. 8 REAL :: X0(2,-2:NM) 9 10 DO i = 1 , MOLsp 11 boxjmp = PBCx*INT(X0(1,i)+SIGN(HALf,X0(1,i))) 12 X0(1,i) = X0(1,i) - boxjmp 13 boxjmp = PBCy*INT(X0(2,i)+SIGN(HALf,X0(2,i))) 14 X0(2,i) = X0(2,i) - boxjmp 15 ENDDO 16 END 17 18 19