program example13 implicit none integer i, imax parameter (imax = 3) real a(0:imax) c= region :: r1 = pointed(dim=1) do i = 0, imax c= stencil readOnce, atLeast, r1 :: a c= stencil readOnce, atMost, forward(depth=2, dim=1) :: a a(i) = a(i) + a(i+2) end do end