From d4f6e184beaec65aa4d16ab8ca9cdc40ddb9f0a3 Mon Sep 17 00:00:00 2001
From: Henry <Henry>
Date: Sat, 18 Apr 2015 19:31:31 +0100
Subject: [PATCH] prghPressureFvPatchScalarField: Add support for hRef

---
 .../prghPressure/prghPressureFvPatchScalarField.C  | 14 ++++++++++++--
 .../prghPressure/prghPressureFvPatchScalarField.H  |  7 ++++---
 2 files changed, 16 insertions(+), 5 deletions(-)

diff --git a/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.C b/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.C
index 2338247e7..c86986f71 100644
--- a/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.C
+++ b/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.C
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2013-2015 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -152,7 +152,17 @@ void Foam::prghPressureFvPatchScalarField::updateCoeffs()
     const uniformDimensionedVectorField& g =
         db().lookupObject<uniformDimensionedVectorField>("g");
 
-    operator==(p_ - rhop*((g.value() & patch().Cf())));
+    const uniformDimensionedScalarField& hRef =
+        db().lookupObject<uniformDimensionedScalarField>("hRef");
+
+    dimensionedScalar ghRef
+    (
+        mag(g.value()) > SMALL
+      ? g & (cmptMag(g.value())/mag(g.value()))*hRef
+      : dimensionedScalar("ghRef", g.dimensions()*dimLength, 0)
+    );
+
+    operator==(p_ - rhop*((g.value() & patch().Cf()) - ghRef.value()));
 
     fixedValueFvPatchScalarField::updateCoeffs();
 }
diff --git a/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.H b/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.H
index 9431ae325..bad53659f 100644
--- a/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.H
+++ b/src/finiteVolume/fields/fvPatchFields/derived/prghPressure/prghPressureFvPatchScalarField.H
@@ -2,7 +2,7 @@
   =========                 |
   \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
    \\    /   O peration     |
-    \\  /    A nd           | Copyright (C) 2013 OpenFOAM Foundation
+    \\  /    A nd           | Copyright (C) 2013-2015 OpenFOAM Foundation
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
 License
@@ -32,7 +32,7 @@ Description
     calculated as:
 
         \f[
-            p_rgh = p - \rho g h
+            p_rgh = p - \rho g (h - hRef)
         \f]
 
     where
@@ -40,8 +40,9 @@ Description
         p_rgh   | Pseudo hydrostatic pressure [Pa]
         p       | Static pressure [Pa]
         h       | Height in the opposite direction to gravity
+        hRef    | Reference height in the opposite direction to gravity
         \rho    | density
-        g       | acceleration due to gravity [m/s2]
+        g       | acceleration due to gravity [m/s^2]
     \endtable
 
     \heading Patch usage