compass.landice.tests.ismip7_calibration.ais.melt_model.initial_draft

compass.landice.tests.ismip7_calibration.ais.melt_model.initial_draft(mesh_file, constants=None)[source]

Reconstruct the ice draft from the mesh file, before the timestep.

The diagnostic run evolves the geometry. Melt is applied over the single timestep and thins the ice by up to a metre or so, so the lowerSurface and thickness written to the output are post-step while TFdraft was computed pre-step. Any check that pairs output geometry with output melt is therefore inconsistent by up to a metre of draft; the initial draft has to come from the mesh file instead.

Parameters:

  • mesh_file (str) – The MALI mesh file the run started from

  • constants (compass.landice.tests.ismip7_calibration.quadratic.Constants, optional) – The densities to use; defaults to MALI’s

Returns:

draft (xarray.DataArray) – The ice draft, negative below sea level