Merge pull request #95 from tldr-group/feature---change-err-for-negat…

…ive-flux

Feature   change err for negative flux resolves #94

taufactor/taufactor.py

@@ -197,7 +197,7 @@ def check_vertical_flux(self, conv_crit):
         vert_flux[self.conc[:, :-2, 1:-1, 1:-1] == 0] = 0
         vert_flux[self.conc[:, 1:-1, 1:-1, 1:-1] == 0] = 0
         fl = torch.sum(vert_flux, (0, 2, 3))[1:-1]
-        err = (fl.max() - fl.min())*2/(fl.max() + fl.min())
+        err = (fl.max() - fl.min())/(fl.max())
         if fl.min() == 0:
             return 'zero_flux', torch.mean(fl), err
         if err < conv_crit or torch.isnan(err).item():