diff --git a/xcube/core/resampling/spatial.py b/xcube/core/resampling/spatial.py
index 2ce5e2f2a..22b3e79b3 100644
--- a/xcube/core/resampling/spatial.py
+++ b/xcube/core/resampling/spatial.py
@@ -208,34 +208,23 @@ def resample_in_space(
     # If CRSes are not both geographic and their CRSes are different
     # transform the source_gm so its CRS matches the target CRS:
     transformed_source_gm = source_gm.transform(target_gm.crs, xy_res=target_gm.xy_res)
-    # source_ds = source_ds.drop_vars(source_gm.xy_dim_names)
-    # if "crs" in source_ds:
-    #     source_ds = source_ds.drop_vars("crs")
-    # if "spatial_ref" in source_ds:
-    #     source_ds = source_ds.drop_vars("spatial_ref")
-    # source_ds = source_ds.copy()
-    # for var in source_ds.data_vars:
-    #     if "grid_mapping" in source_ds[var].attrs:
-    #         attrs = source_ds[var].attrs
-    #         del attrs["grid_mapping"]
-    #         source_ds[var] = source_ds[var].assign_attrs(attrs)
-    # transformed_x, transformed_y = transformed_source_gm.xy_coords
-    # attrs = dict(grid_mapping="spatial_ref")
-    # transformed_x.attrs = attrs
-    # transformed_y.attrs = attrs
-    # source_ds = source_ds.assign_coords(
-    #     spatial_ref=xr.DataArray(0, attrs=transformed_source_gm.crs.to_cf()),
-    #     transformed_x=transformed_x,
-    #     transformed_y=transformed_y,
-    # )
+    source_ds = source_ds.drop_vars(source_gm.xy_dim_names)
+    if "crs" in source_ds:
+        source_ds = source_ds.drop_vars("crs")
+    if "spatial_ref" in source_ds:
+        source_ds = source_ds.drop_vars("spatial_ref")
+    source_ds = source_ds.copy()
+    for var in source_ds.data_vars:
+        if "grid_mapping" in source_ds[var].attrs:
+            attrs = source_ds[var].attrs
+            del attrs["grid_mapping"]
+            source_ds[var] = source_ds[var].assign_attrs(attrs)
     transformed_x, transformed_y = transformed_source_gm.xy_coords
-    attrs = dict(grid_mapping="transformed_spatial_ref")
+    attrs = dict(grid_mapping="spatial_ref")
     transformed_x.attrs = attrs
     transformed_y.attrs = attrs
     source_ds = source_ds.assign_coords(
-        transformed_spatial_ref=xr.DataArray(
-            0, attrs=transformed_source_gm.crs.to_cf()
-        ),
+        spatial_ref=xr.DataArray(0, attrs=transformed_source_gm.crs.to_cf()),
         transformed_x=transformed_x,
         transformed_y=transformed_y,
     )