eusoubrasileiro
commented
9 years ago - indexes must be correct by +2/-2 due unused region due 4th order space.
- the source or any indexed position must correct for unused regions [0:2,:], [:-2,:], [:,0:2], [:,-2:]
- the real calculation grid is (nx+2*pad-4, nz+pad-4)