Assign heat dist is modification fails.

bigladder · Dec 2, 2023 · c89b426 · c89b426
1 parent a405503
commit c89b426
Showing 1 changed file with 7 additions and 1 deletion.
diff --git a/src/HPWH.cc b/src/HPWH.cc
@@ -2762,6 +2762,7 @@ void HPWH::modifyHeatDistribution(std::vector<double> &heatDistribution)
 	std::vector<double> extTankTemps_C(heatDistribution.size());
 	resampleIntensive(extTankTemps_C,tankTemps_C);
 
+	double totDist = 0.;
 	for(std::size_t iNode = 0; iNode < heatDistribution.size(); ++iNode) {
 		double dist = 0.;
 		if(iNode >= lowestNode){
@@ -2773,8 +2774,13 @@ void HPWH::modifyHeatDistribution(std::vector<double> &heatDistribution)
 				dist = 0.;
 		}
 		heatDistribution[iNode] = dist;
+		totDist += dist;
 	}
-	normalize(heatDistribution);
+
+	if(totDist > 0.)
+		normalize(heatDistribution);
+	else
+		heatDistribution.assign(heatDistribution.size(), 1./static_cast<double>(heatDistribution.size()));
 
 	for(auto &heat: heatDistribution)
 		heat *= totalHeat;