def translate_3D_coordinates_along_axes(x, y, z, shift_x, shift_y, shift_z):
    """
    Translate coordinates along one or more axes.

    Translate or shift coordinates along one or more axes.

    Parameters
    ----------
    x: Numpy array of numbers
        A 1D array representing x axis coordinates. Must have the same length as the `y` and `z` arrays.

    y: Numpy array of numbers
        A 1D array representing y axis coordinates. Must have the same length as the `x` and `z` arrays.

    z: Numpy array of numbers
        A 1D array, representing z axis coordinates. Must have the same length as the `x` and `y` arrays.

    shift_x: number
        A single number, representing the shift along the x axis.

    shift_y: number
        A single number, representing the shift along the y axis.

    shift_z: number
        A single number, representing the shift along the z axis.

    Returns
    -------
    x_shifted: Numpy array of numbers
        The shifted x coordinates.

    y_shifted: Numpy array of numbers
        The shifted y coordinates.

    z_shifted: Numpy array of numbers
        The shifted z coordinates.

    Examples
    --------
    >>> import brainload.spatial as st; import numpy as np
    >>> x = np.array([5, 6])
    >>> y = np.array([7, 8])
    >>> z = np.array([9, 10])
    >>> xt, yt, zt = st.translate_3D_coordinates_along_axes(x, y, z, 2, -4, 0)
    >>> print "%d %d %d" % (xt[0], yt[0], zt[0])     # 7 3 9
    >>> print "%d %d %d" % (xt[1], yt[1], zt[1])     # 8 4 10
    """
    x_shifted = x + shift_x
    y_shifted = y + shift_y
    z_shifted = z + shift_z
    return x_shifted, y_shifted, z_shifted
