compass.ocean.tests.isomip_plus.geom.interpolate_geom

compass.ocean.tests.isomip_plus.geom.interpolate_geom(ds_mesh, ds_geom, min_ocean_fraction, thin_film_present)

Interpolate the ice geometry from the original BISICLES grid to the MPAS mesh.

Parameters:

• ds_mesh (xarray.Dataset) – An MPAS-Ocean mesh

• ds_geom (xarray.Dataset) – Ice-sheet topography produced by compass.ocean.tests.isomip_plus.process_geom.ProcessGeom

• min_ocean_fraction (float) – The minimum ocean fraction after interpolation, below which the cell is masked as land (which is not distinguished from grounded ice)

• thin_film_present (bool) – Whether domain contains a thin film below grounded ice

Returns:

ds_out (xarray.Dataset) – A dataset containing :

• bottomDepthObserved – the bedrock elevation (positive up)

• ssh – the sea surface height

• oceanFracObserved – the fraction of the mesh cell that is ocean

• landIceFraction – the fraction of the mesh cell that is covered by land ice, grounded or floating

• landIceFloatingFraction – the fraction of the mesh cell that is covered by an ice shelf

• smoothedDraftMask – a smoothed version of the floating mask that may be useful for determining where to alter the vertical coordinate to accommodate ice-shelf cavities