docs: def zoom2D(xi, yi, zi, xi_zoom=3., yi_zoom=3., order=3, mode="nearest", cval=0.):
      """Zoom a 2D array, with axes.
  
      Parameters
      ----------
      xi : 1D array
          x axis points.
      yi : 1D array
          y axis points.
      zi : 2D array
          array values. Shape of (x, y).
      xi_zoom : float (optional)
          Zoom factor along x axis. Default is 3.
      yi_zoom : float (optional)
          Zoom factor along y axis. Default is 3.
      order : int (optional)
          The order of the spline interpolation, between 0 and 5. Default is 3.
      mode : {'constant', 'nearest', 'reflect', or 'wrap'}
          Points outside the boundaries of the input are filled according to the
          given mode. Default is nearest.
      cval : scalar (optional)
          Value used for constant mode. Default is 0.0.
      """
      xi = ndimage.interpolation.zoom(xi, xi_zoom, order=order, mode="nearest")
      yi = ndimage.interpolation.zoom(yi, yi_zoom, order=order, mode="nearest")
      zi = ndimage.interpolation.zoom(zi, (xi_zoom, yi_zoom), order=order, mode=mode, cval=cval)
      return xi, yi, zi