diff --git a/.integrated_tests.yaml b/.integrated_tests.yaml
index 3de3368090d..3a15c942b22 100644
--- a/.integrated_tests.yaml
+++ b/.integrated_tests.yaml
@@ -1,6 +1,6 @@
 baselines:
   bucket: geosx
-  baseline: integratedTests/baseline_integratedTests-pr3836-17046-2e89f64
+  baseline: integratedTests/baseline_integratedTests-pr4100-17125-aafd403
 
 allow_fail:
   all: ''
diff --git a/BASELINE_NOTES.md b/BASELINE_NOTES.md
index 6f77c74b57b..274d149e606 100644
--- a/BASELINE_NOTES.md
+++ b/BASELINE_NOTES.md
@@ -5,8 +5,10 @@ This file is designed to track changes to the integrated test baselines.
 Any developer who updates the baseline ID in the .integrated_tests.yaml file is expected to create an entry in this file with the pull request number, date, and their justification for rebaselining.
 These notes should be in reverse-chronological order, and use the following time format: (YYYY-MM-DD).
 
+PR #4100 (2026-07-22) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr4100-17125-aafd403.tar.gz>
+Fix to include all well residuals in convergence check.
+
 PR #3836 (2026-05-20) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr3836-17046-2e89f64.tar.gz>
-=====================
 Added statistics `Group` objects for each statistics `Task` instance
 
 PR #4040 (2026-06-16) <https://storage.googleapis.com/geosx/integratedTests/baseline_integratedTests-pr4040-16993-1393f80.tar.gz>
diff --git a/src/coreComponents/physicsSolvers/fluidFlow/wells/CompositionalMultiphaseWell.cpp b/src/coreComponents/physicsSolvers/fluidFlow/wells/CompositionalMultiphaseWell.cpp
index 6a72af7a90c..8870bddcd7d 100644
--- a/src/coreComponents/physicsSolvers/fluidFlow/wells/CompositionalMultiphaseWell.cpp
+++ b/src/coreComponents/physicsSolvers/fluidFlow/wells/CompositionalMultiphaseWell.cpp
@@ -1443,7 +1443,6 @@ CompositionalMultiphaseWell::calculateResidualNorm( real64 const & time_n,
       {
         for( integer i = 0; i < numNorm; ++i )
         {
-          localResidualNorm[i] = 0.0;
           localResidualNormalizer[i] = m_nonlinearSolverParameters.m_minNormalizer;
         }
       }