Added error message when input redshifts go outside the grid

Author: cschreib

src/fast++-read_input.cpp

@@ -534,6 +534,18 @@ bool read_fluxes(const options_t& opts, input_state_t& state) {
     if (col_zspec == npos) {
         // If no zspec column, give no zspec to all galaxies
         state.zspec = replicate(fnan, state.id.size());
+    } else {
+        // Check that zspecs are covered by the redshift grid
+        if (min(state.zspec) < opts.z_min) {
+            error("the smallest z_spec is outside of the grid ", min(state.zspec),
+                " vs. ", opts.z_min, ")");
+            return false;
+        }
+        if (max(state.zspec) > opts.z_max) {
+            error("the largest z_spec is outside of the grid ", max(state.zspec),
+                " vs. ", opts.z_max, ")");
+            return false;
+        }
     }
 
     // Convert photometry from [catalog unit] to [uJy]
@@ -1014,6 +1026,30 @@ bool read_photoz(const options_t& opts, input_state_t& state) {
         return false;
     }
 
+    // Check that zspecs are covered by the redshift grid
+    if (min(state.zspec) < opts.z_min) {
+        error("the smallest z_spec is outside of the grid ", min(state.zspec),
+            " vs. ", opts.z_min, ")");
+        return false;
+    }
+    if (max(state.zspec) > opts.z_max) {
+        error("the largest z_spec is outside of the grid ", max(state.zspec),
+            " vs. ", opts.z_max, ")");
+        return false;
+    }
+
+    // Check that zphots are covered by the redshift grid
+    if (min(state.zphot(_,0)) < opts.z_min) {
+        error("the smallest z_phot is outside of the grid ", min(state.zphot(_,0)),
+            " vs. ", opts.z_min, ")");
+        return false;
+    }
+    if (max(state.zphot(_,0)) > opts.z_max) {
+        error("the largest z_phot is outside of the grid ", max(state.zphot(_,0)),
+            " vs. ", opts.z_max, ")");
+        return false;
+    }
+
     return true;
 }