Initial (probably buggy) implementation of the 3D d_to_a grid interpolation.
Committing latest m_operators.f90