fixed logic for xesmf

Author: kthyng

extract_model/extract_model.py

@@ -393,6 +393,10 @@ def select(
             )
             kwargs_out["weights"] = weights
             kwargs_out["interp_coords"] = {}
+        elif not XESMF_AVAILABLE and horizontal_interp_code == "xesmf":
+            raise ModuleNotFoundError(
+                "xESMF is not available so horizontal interpolation in 2D cannot be performed."
+            )
         elif horizontal_interp_code == "delaunay":
 
             # This might be faster for interpolation between a few points and many points
@@ -542,11 +546,8 @@ def pt_in_itriangle_proj(ix, iy):
             # interp_coords are the coords and indices that went into the interpolation
             da, interp_coords = interp_with_barycentric(da, ixs, iys, lam)
             kwargs_out["interp_coords"] = interp_coords
-
         else:
-            raise ModuleNotFoundError(
-                "xESMF is not available so horizontal interpolation in 2D cannot be performed."
-            )
+            raise RuntimeError("Horizontal interpolation cannot be run.")
 
         end_time = time()
         print("time: ", end_time - start_time)