def _parse_bounds_vertex(da, dim="bnds", position=[0, 1]):
    """Convenience function to extract positions from bounds/verticies"""
    return tuple([da.isel({dim: i}).load().data for i in position])