add comment

src/standalone/Vegetation/PlantHydraulics.jl

@@ -342,10 +342,11 @@ Computes the water flux given the absolute potential (pressure/(ρg))
 and the conductivity along
 the flow path between these two points.
 
-We currently assuming an arithmetic
-mean for mean K_sat between the two points (Bonan, 2019; Zhu, 2008)
-to take into account the change in K_sat halfway
-between z1 and z2; this is incorrect for compartments of differing sizes.
+We currently assuming an geometric
+mean for mean K_sat between the two points. Previously,
+    we used the arithmetic mean (Bonan, 2019; Zhu, 2008),
+    but then when the soil K was very low, root extraction would
+    continue. This should be modified for compartments of differing sizes.
 """
 function flux(z1, z2, ψ1, ψ2, K1, K2)
     K_eff = K1 * K2 / (K1 + K2)