code: def slice_orthogonal(dataset, x=None, y=None, z=None,
                           generate_triangles=False, contour=False):
          
          output = vtki.MultiBlock()
          # Create the three slices
          if x is None:
              x = dataset.center[0]
          if y is None:
              y = dataset.center[1]
          if z is None:
              z = dataset.center[2]
          output[0, 'YZ'] = dataset.slice(normal='x', origin=[x,y,z], generate_triangles=generate_triangles)
          output[1, 'XZ'] = dataset.slice(normal='y', origin=[x,y,z], generate_triangles=generate_triangles)
          output[2, 'XY'] = dataset.slice(normal='z', origin=[x,y,z], generate_triangles=generate_triangles)
          return output
docs: Creates three orthogonal slices through the dataset on the three
          caresian planes. Yields a MutliBlock dataset of the three slices
  
          Parameters
          ----------
          x : float
              The X location of the YZ slice
  
          y : float
              The Y location of the XZ slice
  
          z : float
              The Z location of the XY slice
  
          generate_triangles: bool, optional
              If this is enabled (``False`` by default), the output will be
              triangles otherwise, the output will be the intersection polygons.
  
          contour : bool, optional
              If True, apply a ``contour`` filter after slicing