<SYSTEM_TASK:>
Solve the following problem using Python, implementing the functions described below, one line at a time
<END_TASK>
<USER_TASK:>
Description:
 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