From 67e2f83134b6d042cdf4058fbc19c078571c3f36 Mon Sep 17 00:00:00 2001
From: Hrvoje Jasak <h.jasak@wikki.co.uk>
Date: Tue, 29 Jan 2019 09:53:30 +0000
Subject: [PATCH] Bugfix: resolving read problem in initialisation of IB patch
 fields

---
 .../fixedValueIbFvPatchField.C                        |  2 +-
 .../basic/mixedIbFvPatchField/mixedIbFvPatchField.C   |  4 ++--
 ...movingImmersedBoundaryVelocityFvPatchVectorField.C |  2 +-
 .../immersedBoundaryFvsPatchField.C                   | 11 +++++++++--
 4 files changed, 13 insertions(+), 6 deletions(-)

diff --git a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/basic/fixedValueIbFvPatchField/fixedValueIbFvPatchField.C b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/basic/fixedValueIbFvPatchField/fixedValueIbFvPatchField.C
index 9cdc3e107..ba21b5e4c 100644
--- a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/basic/fixedValueIbFvPatchField/fixedValueIbFvPatchField.C
+++ b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/basic/fixedValueIbFvPatchField/fixedValueIbFvPatchField.C
@@ -66,7 +66,7 @@ fixedValueIbFvPatchField<Type>::fixedValueIbFvPatchField
     const dictionary& dict
 )
 :
-    fixedValueFvPatchField<Type>(p, iF, dict),
+    fixedValueFvPatchField<Type>(p, iF),   // Do not read mixed data
     immersedBoundaryFieldBase<Type>
     (
         p,
diff --git a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/basic/mixedIbFvPatchField/mixedIbFvPatchField.C b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/basic/mixedIbFvPatchField/mixedIbFvPatchField.C
index c8abaaea7..5f688461e 100644
--- a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/basic/mixedIbFvPatchField/mixedIbFvPatchField.C
+++ b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/basic/mixedIbFvPatchField/mixedIbFvPatchField.C
@@ -70,7 +70,7 @@ mixedIbFvPatchField<Type>::mixedIbFvPatchField
     const dictionary& dict
 )
 :
-    mixedFvPatchField<Type>(p, iF),
+    mixedFvPatchField<Type>(p, iF),   // Do not read mixed data
     immersedBoundaryFieldBase<Type>
     (
         p,
@@ -109,7 +109,7 @@ mixedIbFvPatchField<Type>::mixedIbFvPatchField
     // the patch is active
     // Initialise the value to avoid errors
     // HJ, 1/Dec/2017
-    if (this->ibPatch().ibPolyPatch().active())
+//    if (this->ibPatch().ibPolyPatch().active())
     {
         // Re-interpolate the data related to immersed boundary
         this->updateIbValues();
diff --git a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/derived/movingImmersedBoundaryVelocity/movingImmersedBoundaryVelocityFvPatchVectorField.C b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/derived/movingImmersedBoundaryVelocity/movingImmersedBoundaryVelocityFvPatchVectorField.C
index 1a61704a5..c79ac392d 100644
--- a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/derived/movingImmersedBoundaryVelocity/movingImmersedBoundaryVelocityFvPatchVectorField.C
+++ b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvPatchFields/derived/movingImmersedBoundaryVelocity/movingImmersedBoundaryVelocityFvPatchVectorField.C
@@ -65,7 +65,7 @@ movingImmersedBoundaryVelocityFvPatchVectorField
     const dictionary& dict
 )
 :
-    fixedValueFvPatchVectorField(p, iF, dict),
+    fixedValueFvPatchVectorField(p, iF),   // Do not read data
     immersedBoundaryFieldBase<vector>
     (
         p,
diff --git a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvsPatchFields/immersedBoundaryFvsPatchField/immersedBoundaryFvsPatchField.C b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvsPatchFields/immersedBoundaryFvsPatchField/immersedBoundaryFvsPatchField.C
index 255d29b3d..a13df8432 100644
--- a/src/immersedBoundary/immersedBoundary/immersedBoundaryFvsPatchFields/immersedBoundaryFvsPatchField/immersedBoundaryFvsPatchField.C
+++ b/src/immersedBoundary/immersedBoundary/immersedBoundaryFvsPatchFields/immersedBoundaryFvsPatchField/immersedBoundaryFvsPatchField.C
@@ -66,9 +66,11 @@ immersedBoundaryFvsPatchField<Type>::immersedBoundaryFvsPatchField
     const dictionary& dict
 )
 :
-    fvsPatchField<Type>(p, iF, dict),
+    fvsPatchField<Type>(p, iF),   // Do not read base data
     ibPatch_(refCast<const immersedBoundaryFvPatch>(p))
-{}
+{
+    operator=(pTraits<Type>::zero);
+}
 
 
 template<class Type>
@@ -149,6 +151,11 @@ void immersedBoundaryFvsPatchField<Type>::evaluate
 )
 {
     this->updateSize();
+
+    Field<Type>::operator=
+    (
+        Field<Type>(this->patch().size(), pTraits<Type>::zero)
+    );
 }