Scalable robust solvers for unstructured FE geodynamic models, with an application to geoid calculations from high resolution tomography models in a subduction setting