def METHOD_NAME(self):
    level = 1
    scheme = 'nearest'
    result = extract_levels(self.cube, level, scheme)
    expected = np.array([[[4., 5.], [6., 7.]], [[16., 17.], [18., 19.]]])
    self.assert_array_equal(result.data, expected)
    del self.shape[self.z_dim]
    self.assertEqual(result.shape, tuple(self.shape))