BUG: fix weertman exponent

Author: enigne

Project.toml

@@ -30,8 +30,8 @@ ColorSchemes = "3.25"
 Enzyme = "0.13"
 JLD2 = "0.6"
 JSOSolvers = "0.14.1"
-Lux = "1.21"
-MAT = "0.10"
+Lux = "1.26"
+MAT = "0.11"
 MLUtils = "0.4.8"
 MadNLP = "0.8"
 Makie = "0.24"

src/core/friction.jl

@@ -125,10 +125,10 @@ end #}}}
 	m = GetInputValue(friction.m_input, gauss, i)
 	vmag = VelMag(friction, gauss, i)
 
-	if vmag==0.0 && m<1.0
+	if vmag==0.0 && (1.0/m)<1.0
 		return 0.0
 	else
-		return c^2*vmag^(m-1)
+		return c^2*vmag^(1.0/m-1.0)
 	end
 end#}}}
 @inline function Alpha2(friction::CoreSchoofFriction, gauss::GaussTria, i::Int64) #{{{