fix Cp calculation to use sourceside temp, add MyCp output var for de…

…bugging purpose
NREL · Aug 4, 2023 · 24963b4 · 24963b4 · nrel-bot-2b · Aug 4, 2023
1 parent 62ae824
commit 24963b4
Show file tree

Hide file tree

Showing 2 changed files with 10 additions and 1 deletion.
diff --git a/src/EnergyPlus/PlantLoopHeatPumpEIR.cc b/src/EnergyPlus/PlantLoopHeatPumpEIR.cc
@@ -354,12 +354,13 @@ void EIRPlantLoopHeatPump::doPhysics(EnergyPlusData &state, Real64 currentLoad)
     if (this->waterSource) {
         CpSrc = FluidProperties::GetSpecificHeatGlycol(state,
                                                        thisSourcePlantLoop.FluidName,
-                                                       state.dataLoopNodes->Node(this->loadSideNodes.inlet).Temp,
+                                                       state.dataLoopNodes->Node(this->sourceSideNodes.inlet).Temp,
                                                        thisSourcePlantLoop.FluidIndex,
                                                        "PLHPEIR::simulate()");
     } else if (this->airSource) {
         CpSrc = Psychrometrics::PsyCpAirFnW(state.dataEnvrn->OutHumRat);
     }
+    this->sourceSideCp = CpSrc;
     Real64 const sourceMCp = this->sourceSideMassFlowRate * CpSrc;
     this->sourceSideOutletTemp = this->calcSourceOutletTemp(this->sourceSideInletTemp, this->sourceSideHeatTransfer / sourceMCp);
 }
@@ -1225,6 +1226,13 @@ void EIRPlantLoopHeatPump::oneTimeInit(EnergyPlusData &state)
                             OutputProcessor::SOVTimeStepType::System,
                             OutputProcessor::SOVStoreType::Average,
                             this->name);
+        SetupOutputVariable(state,
+                            "MyCp",
+                            OutputProcessor::Unit::J_kgK,
+                            this->sourceSideCp,
+                            OutputProcessor::SOVTimeStepType::System,
+                            OutputProcessor::SOVStoreType::Average,
+                            this->name);
 
         // find this component on the plant
         bool thisErrFlag = false;

diff --git a/src/EnergyPlus/PlantLoopHeatPumpEIR.hh b/src/EnergyPlus/PlantLoopHeatPumpEIR.hh
@@ -119,6 +119,7 @@ namespace EIRPlantLoopHeatPumps {
         Real64 sourceSideEnergy = 0.0;
         Real64 powerEnergy = 0.0;
         bool running = false;
+        Real64 sourceSideCp = 0.0;
 
         // topology variables
         PlantLocation loadSidePlantLoc;