Support for Boundaries
of automatically extends arrays to have borders
Borders seamlessly participate in computation
Compare boundary code from
/* Periodic Continuation */ [e of I] wrap U,Uold,V,Vold,P,Pold; [s of I] wrap U,Uold,V,Vold,P,Pold; [se of I] wrap U,Uold,V,Vold,P,Pold;
C C PERIODIC CONTINUATION C uold(m + 1,:n) = uold(1,:n) vold(m + 1,:n) = vold(1,:n) pold(m + 1,:n) = pold(1,:n) u(m + 1,:n) = u(1,:n) v(m + 1,:n) = v(1,:n) p(m + 1,:n) = p(1,:n) uold(:m,n + 1) = uold(:m,1) vold(:m,n + 1) = vold(:m,1) pold(:m,n + 1) = pold(:m,1) u(:m,n + 1) = u(:m,1) v(:m,n + 1) = v(:m,1) p(:m,n + 1) = p(:m,1) uold(m + 1,n + 1) = uold(1,1) vold(m + 1,n + 1) = vold(1,1) pold(m + 1,n + 1) = pold(1,1) u(m + 1,n + 1) = u(1,1) v(m + 1,n + 1) = v(1,1) p(m + 1,n + 1) = p(1,1)