From 2b19bc57335b0d777691b0d423e63a0f688c5d1a Mon Sep 17 00:00:00 2001
From: wpbonelli <wbonelli@ucar.edu>
Date: Wed, 13 Dec 2023 11:56:31 -0500
Subject: [PATCH] switch to is_close for csub brg comparison with 0

---
 src/Model/GroundWaterFlow/gwf3csub8.f90 | 3 ++-
 1 file changed, 2 insertions(+), 1 deletion(-)

diff --git a/src/Model/GroundWaterFlow/gwf3csub8.f90 b/src/Model/GroundWaterFlow/gwf3csub8.f90
index 4f200c26d04..db87b8ff411 100644
--- a/src/Model/GroundWaterFlow/gwf3csub8.f90
+++ b/src/Model/GroundWaterFlow/gwf3csub8.f90
@@ -18,6 +18,7 @@ module GwfCsubModule
                              TABLEFT, TABCENTER, TABRIGHT, &
                              TABSTRING, TABUCSTRING, TABINTEGER, TABREAL
   use MemoryHelperModule, only: create_mem_path
+  use MathUtilModule, only: is_close
   use MessageModule, only: write_message
   use SmoothingModule, only: sQuadraticSaturation, &
                              sQuadraticSaturationDerivative, &
@@ -2874,7 +2875,7 @@ subroutine csub_fc(this, kiter, hold, hnew, matrix_sln, idxglo, rhs)
           rhs(node) = rhs(node) + rhsterm
           !
           ! -- calculate interbed water compressibility terms
-          if (this%brg /= DZERO .and. idelay == 0) then
+          if (.not. is_close(this%brg, DZERO) .and. idelay == 0) then
             call this%csub_nodelay_wcomp_fc(ib, node, tled, area, &
                                             hnew(node), hold(node), &
                                             hcof, rhsterm)