Fix non-equilibrium scaling factor

1 files changed, 2 insertions, 1 deletions

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>