We present the technical implementation of a coarse-resolution coupling between an ice sheet model and an ocean model that allows one to simulate ice–ocean interactions at timescales from centuries to millennia. As ice shelf cavities cannot be resolved in the ocean model at coarse resolution, we bridge the gap using an sub-shelf cavity module. It is shown that the framework is computationally efficient, conserves mass and energy, and can produce a stable coupled state under present-day forcing.