docs: def xyz2lonlat(x__, y__, z__):
      """Get longitudes from cartesian coordinates.
      """
      R = 6370997.0
      lons = da.rad2deg(da.arccos(x__ / da.sqrt(x__ ** 2 + y__ ** 2))) * da.sign(y__)
      lats = da.sign(z__) * (90 - da.rad2deg(da.arcsin(da.sqrt(x__ ** 2 + y__ ** 2) / R)))
  
      return lons, lats