From 936df25bcb5cd637f9c8d4849ad84d0f6a6cbdf8 Mon Sep 17 00:00:00 2001
From: Chloe Mimeau <Chloe.Mimeau@imag.fr>
Date: Tue, 22 Apr 2014 08:28:03 +0000
Subject: [PATCH] integrate only on the inlet and outlet surfaces of the
control box
---
HySoP/hysop/operator/monitors/compute_forces.py | 17 +++++++++++------
1 file changed, 11 insertions(+), 6 deletions(-)
diff --git a/HySoP/hysop/operator/monitors/compute_forces.py b/HySoP/hysop/operator/monitors/compute_forces.py
index 44694fdb1..6abafacc0 100644
--- a/HySoP/hysop/operator/monitors/compute_forces.py
+++ b/HySoP/hysop/operator/monitors/compute_forces.py
@@ -149,12 +149,17 @@ class DragAndLift(Monitoring):
# Update previous for next time step ...
self._previous[...] = self._buffer[...]
# -- Integrals on surfaces --
- for s in self._voc.upperS:
- self._buffer = self._integrateOnSurface(s, self._buffer)
- self.force += self._buffer
- for s in self._voc.lowerS:
- self._buffer = self._integrateOnSurface(s, self._buffer)
- self.force += self._buffer
+# for s in self._voc.upperS:
+# self._buffer = self._integrateOnSurface(s, self._buffer)
+# self.force += self._buffer
+# for s in self._voc.lowerS:
+# self._buffer = self._integrateOnSurface(s, self._buffer)
+# self.force += self._buffer
+
+ self._buffer = self._integrateOnSurface(self._voc.upperS[0], self._buffer)
+ self.force += self._buffer
+ self._buffer = self._integrateOnSurface(self._voc.lowerS[0], self._buffer)
+ self.force += self._buffer
# Reduce results over all MPI process in topo
self.mpi_sum()
--
GitLab