Numerical simulations for the non-linear Molodensky problem