add check to insure that all forcing values are non-fill-values inside routing mask

Author: Joe Hamman

rvic/convolution.py

@@ -151,7 +151,7 @@ def convolution_init(config_file):
                         options['CALENDAR'], data_model.secs_per_step)
     time_handle.end = data_model.end
 
-    data_model.start(time_handle.timestamp)
+    data_model.start(time_handle.timestamp, rout_var)
     # ---------------------------------------------------------------- #
 
     # ---------------------------------------------------------------- #
@@ -353,4 +353,3 @@ def convolution_final(time_handle, hist_tapes):
     # ---------------------------------------------------------------- #
     return
 # -------------------------------------------------------------------- #
-

rvic/core/read_forcing.py

@@ -158,7 +158,7 @@ def __init__(self, path, file_str, time_fld, lat_fld, liq_flds,
     # ---------------------------------------------------------------- #
 
     # ---------------------------------------------------------------- #
-    def start(self, timestamp):
+    def start(self, timestamp, rout_var):
         """ Initialize the first files inputs"""
         # ------------------------------------------------------------ #
         # find and open first file
@@ -205,6 +205,26 @@ def start(self, timestamp):
             else:
                 raise ValueError('unknown forcing units')
         # ------------------------------------------------------------ #
+
+        # ------------------------------------------------------------ #
+        # Compare rout_var mask to forcing mask
+        for fld in self.liq_flds:
+            self.current_fhdl.variables[fld].set_auto_maskandscale(False)
+            try:
+                fill_val = getattr(self.current_fhdl.variables[fld],
+                                   '_FillValue')
+                fmask = np.squeeze(self.current_fhdl.variables[fld][0]) == fill_val
+
+                if np.any(fmask[rout_var.source_y_ind, rout_var.source_x_ind]):
+                    log.error('There are fill values in the routing domain')
+                    log.error('Can not continue...')
+                    raise ValueError('Exiting now due to fill values being '
+                                     'inside the domain defined by the RVIC '
+                                     'parameter file')
+
+            except AttributeError:
+                pass
+        # ------------------------------------------------------------ #
     # ---------------------------------------------------------------- #
 
     # ---------------------------------------------------------------- #