YelmoX Bipolar

This flavour is intended to communicate two yelmo domains using an ocean box model in between, i.e., reproducing the bipolar seesaw with a simple and computationally efficient approach.

Some notes: - An additional type is built, yelmox. It contains everything that an individual yelmox domain needs to work. - The coding principle is to reproduce 1:1 yelmox.f90 in those groups in order to be completely analogue to 1-domain yelmox executables.

Dependencies

yelmox_bipolar.f90 mainly needs snapclim.f90 and obm.f90.

Structure

This script is organized in "groups" in order to make the code more clear. Each group is a subroutine that contains the specified lines of the original yelmox.f90 code. Some exceptions are the groups that involve actions related with the coupling between domains and/or the ocean box model.

Groups description:

1. Parameters reading.
2. Output definitions.
3. Initialization: ice sheet model, external models, boundary conditions, output files.
4. Loop: spinup procedure, ice sheet, external models, output.