SCRIP reader does not respect units (in particular, radians) #1564

import uxarray as ux
import numpy as np
import xarray as xr
import os

# Define a tiny 2-cell grid in degrees
center_lon_deg = np.array([10.0, 20.0])
center_lat_deg = np.array([30.0, 40.0])
corner_lon_deg = np.array([[5., 15., 15., 5.], [15., 25., 25., 15.]])
corner_lat_deg = np.array([[25., 25., 35., 35.], [35., 35., 45., 45.]])

# Build a SCRIP dataset with coordinates stored in rad
ds = xr.Dataset()
ds["grid_center_lon"] = xr.DataArray(np.deg2rad(center_lon_deg), dims=["grid_size"], attrs={"units": "radians"})
ds["grid_center_lat"] = xr.DataArray(np.deg2rad(center_lat_deg), dims=["grid_size"], attrs={"units": "radians"})
ds["grid_corner_lon"] = xr.DataArray(np.deg2rad(corner_lon_deg), dims=["grid_size", "grid_corners"], attrs={"units": "radians"})
ds["grid_corner_lat"] = xr.DataArray(np.deg2rad(corner_lat_deg), dims=["grid_size", "grid_corners"], attrs={"units": "radians"})
ds["grid_imask"] = xr.DataArray(np.ones(2, dtype=np.int32), dims=["grid_size"])
ds["grid_area"] = xr.DataArray(np.ones(2), dims=["grid_size"])
ds["grid_dims"] = xr.DataArray(np.array([2], dtype=np.int32), dims=["grid_rank"])

# Write to a small test file to mimic reading a SCRIP nc
out_path = os.path.join(os.path.dirname(os.path.abspath(__file__)), "scrip_radians.nc")
ds.to_netcdf(out_path)

grid = ux.open_grid(out_path)

print("face_lon (got):     ", grid.face_lon.values)
print("face_lon (expected):", center_lon_deg)
# Bug: without the fix, face_lon will be ~[0.17, 0.35] (radian values treated as degrees)
# With fix:             face_lon will be ~[10.0, 20.0]