find drifts

src/glide/cli.py

@@ -228,7 +228,9 @@ def l2(
 
     merged = process_l1.calculate_thermodynamics(merged, conf)
 
-    out = profiles.get_profiles(merged, shallowest_profile, min_surface_time)
+    out = profiles.get_profiles(
+        merged, shallowest_profile, min_surface_time, flt=flt_raw
+    )
 
     out = profiles.assign_surface_state(out, flt=flt_raw)
 

src/glide/process_l2.py

@@ -231,6 +231,9 @@ def bin_l2(
 
     for row in profile_indexes:
         profile = ds.isel(i=slice(row[0], row[1]))
+        if not np.isfinite(profile.depth.values).any():
+            _log.debug("Skipping profile [%d:%d] with no finite depth", row[0], row[1])
+            continue
         time = profile.time.values
         idx_mid = np.searchsorted(time, 0.5 * (time[0] + time[-1]))