diff --git a/uxarray/grid/coordinates.py b/uxarray/grid/coordinates.py index bb408fdec..2ace69d97 100644 --- a/uxarray/grid/coordinates.py +++ b/uxarray/grid/coordinates.py @@ -121,10 +121,10 @@ def _xyz_to_lonlat_rad( if normalize: x, y, z = _normalize_xyz(x, y, z) - denom = np.abs(x * x + y * y + z * z) - x /= denom - y /= denom - z /= denom + denom = np.hypot(np.hypot(x, y), z) + x = x / denom + y = y / denom + z = z / denom lon = np.arctan2(y, x) lat = np.arcsin(z)