From 15c481204dffa7d418f7063390712e72c3c4f259 Mon Sep 17 00:00:00 2001
From: Tobias Holzmann <>
Date: Mon, 17 Jan 2022 09:39:58 +0000
Subject: [PATCH] ENH: comfort function object - added operative temperature.
 See #2325

---
 src/functionObjects/field/comfort/comfort.C | 10 ++++++++--
 1 file changed, 8 insertions(+), 2 deletions(-)

diff --git a/src/functionObjects/field/comfort/comfort.C b/src/functionObjects/field/comfort/comfort.C
index f073ea372d..56423dc158 100644
--- a/src/functionObjects/field/comfort/comfort.C
+++ b/src/functionObjects/field/comfort/comfort.C
@@ -455,15 +455,20 @@ bool Foam::functionObjects::comfort::execute()
     tmp<volScalarField> DR =
         correctUnit*(factor12 - T)*pow(Umag - Umin, 0.62)*(pre*Umag*TI + C1);
 
+    // Calculate the operative temperature
+    tmp<volScalarField> Top = (T + Trad)/2;
+
     // Workaround
     word fieldNamePMV = "PMV";
     word fieldNamePPD = "PPD";
     word fieldNameDR = "DR";
+    word fieldNameTop = "Top";
 
     return
         store(fieldNamePMV, PMV)
      && store(fieldNamePPD, PPD)
-     && store(fieldNameDR, DR);
+     && store(fieldNameDR, DR)
+     && store(fieldNameDR, Top);
 }
 
 bool Foam::functionObjects::comfort::write()
@@ -471,7 +476,8 @@ bool Foam::functionObjects::comfort::write()
     return
         writeObject("PMV")
      && writeObject("PPD")
-     && writeObject("DR");
+     && writeObject("DR")
+     && writeObject("Top");
 }