From e8766829cdc39f0903f6b4bf7f4b2228d16991c5 Mon Sep 17 00:00:00 2001
From: Daniel <2929574+danielque@users.noreply.github.com>
Date: Wed, 29 Jul 2020 17:26:57 +0200
Subject: [PATCH] fix incompatible dimensions in volWeightedAverage

ref() gives access to dimensioned field reference
primitiveFieldRef() gives access to primitive type field reference
(both actually return a 'this'-pointer reference)
---
 .../forceModel/volWeightedAverage/volWeightedAverage.C          | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/src/lagrangian/cfdemParticle/subModels/forceModel/volWeightedAverage/volWeightedAverage.C b/src/lagrangian/cfdemParticle/subModels/forceModel/volWeightedAverage/volWeightedAverage.C
index c5cffe6f..9f76f1a7 100644
--- a/src/lagrangian/cfdemParticle/subModels/forceModel/volWeightedAverage/volWeightedAverage.C
+++ b/src/lagrangian/cfdemParticle/subModels/forceModel/volWeightedAverage/volWeightedAverage.C
@@ -171,7 +171,7 @@ void volWeightedAverage::setForce() const
             MPI_Allreduce(&totVol, &totVol_all, 3, MPI_DOUBLE, MPI_SUM, MPI_COMM_WORLD);
             integralValue = gSum(scalarFields_[i]);
             volWeightedAverage = integralValue / (totVol_all+SMALL);
-            scalarFields_[i].ref() = volWeightedAverage;
+            scalarFields_[i].primitiveFieldRef() = volWeightedAverage;
 
             if(verbose_)
             {