# HG changeset patch
# User mschaefer
# Date 1528205854 -7200
# Node ID 8675912c7f4881249daa62c285aed537371bb21e
# Parent  2fb72a0aa313e363ac46d0e54f44f9cdd6040be9
Fixed tkh calculation (sign of an exponent of the critical shields term)

diff -r 2fb72a0aa313 -r 8675912c7f48 artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/tkhcalculation/TkhCalculator.java
--- a/artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/tkhcalculation/TkhCalculator.java	Tue Jun 05 11:44:05 2018 +0200
+++ b/artifacts/src/main/java/org/dive4elements/river/artifacts/sinfo/tkhcalculation/TkhCalculator.java	Tue Jun 05 15:37:34 2018 +0200
@@ -209,7 +209,7 @@
 
         final double froude = vm / Math.sqrt(PHYS_G * h);
         final double partReynolds = Math.sqrt((PHYS_SPECGRAV_S - 1) * PHYS_G * d50) / PHYS_VISCOSITY_NUE * d50;
-        final double critShields = 0.22 * Math.pow(partReynolds, -0.6) + 0.06 * Math.pow(10, 7.7 * Math.pow(partReynolds, -0.6));
+        final double critShields = 0.22 * Math.pow(partReynolds, -0.6) + 0.06 * Math.pow(10, -7.7 * Math.pow(partReynolds, -0.6));
         final double critTau = critShields * (PHYS_GRAIN_DENSITY_RHOS - PHYS_WATER_DENSITY_RHO) * PHYS_G * d50;
         final double tkh = 100 * h * (1 - Math.pow(froude, 2)) / (2 * PHYS_VELOCCOEFF_N * PHYS_FORMCOEFF_ALPHA) * (1 - critTau / tau);
         // Some regular input values may give a negative calculation result; that is unwanted