Earth system models consist of many components, coupled in time and space. Standard coupling algorithms introduce a numerical error, which one can compute with iterative coupling methods. We use such a method for a coupled model of a single vertical column of the atmosphere, ocean, and sea ice. We find that coupling errors in the atmosphere and at the ice surface can be substantial and that discontinuous physics parameterizations lead to convergence issues of the iteration.