| ! { dg-do compile } | |
| SUBROUTINE MLIST(MOLsp,PBCx,PBCy,PBCz, X0) | |
| IMPLICIT NONE | |
| INTEGER, PARAMETER :: NM=16384 | |
| INTEGER :: MOLsp, i | |
| REAL :: PBCx, PBCy, PBCz, boxjmp, HALf=1./2. | |
| REAL :: X0(2,-2:NM) | |
| DO i = 1 , MOLsp | |
| boxjmp = PBCx*INT(X0(1,i)+SIGN(HALf,X0(1,i))) | |
| X0(1,i) = X0(1,i) - boxjmp | |
| boxjmp = PBCy*INT(X0(2,i)+SIGN(HALf,X0(2,i))) | |
| X0(2,i) = X0(2,i) - boxjmp | |
| ENDDO | |
| END | |