From 246bbfa7fb88da2b30bbbfb87cde90a55d7f122f Mon Sep 17 00:00:00 2001
From: Ray Speth <speth@mit.edu>
Date: Tue, 5 Mar 2013 17:02:50 +0000
Subject: [PATCH] [1D] Prevent removal of the grid point used for the fixed
 temperature

---
 src/oneD/refine.cpp | 5 +++++
 1 file changed, 5 insertions(+)

diff --git a/src/oneD/refine.cpp b/src/oneD/refine.cpp
index 042b2dac0..b0d6fcbd8 100644
--- a/src/oneD/refine.cpp
+++ b/src/oneD/refine.cpp
@@ -149,6 +149,11 @@ int Refiner::analyze(size_t n, const doublereal* z,
         if (j < n-2 && z[j+1]-z[j] > m_ratio * dz[j+1]) {
             m_keep[j] = 1;
         }
+
+        // Keep the point where the temperature is fixed
+        if (z[j] == m_domain->m_zfixed) {
+            m_keep[j] = 1;
+        }
     }
 
     // Don't allow pruning to remove multiple adjacent grid points