Fix a typo which introduced a bug in PR #1522.

Found through use of some new integration tests in `sopipetest`,
PR will be opened soon with those.

Author: tskisner

sotodlib/toast/ops/load_context.py

@@ -1,6 +1,7 @@
 # Copyright (c) 2022-2024 Simons Observatory.
 # Full license can be found in the top level "LICENSE" file.
 
+import copy
 import re
 import sqlite3
 import yaml
@@ -596,7 +597,7 @@ def _exec(self, data, detectors=None, **kwargs):
 
             # If we have an observation with a specific wafer, modify the detector
             # selection dictionary to include only that wafer.
-            obs_dets_select = dets_select
+            obs_dets_select = copy.deepcopy(dets_select)
             if obs_props[obindx]["wafer"] != "all":
                 if obs_dets_select is None:
                     obs_dets_select = dict()
@@ -722,6 +723,8 @@ def _load_metadata(self, obs_name, session_name, gcomm, dets_select, preproc_con
         # hence one reader).
         if rank == 0:
             # Load metadata
+            msg = f"LoadContext {obs_name} metadata using dets_select={dets_select}"
+            log.debug(msg)
             ctx = open_context(context=self.context, context_file=self.context_file)
             meta = ctx.get_meta(session_name, dets=dets_select)
             if self.context_file is not None:
@@ -755,7 +758,11 @@ def _load_metadata(self, obs_name, session_name, gcomm, dets_select, preproc_con
                     band = fp_cols[
                         f"det_info{self.ax_pathsep}wafer{self.ax_pathsep}bandpass"
                     ].data
-                freq = [float(b[1:]) for b in band]
+
+                fpat = re.compile(r"^f[\d]+$")
+                freq = [
+                    float(b[1:]) if fpat.match(b) is not None else 0.0 for b in band
+                ]
                 bandcenter = np.array(freq) * u.GHz
                 bandwidth = bandcenter * self.bandwidth
                 fp_cols["bandcenter"] = Column(name="bandcenter", data=bandcenter)

sotodlib/toast/ops/load_context_utils.py

@@ -636,14 +636,14 @@ def _process_flag(mask, inbuf, detdata, det_begin, det_end, invert):
                 # Just assign
                 for idet in range(n_recv_det):
                     if daq_units:
-                        nan_mask = np.isnan(recv_data[idet])
+                        nan_mask = np.isnan(recv_2d[idet])
                         if np.count_nonzero(nan_mask) > 0:
                             det_flags[idet] = 1
-                        obs.detdata[field][idet + recv_dets[0], :] = recv_data[
+                        obs.detdata[field][idet + recv_dets[0], :] = recv_2d[
                             idet
                         ].astype(np.int32, casting="unsafe")
                     else:
-                        obs.detdata[field][idet + recv_dets[0], :] = recv_data[idet]
+                        obs.detdata[field][idet + recv_dets[0], :] = recv_2d[idet]
             # Update per-detector flags
             dflags = {
                 obs.local_detectors[recv_dets[0] + x]: defaults.det_mask_invalid