Update leakage.c

Adjustment to reported pipe leakage when one node is fixed grade was corrected.
This commit is contained in:
Lew Rossman
2025-06-27 09:44:32 -04:00
parent d9eba03a9b
commit e9cf8d7d3f

View File

@@ -307,10 +307,10 @@ double findlinkleakage(Project *pr, int i)
hsqrt = sqrt(h2); hsqrt = sqrt(h2);
q2 = c * (a + m * h2) * hsqrt; q2 = c * (a + m * h2) * hsqrt;
} }
// Adjust leakage flows to account for one node being fixed grade // Adjust leakage flows to account for one node being fixed grade
if (q2 == 0.0) q1 *= 2.0; if (n2 > net->Njuncs) q1 *= 2.0;
if (q1 == 0.0) q2 *= 2.0; if (n1 > net->Njuncs) q2 *= 2.0;
return q1 + q2; return q1 + q2;
} }