In order to predict how much the sea level will rise in the future, it is crucial to accurately model ice sheets and their heat exchange with the surrounding ocean. In this project, we use the software Elmer/Ice as well as FEniCS to simulate ice sheets by solving the p-Stokes equations. We also model the ice-ocean interactions and the resulting ocean circulation in Greenland Fjords by solving the Boussinesq equations. These problems are both non-linear and involve moving interfaces. We develop numerical methods to run efficient and stable simulations and conduct sets of experiments to investigate future changes of marine glacial melt in a warming climate. In the previous Medium project NAISS 2024/5-210 we did initial simulations of the ocean using FEniCS in 2D and some first trial runs in 3D. We now plan to run multiple sets of model experiments in 2D and in 3D, which will lead to a higher number of simulations and be more computationally expensive. Furthermore, we simulated grounding line migration and ice shelves in 2D and 3D, with a high resolution with Elmer/Ice. These simulations are particularly expensive as a high spatial resolution and the full p-Stokes equations (rather than ice-flow approximations) are needed due to complicated ice-flow physics at the grounding line.