Disable seasonable variance normalization when variance is 0

metno · Apr 25, 2017 · bfde0da · bfde0da
1 parent 28caa1d
commit bfde0da
Showing 1 changed file with 6 additions and 3 deletions.
diff --git a/wxgen/plot.py b/wxgen/plot.py
@@ -449,9 +449,12 @@ def compute_sim_variance(self, array, scales):
          The timeseries is scaled up again to match the average anomaly variance in the timeseries.
          """
          std = np.nanstd(array, axis=1)
-         meanstd = np.nanmean(std)
-         for i in range(0, N):
-            values[:, i] = values[:, i] / std * meanstd
+         if np.min(std) == 0:
+            wxgen.util.warning("Standard deviation of 0 at one or more days. Not normalizing variance")
+         else:
+            meanstd = np.nanmean(std)
+            for i in range(0, N):
+               values[:, i] = values[:, i] / std * meanstd
 
       """
       Compute the variance at different time scales. This is done using a convolution with