Fixed severed bug in how the ICON input data was thresholded

cctrigrid/prepare.py

@@ -49,8 +49,8 @@ def field_icon(path, file, var, threshold, cubulation):
     
     # field on triangular grid
     field = xr.open_dataset(path + file)[var].squeeze().values
-    field[field>=threshold] = 1
-    field[field< threshold] = 0
+    # threshold field
+    field = np.where(field<threshold, 0, 1)
     
     # determine size of array that holds cubulation
     array_size=0
@@ -75,8 +75,8 @@ def field_icon_lev(path, file, var, threshold, cubulation):
     
     # field on triangular grid
     field = xr.open_dataset(path + file)[var].squeeze().values
-    field[field>=threshold] = 1
-    field[field< threshold] = 0
+    # threshold field
+    field = np.where(field<threshold, 0, 1)
     
     nlev = np.shape(field)[0]