polaris.ocean.init_state.layer_thickness_from_geom_interfaces

polaris.ocean.init_state.layer_thickness_from_geom_interfaces(ds)

Compute restingThickness and layerThickness from geometric interface heights in ds.

Both fields equal the geometric layer thickness (appropriate for a quiescent initialisation). Layers below the seafloor are zeroed out using cellMask.

Parameters:

ds (xarray.Dataset) – Dataset containing GeomZInterface and cellMask, for example produced by polaris.ocean.vertical.pstar_init.PStarInitStep.

Returns:

xarray.Dataset – Dataset with restingThickness and layerThickness added.