From f5d256cbe9c5ed71a08774da1dfa32ba1e11f674 Mon Sep 17 00:00:00 2001
From: Andreas Lauser <and@poware.org>
Date: Fri, 26 Feb 2016 15:03:41 +0100
Subject: [PATCH] ebos: minor changes to the EclWellManager

the only non-cosmetic change is the reduction of the magnitude of the
tubing head pressure from 10^100 to 10^30. this does not matter
normally, but if single precision floating point values are used
10^100 cannot be represented by such variables (and also 10^30 is
large enough by a fair margin: I doubt that one ever wants to simulate
the conditions in the core of a star using this code.)
---
 applications/ebos/eclwellmanager.hh | 10 +++++-----
 1 file changed, 5 insertions(+), 5 deletions(-)

diff --git a/applications/ebos/eclwellmanager.hh b/applications/ebos/eclwellmanager.hh
index 063cedcd7..b1f4bf061 100644
--- a/applications/ebos/eclwellmanager.hh
+++ b/applications/ebos/eclwellmanager.hh
@@ -237,7 +237,7 @@ public:
                 well->setTargetBottomHolePressure(injectProperties.BHPLimit);
 
                 // TODO
-                well->setTargetTubingHeadPressure(1e100);
+                well->setTargetTubingHeadPressure(1e30);
                 //well->setTargetTubingHeadPressure(injectProperties.THPLimit);
             }
 
@@ -304,7 +304,7 @@ public:
                 well->setTargetBottomHolePressure(producerProperties.BHPLimit);
 
                 // TODO
-                well->setTargetTubingHeadPressure(-1e100);
+                well->setTargetTubingHeadPressure(-1e30);
                 //well->setTargetTubingHeadPressure(producerProperties.THPLimit);
             }
         }
@@ -612,10 +612,11 @@ protected:
         // first, remove all wells from the reservoir
         model.clearAuxiliaryModules();
         auto wellIt = wells_.begin();
-        const auto wellEndIt = wells_.end();
+        const auto& wellEndIt = wells_.end();
         for (; wellIt != wellEndIt; ++wellIt)
             (*wellIt)->clear();
 
+        //////
         // tell the active wells which DOFs they contain
         const auto gridView = simulator_.gridManager().gridView();
 
@@ -648,15 +649,14 @@ protected:
 
                 wells.insert(eclWell);
             }
+            //////
         }
 
         // register all wells at the model as auxiliary equations
         auto wellIt2 = wells.begin();
         const auto& wellEndIt2 = wells.end();
         for (; wellIt2 != wellEndIt2; ++wellIt2)
-        {
             model.addAuxiliaryModule(*wellIt2);
-        }
     }
 
     void computeWellCompletionsMap_(unsigned reportStepIdx, WellCompletionsMap& cartesianIdxToCompletionMap)