3D anisotropic Finite-Difference modelling on the standard staggered grid

3D anisotropic modelling is becoming more important to resolve more complex structures.
Therefore, I implemented a new algorithms into an already existing program to be able to
simulate anisotropy. This implementation models with finite differences on the standard stag-
gered grid. For complex anisotropic structures a wavefield interpolation is necessary. I have
implemented to interpolation options and compared them. There are the common sinc-function
interpolation and the FD-consistent interpolation which was not used in this context, yet.
I look at the stability of the algorithm and compare different modelling methods. Most results
are as expected and indicate the algorithm works. The FD-consistent interpolation seems to
be a lot more consistent than the sinc-function interpolation and it requires a lower spatial FD
order to do so.