From 583fc4fb0dc9fbbdf0aa58c115b6ee05361ee413 Mon Sep 17 00:00:00 2001
From: Nima Samkhaniani <>
Date: Wed, 30 Mar 2022 11:37:49 +0100
Subject: [PATCH] ENH: atmPlantCanopyUSource: improve implicit behaviour (fixes
 #2343)

---
 .../atmPlantCanopyUSource/atmPlantCanopyUSource.C        | 9 +++++----
 1 file changed, 5 insertions(+), 4 deletions(-)

diff --git a/src/atmosphericModels/fvOptions/atmPlantCanopyUSource/atmPlantCanopyUSource.C b/src/atmosphericModels/fvOptions/atmPlantCanopyUSource/atmPlantCanopyUSource.C
index f1044030bd..15f448ff57 100644
--- a/src/atmosphericModels/fvOptions/atmPlantCanopyUSource/atmPlantCanopyUSource.C
+++ b/src/atmosphericModels/fvOptions/atmPlantCanopyUSource/atmPlantCanopyUSource.C
@@ -6,7 +6,7 @@
      \\/     M anipulation  |
 -------------------------------------------------------------------------------
     Copyright (C) 2020 ENERCON GmbH
-    Copyright (C) 2020-2021 OpenCFD Ltd.
+    Copyright (C) 2020-2022 OpenCFD Ltd.
 -------------------------------------------------------------------------------
 License
     This file is part of OpenFOAM.
@@ -28,6 +28,7 @@ License
 
 #include "atmPlantCanopyUSource.H"
 #include "addToRunTimeSelectionTable.H"
+#include "fvmSup.H"
 
 // * * * * * * * * * * * * * Static Member Functions * * * * * * * * * * * * //
 
@@ -99,7 +100,7 @@ void Foam::fv::atmPlantCanopyUSource::addSup
     if (V_ > VSMALL)
     {
         // (SP:Eq. 42)
-        eqn -= (plantCd_*leafAreaDensity_*mag(U))*U;
+        eqn -= fvm::Sp(plantCd_*leafAreaDensity_*mag(U), U);
     }
 }
 
@@ -115,7 +116,7 @@ void Foam::fv::atmPlantCanopyUSource::addSup
 
     if (V_ > VSMALL)
     {
-        eqn -= rho*(plantCd_*leafAreaDensity_*mag(U))*U;
+        eqn -= fvm::Sp(rho*plantCd_*leafAreaDensity_*mag(U), U);
     }
 }
 
@@ -132,7 +133,7 @@ void Foam::fv::atmPlantCanopyUSource::addSup
 
     if (V_ > VSMALL)
     {
-        eqn -= alpha*rho*(plantCd_*leafAreaDensity_*mag(U))*U;
+        eqn -= fvm::Sp(alpha*rho*plantCd_*leafAreaDensity_*mag(U), U);
     }
 }