1! { dg-do compile }
2! { dg-additional-options "-std=legacy" }
3! { dg-require-effective-target vect_double }
4      SUBROUTINE CALC2
5      IMPLICIT REAL*8	(A-H, O-Z)
6      PARAMETER (N1=1335, N2=1335)
7
8      COMMON  U(N1,N2), V(N1,N2), P(N1,N2),
9     *        UNEW(N1,N2), VNEW(N1,N2),
10     1        PNEW(N1,N2), UOLD(N1,N2),
11     *        VOLD(N1,N2), POLD(N1,N2),
12     2        CU(N1,N2), CV(N1,N2),
13     *        Z(N1,N2), H(N1,N2), PSI(N1,N2)
14      COMMON /CONS/ DT,TDT,DX,DY,A,ALPHA,ITMAX,MPRINT,M,N,MP1,
15     1              NP1,EL,PI,TPI,DI,DJ,PCF
16      TDTS8 = TDT/8.D0
17      TDTSDX = TDT/DX
18      TDTSDY = TDT/DY
19
20      DO 200 J=1,N
21      DO 200 I=1,M
22      UNEW(I+1,J) = UOLD(I+1,J)+
23     1    TDTS8*(Z(I+1,J+1)+Z(I+1,J))*(CV(I+1,J+1)+CV(I,J+1)+CV(I,J)
24     2       +CV(I+1,J))-TDTSDX*(H(I+1,J)-H(I,J))
25      VNEW(I,J+1) = VOLD(I,J+1)-TDTS8*(Z(I+1,J+1)+Z(I,J+1))
26     1       *(CU(I+1,J+1)+CU(I,J+1)+CU(I,J)+CU(I+1,J))
27     2       -TDTSDY*(H(I,J+1)-H(I,J))
28      PNEW(I,J) = POLD(I,J)-TDTSDX*(CU(I+1,J)-CU(I,J))
29     1       -TDTSDY*(CV(I,J+1)-CV(I,J))
30  200 CONTINUE
31      RETURN
32      END
33! { dg-final { scan-tree-dump-times "LOOP VECTORIZED" 1 "vect" } }
34