From 82de75b26dc6537dc084b0a14fd96d6576c1e6e3 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Wed, 29 May 2013 10:39:14 +0100
Subject: [PATCH 1/2] BUG: Updated particle SRF force - mantis #716

---
 .../submodels/Kinematic/ParticleForces/SRF/SRFForce.C       | 6 ++++--
 1 file changed, 4 insertions(+), 2 deletions(-)

diff --git a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C
index 758e58743b..eed0c800d0 100644
--- a/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C
+++ b/src/lagrangian/intermediate/submodels/Kinematic/ParticleForces/SRF/SRFForce.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2011 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2011-2013 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -95,7 +95,9 @@ Foam::forceSuSp Foam::SRFForce<CloudType>::calcNonCoupled
     const vector& r = p.position();
 
     // Coriolis and centrifugal acceleration terms
-    value.Su() = mass*(2.0*(p.U() ^ omega) + (omega ^ (r ^ omega)));
+    value.Su() =
+        mass*(1.0 - p.rhoc()/p.rho())
+       *(2.0*(p.U() ^ omega) + (omega ^ (r ^ omega)));
 
     return value;
 }

From 5ed77670542df5c22fa89c845a909fba094719f3 Mon Sep 17 00:00:00 2001
From: andy <andy>
Date: Wed, 29 May 2013 11:00:40 +0100
Subject: [PATCH 2/2] BUG: Updated fieldAverage function object for restart -
 mantis 854

---
 .../fieldAverage/fieldAverage/fieldAverage.C  | 30 +++++++------------
 1 file changed, 10 insertions(+), 20 deletions(-)

diff --git a/src/postProcessing/functionObjects/field/fieldAverage/fieldAverage/fieldAverage.C b/src/postProcessing/functionObjects/field/fieldAverage/fieldAverage/fieldAverage.C
index d7ba9df79f..5a73765f93 100644
--- a/src/postProcessing/functionObjects/field/fieldAverage/fieldAverage/fieldAverage.C
+++ b/src/postProcessing/functionObjects/field/fieldAverage/fieldAverage/fieldAverage.C
@@ -144,11 +144,21 @@ void Foam::fieldAverage::initialize()
             }
         }
     }
+
+    // ensure first averaging works unconditionally
+    prevTimeIndex_ = -1;
+
+    initialised_ = true;
 }
 
 
 void Foam::fieldAverage::calcAverages()
 {
+    if (!initialised_)
+    {
+        initialize();
+    }
+
     const label currentTimeIndex =
         static_cast<const fvMesh&>(obr_).time().timeIndex();
 
@@ -374,16 +384,6 @@ void Foam::fieldAverage::execute()
 {
     if (active_)
     {
-        if (!initialised_)
-        {
-            initialize();
-
-            // ensure first averaging works unconditionally
-            prevTimeIndex_ = -1;
-
-            initialised_ = true;
-        }
-
         calcAverages();
     }
 }
@@ -397,16 +397,6 @@ void Foam::fieldAverage::write()
 {
     if (active_)
     {
-        if (!initialised_)
-        {
-            initialize();
-
-            // ensure first averaging works unconditionally
-            prevTimeIndex_ = -1;
-
-            initialised_ = true;
-        }
-
         calcAverages();
         writeAverages();
         writeAveragingProperties();