From 5b8bad24546c0c01b11113fcd9b212bc633b6757 Mon Sep 17 00:00:00 2001
From: Adrian Kummerlaender
Date: Thu, 17 Jan 2019 12:45:25 +0100
Subject: Fix non-equilibrium scaling factor

---
 src/refinement/coupler2D.hh | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/refinement/coupler2D.hh b/src/refinement/coupler2D.hh
index 2218c28..91a7db8 100644
--- a/src/refinement/coupler2D.hh
+++ b/src/refinement/coupler2D.hh
@@ -33,7 +33,8 @@ namespace olb {
 template <typename T, template<typename> class DESCRIPTOR>
 T Coupler2D<T,DESCRIPTOR>::getScalingFactor() const
 {
-  return 4. - _coarse.getConverter().getLatticeRelaxationFrequency();
+  const T coarseTau = _coarse.getConverter().getLatticeRelaxationTime();
+  return (coarseTau - 0.25) / coarseTau;
 }
 
 template <typename T, template<typename> class DESCRIPTOR>
-- 
cgit v1.2.3