compass.landice.tests.ismip7_calibration.ais.run_state.RunState.runtime_setup

RunState.runtime_setup()

Set the number of ocean layers from the forcing file

MPAS takes dimension sizes from the input stream only, and the mesh file has no nISMIP6OceanLayers. Without this the 3-D forcing fields are allocated against a zero-length dimension and the run dies trying to allocate hundreds of GB. It is read from the forcing rather than hard-coded, so it stays right if ISMIP7 changes the vertical grid.