Implement Absorbing Boundary Conditions using PML

[thecodeboss]

To simulate absorbing boundaries, we need to create a small pressure layer at the boundaries. On the boundary of this pressure layer, the forcing terms are multiplied by a reflectivity constant between 0 and 1. See paper for more details.