@staticmethod
def section_shape(positions, field):
    """defines shape of bilinear interpolated section/surface"""
    if field.ndim == 2:
        return len(positions),
    return len(positions), field.shape[2]
