From 2bd05755232e9e7a0c2ece99cba5f5e07f6ed785 Mon Sep 17 00:00:00 2001
From: Tor Harald Sandve <tosa@norceresearch.no>
Date: Wed, 19 Oct 2022 09:16:35 +0200
Subject: [PATCH] implement dirichlet boundary conditions

---
 opm/models/discretization/common/tpfalinearizer.hh | 2 +-
 1 file changed, 1 insertion(+), 1 deletion(-)

diff --git a/opm/models/discretization/common/tpfalinearizer.hh b/opm/models/discretization/common/tpfalinearizer.hh
index ae13410c2..08c5e6db8 100644
--- a/opm/models/discretization/common/tpfalinearizer.hh
+++ b/opm/models/discretization/common/tpfalinearizer.hh
@@ -352,7 +352,7 @@ private:
                     }
                     const bool nonzero_massrate = massrate != VectorBlock(0.0);
                     if (free || nonzero_massrate) {
-                        const auto& exFluidState = problem_().initialFluidState(myIdx);
+                        const auto& exFluidState = problem_().boundaryFluidState(myIdx, dir_id);
                         BoundaryConditionData bcdata{free ? BCType::FREE : BCType::RATE,
                                                      massrate,
                                                      exFluidState.pvtRegionIndex(),