Since receiving my PhD from the University of Toronto in 2000. I joined Prof. Jerry Mitrovica's group in Canada. After he left for Harvard, the collaboration continued and I became his Harvard group member in October 2017. The focus of my research is on development and application of numerical tools to examine problems related to the ice age geodynamics: visco-elastic deformation, sea level change, tidal response, Earth rotation, stress and gravity field calculation. Generally, we are trying to understand how Earth (or similarly structured planets) would respond to surface and/or potential forcing, given a loading history and an assumed visco-elastic structure, particularly three-dimensional, heterogeneous viscosity, as may be inferred from e.g. seismic tomography. Other important 3-D effects which we include in the models are variations in the lithospheric thickness, plate boundaries, slabs, low viscosity wedges, etc. This list is open-ended. The results are ultimately compared to observables, such as sea level markers or present-day deformation rates, available from satellite measuremets. We can, for example, provide a correction for the glacial isostatic adjustment contribution or assess the impact of a 3-D structure on such predictions, since conventional, radially stratified models are handled in a much more efficient way. The 3-D models require an excessive computer power, even for the forward problem. Currently, the simulations are performed on the Odyssey cluster (Harvard) with a finite volume MPI code, developed in Toronto in the early 2000-s. The latter is an ongoing project, including coding, maintenance, consulting and working with interested researchers on improvements. As a side product, this development has stimulated a strong interest in interpolation techniques and adaprive grid generation.