Introduce preprocessor variable PROJECTED_INTERPOLATION to switch between geodesic FE and projection-based FE

This is a temporary solution; we need something smarter eventually.