docs: def unorm(data, ord=None, axis=None, keepdims=False):
      """Matrix or vector norm that preserves units
  
      This is a wrapper around np.linalg.norm that preserves units. See
      the documentation for that function for descriptions of the keyword
      arguments.
  
      Examples
      --------
      >>> from unyt import km
      >>> data = [1, 2, 3]*km
      >>> print(unorm(data))
      3.7416573867739413 km
      """
      norm = np.linalg.norm(data, ord=ord, axis=axis, keepdims=keepdims)
      if norm.shape == ():
          return unyt_quantity(norm, data.units)
      return unyt_array(norm, data.units)