o Fix errors found in new test_ugrid changes, bug was in _ugrid.py - shouldn't add face/edge coordinates when they don't exist

Author: rajeeja

uxarray/io/_exodus.py

@@ -150,15 +150,15 @@ def _encode_exodus(ds, outfile=None):
     now = datetime.now()
     date = now.strftime("%Y:%m:%d")
     time = now.strftime("%H:%M:%S")
-    fp_word = np.int32(8)  # Assuming INT_DTYPE is int32 based on usage
+    fp_word = INT_DTYPE(8)
     exo_version = np.float32(5.0)
     api_version = np.float32(5.0)
 
     exo_ds.attrs = {
         "api_version": api_version,
         "version": exo_version,
         "floating_point_word_size": fp_word,
-        "file_size": 0,  # Will be 0 for an in-memory representation
+        "file_size": 0,  
     }
 
     if outfile:

uxarray/io/_ugrid.py

@@ -20,18 +20,18 @@ def _read_ugrid(ds):
         node_lon_name: ugrid.NODE_COORDINATES[0],
         node_lat_name: ugrid.NODE_COORDINATES[1],
     }
-
-    if "edge_coordinates" in ds["grid_topology"].attrs:
-        # get the names of edge_lon and edge_lat, if they exist
-        edge_lon_name, edge_lat_name = ds["grid_topology"].edge_coordinates.split()
-        coord_dict[edge_lon_name] = ugrid.EDGE_COORDINATES[0]
-        coord_dict[edge_lat_name] = ugrid.EDGE_COORDINATES[1]
-
-    if "face_coordinates" in ds["grid_topology"].attrs:
-        # get the names of face_lon and face_lat, if they exist
-        face_lon_name, face_lat_name = ds["grid_topology"].face_coordinates.split()
-        coord_dict[face_lon_name] = ugrid.FACE_COORDINATES[0]
-        coord_dict[face_lat_name] = ugrid.FACE_COORDINATES[1]
+    if "edge_lon" in ds and "edge_lat" in ds:
+        if "edge_coordinates" in ds["grid_topology"].attrs:
+            # get the names of edge_lon and edge_lat, if they exist
+            edge_lon_name, edge_lat_name = ds["grid_topology"].edge_coordinates.split()
+            coord_dict[edge_lon_name] = ugrid.EDGE_COORDINATES[0]
+            coord_dict[edge_lat_name] = ugrid.EDGE_COORDINATES[1]
+    if "face_lon" in ds and "face_lat" in ds:
+        if "face_coordinates" in ds["grid_topology"].attrs:
+            # get the names of face_lon and face_lat, if they exist
+            face_lon_name, face_lat_name = ds["grid_topology"].face_coordinates.split()
+            coord_dict[face_lon_name] = ugrid.FACE_COORDINATES[0]
+            coord_dict[face_lat_name] = ugrid.FACE_COORDINATES[1]
 
     ds = ds.rename(coord_dict)