diff --git a/src/hydcoeffs.c b/src/hydcoeffs.c index 7b2dc47..8fe5426 100644 --- a/src/hydcoeffs.c +++ b/src/hydcoeffs.c @@ -486,7 +486,7 @@ void demandheadloss(Project *pr, int i, double dp, double n, } // Use power head loss function for demand less than full - else if (r <= 1.0) //(2.2.1) + else if (r <= 1.0) { *hgrad = n * dp * pow(r, n - 1.0) / dfull; // ... use linear function for very small gradient diff --git a/src/hydsolver.c b/src/hydsolver.c index 7f80a29..c3b5979 100644 --- a/src/hydsolver.c +++ b/src/hydsolver.c @@ -554,7 +554,7 @@ void newdemandflows(Project *pr, Hydbalance *hbal, double *qsum, double *dqsum) dq = (hloss - dh) / hgrad; dq *= hyd->RelaxFactor; - // Prevent a flow change greater than full demand //(2.2.1) + // Prevent a flow change greater than full demand if (fabs(dq) > hyd->NodeDemand[i]) dq = 0.5 * SGN(dq) * hyd->NodeDemand[i]; hyd->DemandFlow[i] -= dq;