From 769a21d177eea24197a1e7bad9d0cd6625044740 Mon Sep 17 00:00:00 2001
From: Malcolm  Ross <malcolmgti@gmail.com>
Date: Mon, 2 Sep 2024 14:41:00 -0500
Subject: [PATCH] SBT for GEOPHIRES Complete. Added tests from paper
 Techno-Economic Performance of Eavor-Loop 2.0 Koenraad F. Beckers, Henry E.
 Johnston (2022)

---
 src/geophires_x/SBTReservoir.py     |  13 +-
 tests/examples/example_SBT_Hi_T.out | 224 ++++++++++++++++++++++++++++
 tests/examples/example_SBT_Hi_T.txt |  48 ++++++
 tests/examples/example_SBT_Lo_T.out | 224 ++++++++++++++++++++++++++++
 tests/examples/example_SBT_Lo_T.txt |  50 +++++++
 5 files changed, 556 insertions(+), 3 deletions(-)
 create mode 100644 tests/examples/example_SBT_Hi_T.out
 create mode 100644 tests/examples/example_SBT_Hi_T.txt
 create mode 100644 tests/examples/example_SBT_Lo_T.out
 create mode 100644 tests/examples/example_SBT_Lo_T.txt

diff --git a/src/geophires_x/SBTReservoir.py b/src/geophires_x/SBTReservoir.py
index 6821aed4..41856951 100644
--- a/src/geophires_x/SBTReservoir.py
+++ b/src/geophires_x/SBTReservoir.py
@@ -10,7 +10,7 @@
 import geophires_x.Model as Model
 from .CylindricalReservoir import CylindricalReservoir
 from .OptionList import FlowrateModel, InjectionTemperatureModel, Configuration
-from .Parameter import intParameter, floatParameter, OutputParameter, ReadParameter, strParameter
+from .Parameter import intParameter, floatParameter, OutputParameter, ReadParameter, strParameter, boolParameter
 from .Reservoir import Reservoir
 from .Units import *
 import sys
@@ -267,7 +267,14 @@ def __init__(self, model: Model):
             PreferredUnits=TimeUnit.SECOND,
             CurrentUnits=TimeUnit.SECOND,
             ErrMessage='assume default for Initial to Final Timestep Transition (9900 seconds)',
-            ToolTipText='The time in secs at which the time arrays switches from closely spaced linear to logrithmatic'
+            ToolTipText='The time in secs at which the time arrays switches from closely spaced linear to logarithmic'
+        )
+        self.generate_wireframe_graphics = self.ParameterDict[self.generate_wireframe_graphics.Name] = boolParameter(
+            'SBT Generate Wireframe Graphics',
+            DefaultValue=False,
+            UnitType=Units.NONE,
+            ErrMessage='assume default for SBT Generate Wireframe Graphics (False)',
+            ToolTipText='Switch to control the generation of a wireframe drawing of a SBT wells configuration'
         )
 
         sclass = str(__class__).replace("<class \'", "")
@@ -1130,7 +1137,7 @@ def Calculate_Uloop(self, model):
                                                                                            model.wellbores.vertical_section_length.value,
                                                                                            model.wellbores.lateral_inclination_angle.value * np.pi / 180.0,
                                                                                            model.wellbores.vertical_wellbore_spacing.value,
-                                                                                           False)
+                                                                                           self.generate_wireframe_graphics.value)
 
         # Merge x-, y-, and z-coordinates
         x = np.concatenate((xinj, xprod))
diff --git a/tests/examples/example_SBT_Hi_T.out b/tests/examples/example_SBT_Hi_T.out
new file mode 100644
index 00000000..30fc221c
--- /dev/null
+++ b/tests/examples/example_SBT_Hi_T.out
@@ -0,0 +1,224 @@
+                               *****************
+                               ***CASE REPORT***
+                               *****************
+
+Simulation Metadata
+----------------------
+ GEOPHIRES Version: 3.5.2
+ Simulation Date: 2024-09-02
+ Simulation Time:  13:59
+ Calculation Time:     14.837 sec
+
+                           ***SUMMARY OF RESULTS***
+
+      End-Use Option: Electricity
+      Average Net Electricity Production:                     9.62 MW
+      Electricity breakeven price:                           18.23 cents/kWh
+      Number of production wells:                             1
+      Number of injection wells:                              1
+      Flowrate per production well:                          80.0 kg/sec
+      Well depth (or total length, if not vertical):          5.8 kilometer
+      Geothermal gradient:                                     60 degC/km
+
+
+                           ***ECONOMIC PARAMETERS***
+
+      Economic Model  = BICYCLE
+      Accrued financing during construction:                  0.00
+      Project lifetime:                                      30 yr
+      Capacity factor:                                       90.0 %
+      Project NPV:                                           0.96 MUSD
+      Project IRR:                                           6.32 %
+      Project VIR=PI=PIR:                                    1.01
+      Project MOIC:                                          0.65
+      Project Payback Period:                               13.59 yr
+      Estimated Jobs Created:                                 0
+
+                          ***ENGINEERING PARAMETERS***
+
+      Number of Production Wells:                             1
+      Number of Injection Wells:                              1
+      Well depth (or total length, if not vertical):          5.8 kilometer
+      Water loss rate:                                        0.0
+      Pump efficiency:                                       75.0
+      Injection temperature:                                 56.9 degC
+      Production Wellbore heat transmission calculated with Ramey's model
+      Average production well temperature drop:               0.0 degC
+      Flowrate per production well:                          80.0 kg/sec
+      Injection well casing ID:                               8.500 in
+      Production well casing ID:                              8.500 in
+      Number of times redrilling:                             0
+      Power plant type:                                       Supercritical ORC
+
+
+                         ***RESOURCE CHARACTERISTICS***
+
+      Maximum reservoir temperature:                        400.0 degC
+      Number of segments:                                     1
+      Geothermal gradient:                                        60 degC/km
+
+
+                           ***RESERVOIR PARAMETERS***
+
+The AGS models contain an intrinsic reservoir model that doesn't expose values that can be used in extensive reporting.
+
+
+                           ***RESERVOIR SIMULATION RESULTS***
+
+      Maximum Production Temperature:                       313.6 degC
+      Average Production Temperature:                       218.7 degC
+      Minimum Production Temperature:                       205.9 degC
+      Initial Production Temperature:                       313.6 degC
+The AGS models contain an intrinsic reservoir model that doesn't expose values that can be used in extensive reporting.
+
+
+                          ***CAPITAL COSTS (M$)***
+
+         Drilling and completion costs:                     104.87 MUSD
+             Drilling and completion costs per vertical production well:         6.19 MUSD
+             Drilling and completion costs per vertical injection well:          6.19 MUSD
+             Drilling and completion costs per non-vertical sections:           92.48 MUSD
+         Stimulation costs:                                   0.00 MUSD
+         Surface power plant costs:                          52.76 MUSD
+         Field gathering system costs:                        0.97 MUSD
+         Total surface equipment costs:                      53.72 MUSD
+         Exploration costs:                                   0.00 MUSD
+      Total capital costs:                                  158.59 MUSD
+
+
+                ***OPERATING AND MAINTENANCE COSTS (M$/yr)***
+
+         Wellfield maintenance costs:                         1.44 MUSD/yr
+         Power plant maintenance costs:                       1.93 MUSD/yr
+         Water costs:                                         0.00 MUSD/yr
+      Total operating and maintenance costs:                  3.37 MUSD/yr
+
+
+                           ***SURFACE EQUIPMENT SIMULATION RESULTS***
+
+      Initial geofluid availability:                          0.46 MW/(kg/s)
+      Maximum Total Electricity Generation:                  17.60 MW
+      Average Total Electricity Generation:                   9.62 MW
+      Minimum Total Electricity Generation:                   8.40 MW
+      Initial Total Electricity Generation:                  17.60 MW
+      Maximum Net Electricity Generation:                    17.60 MW
+      Average Net Electricity Generation:                     9.62 MW
+      Minimum Net Electricity Generation:                     8.40 MW
+      Initial Net Electricity Generation:                    17.60 MW
+      Average Annual Total Electricity Generation:           75.03 GWh
+      Average Annual Net Electricity Generation:             75.03 GWh
+      Average Pumping Power:                                  0.00 MW
+
+                            ************************************************************
+                            *  HEATING, COOLING AND/OR ELECTRICITY PRODUCTION PROFILE  *
+                            ************************************************************
+  YEAR       THERMAL               GEOFLUID               PUMP               NET               FIRST LAW
+             DRAWDOWN             TEMPERATURE             POWER             POWER              EFFICIENCY
+                                     (degC)               (MW)              (MW)                  (%)
+   1           1.0000                313.57               0.0000           17.6004               19.8758
+   2           0.7703                241.55               0.0000           11.9582               18.7707
+   3           0.7466                234.12               0.0000           11.2166               18.3442
+   4           0.7332                229.89               0.0000           10.7920               18.0811
+   5           0.7237                226.92               0.0000           10.4929               17.8876
+   6           0.7163                224.62               0.0000           10.2615               17.7333
+   7           0.7104                222.76               0.0000           10.0753               17.6061
+   8           0.7054                221.20               0.0000            9.9182               17.4968
+   9           0.7012                219.86               0.0000            9.7845               17.4022
+  10           0.6975                218.71               0.0000            9.6694               17.3197
+  11           0.6942                217.67               0.0000            9.5653               17.2443
+  12           0.6911                216.72               0.0000            9.4701               17.1745
+  13           0.6883                215.83               0.0000            9.3813               17.1088
+  14           0.6856                214.98               0.0000            9.2967               17.0456
+  15           0.6832                214.23               0.0000            9.2222               16.9895
+  16           0.6809                213.52               0.0000            9.1514               16.9358
+  17           0.6788                212.85               0.0000            9.0857               16.8855
+  18           0.6768                212.22               0.0000            9.0230               16.8373
+  19           0.6748                211.60               0.0000            8.9614               16.7895
+  20           0.6729                211.01               0.0000            8.9031               16.7442
+  21           0.6712                210.46               0.0000            8.8491               16.7018
+  22           0.6695                209.94               0.0000            8.7974               16.6611
+  23           0.6679                209.43               0.0000            8.7474               16.6214
+  24           0.6663                208.92               0.0000            8.6975               16.5817
+  25           0.6647                208.42               0.0000            8.6478               16.5420
+  26           0.6631                207.94               0.0000            8.6005               16.5039
+  27           0.6617                207.50               0.0000            8.5570               16.4688
+  28           0.6603                207.05               0.0000            8.5136               16.4335
+  29           0.6590                206.64               0.0000            8.4736               16.4009
+  30           0.6577                206.24               0.0000            8.4338               16.3684
+
+
+                              *******************************************************************
+                              *  ANNUAL HEATING, COOLING AND/OR ELECTRICITY PRODUCTION PROFILE  *
+                              *******************************************************************
+  YEAR             ELECTRICITY                   HEAT                RESERVOIR            PERCENTAGE OF
+                    PROVIDED                   EXTRACTED            HEAT CONTENT        TOTAL HEAT MINED
+                   (GWh/year)                  (GWh/year)            (10^15 J)                 (%)
+   1                 133.2                       673.7                 82.77                   2.85
+   2                  91.1                       491.1                 81.01                   4.92
+   3                  86.7                       476.0                 79.29                   6.93
+   4                  83.8                       466.3                 77.61                   8.90
+   5                  81.8                       459.2                 75.96                  10.84
+   6                  80.2                       453.6                 74.33                  12.76
+   7                  78.8                       449.0                 72.71                  14.66
+   8                  77.7                       445.1                 71.11                  16.54
+   9                  76.7                       441.7                 69.52                  18.40
+  10                  75.8                       438.7                 67.94                  20.26
+  11                  75.0                       436.0                 66.37                  22.10
+  12                  74.3                       433.5                 64.81                  23.93
+  13                  73.6                       431.1                 63.26                  25.75
+  14                  73.0                       429.0                 61.71                  27.57
+  15                  72.4                       427.0                 60.18                  29.37
+  16                  71.9                       425.1                 58.65                  31.17
+  17                  71.4                       423.4                 57.12                  32.96
+  18                  70.9                       421.7                 55.60                  34.74
+  19                  70.4                       420.0                 54.09                  36.51
+  20                  70.0                       418.5                 52.58                  38.28
+  21                  69.6                       417.0                 51.08                  40.04
+  22                  69.2                       415.6                 49.59                  41.80
+  23                  68.8                       414.2                 48.10                  43.55
+  24                  68.4                       412.8                 46.61                  45.29
+  25                  68.0                       411.5                 45.13                  47.03
+  26                  67.6                       410.2                 43.65                  48.77
+  27                  67.3                       409.0                 42.18                  50.49
+  28                  67.0                       407.9                 40.71                  52.22
+  29                  66.6                       406.8                 39.25                  53.94
+  30                  49.8                       304.4                 38.15                  55.22
+
+
+                             ********************************
+                             *  REVENUE & CASHFLOW PROFILE  *
+                             ********************************
+Year            Electricity             |            Heat                  |           Cooling                 |         Carbon                    |          Project
+Since     Price   Ann. Rev.  Cumm. Rev. |   Price   Ann. Rev.   Cumm. Rev. |  Price   Ann. Rev.   Cumm. Rev.   |   Price   Ann. Rev.   Cumm. Rev.  | OPEX    Net Rev.      Net Cashflow
+Start    (cents/kWh)(MUSD/yr) (MUSD)    |(cents/kWh) (MUSD/yr)    (MUSD)   |(cents/kWh) (MUSD/yr)    (MUSD)    |(USD/tonne) (MUSD/yr)    (MUSD)    |(MUSD/yr) (MUSD/yr)    (MUSD)
+________________________________________________________________________________________________________________________________________________________________________________________
+  1      0.00          -158.59   0.00     |    0.00     0.00         0.00    |    0.00     0.00         0.00     |    0.00     0.00         0.00     |  0.00     -158.59     -158.59
+  2     19.00          21.93  25.31     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     21.93     -136.66
+  3     19.00          13.93  42.61     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     13.93     -122.73
+  4     19.00          13.09  59.07     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     13.09     -109.64
+  5     19.00          12.56  75.00     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     12.56     -97.08
+  6     19.00          12.16  90.54     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     12.16     -84.92
+  7     19.00          11.86  105.77     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     11.86     -73.06
+  8     19.00          11.60  120.74     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     11.60     -61.46
+  9     19.00          11.38  135.50     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     11.38     -50.08
+ 10     19.00          11.20  150.07     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     11.20     -38.89
+ 11     19.00          11.03  164.47     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     11.03     -27.85
+ 12     19.00          10.88  178.73     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.88     -16.97
+ 13     19.00          10.75  192.85     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.75     -6.23
+ 14     19.00          10.61  206.83     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.61      4.39
+ 15     19.00          10.50  220.70     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.50     14.88
+ 16     19.00          10.39  234.46     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.39     25.27
+ 17     19.00          10.28  248.12     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.28     35.56
+ 18     19.00          10.19  261.69     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.19     45.75
+ 19     19.00          10.10  275.16     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.10     55.84
+ 20     19.00          10.01  288.53     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37     10.01     65.85
+ 21     19.00           9.92  301.83     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.92     75.77
+ 22     19.00           9.84  315.05     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.84     85.61
+ 23     19.00           9.77  328.19     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.77     95.38
+ 24     19.00           9.69  341.25     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.69     105.07
+ 25     19.00           9.62  354.24     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.62     114.69
+ 26     19.00           9.54  367.16     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.54     124.23
+ 27     19.00           9.48  380.01     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.48     133.71
+ 28     19.00           9.41  392.80     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.41     143.12
+ 29     19.00           9.35  405.52     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.35     152.47
+ 30     19.00           9.29  418.18     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  3.37      9.29     161.76
diff --git a/tests/examples/example_SBT_Hi_T.txt b/tests/examples/example_SBT_Hi_T.txt
new file mode 100644
index 00000000..0aaf7e92
--- /dev/null
+++ b/tests/examples/example_SBT_Hi_T.txt
@@ -0,0 +1,48 @@
+************reservoir*************
+Reservoir Model, 8
+Reservoir Depth, 7.5 kilometer
+Gradient 1, 60
+Reservoir Volume Option, 4
+Reservoir Volume, 8136407202.64, pi * 1 mi * 1 mi * 1000
+Reservoir Heat Capacity, 1112
+Reservoir Density, 2663
+Reservoir Thermal Conductivity, 2.5
+
+Lateral Endpoint Depth, 7.5 kilometer
+Lateral Inclination Angle, 30.5
+Junction Depth, 4.2 kilometer
+Vertical Section Length, 4.2 kilometer
+Number of Multilateral Sections, 12
+SBT Accuracy Desired, 5
+Lateral Spacing, 75
+Discretization Length, 250
+SBT Initial Timestep Count, 5
+SBT Initial to Final Timestep Transition, 10000, seconds = one year
+SBT Final Timestep Count, 120, ~4 per year
+
+****************Wellbore***********
+Is AGS, True
+Well Geometry Configuration, 5
+Number of Production Wells, 1
+Number of Injection Wells, 1
+Production Well Diameter, 8.5
+Injection Well Diameter, 8.5
+Nonvertical Wellbore Diameter, 0.216
+Production Flow Rate per Well, 80
+Reservoir Impedance, 1E-4, there is very little impedance in an AGS system
+Multilaterals Cased, False
+
+****************Surface plant************
+End-Use Option, 1
+Power Plant Type, 2
+Plant Lifetime, 30
+Ambient Temperature, 10
+Surface Temperature, 10
+Injection Temperature, 60
+
+***************Economics****************
+Economic Model, 3
+Reservoir Stimulation Capital Cost, 0
+Exploration Capital Cost, 0
+Starting Electricity Sale Price, 0.19
+Ending Electricity Sale Price, 0.19
diff --git a/tests/examples/example_SBT_Lo_T.out b/tests/examples/example_SBT_Lo_T.out
new file mode 100644
index 00000000..1390f1bd
--- /dev/null
+++ b/tests/examples/example_SBT_Lo_T.out
@@ -0,0 +1,224 @@
+                               *****************
+                               ***CASE REPORT***
+                               *****************
+
+Simulation Metadata
+----------------------
+ GEOPHIRES Version: 3.5.2
+ Simulation Date: 2024-09-02
+ Simulation Time:  14:35
+ Calculation Time:     45.642 sec
+
+                           ***SUMMARY OF RESULTS***
+
+      End-Use Option: Electricity
+      Average Net Electricity Production:                     0.00 MW
+      Electricity breakeven price:                        18339.44 cents/kWh
+      Number of production wells:                             1
+      Number of injection wells:                              1
+      Flowrate per production well:                           6.0 kg/sec
+      Well depth (or total length, if not vertical):          2.5 kilometer
+      Geothermal gradient:                                  31.25 degC/km
+
+
+                           ***ECONOMIC PARAMETERS***
+
+      Economic Model  = BICYCLE
+      Accrued financing during construction:                  0.00
+      Project lifetime:                                      30 yr
+      Capacity factor:                                       90.0 %
+      Project NPV:                                         -47.18 MUSD
+      Project IRR:                                           0.00 %
+      Project VIR=PI=PIR:                                   -0.22
+      Project MOIC:                                         -1.00
+      Project Payback Period:                          N/A
+      Estimated Jobs Created:                                 0
+
+                          ***ENGINEERING PARAMETERS***
+
+      Number of Production Wells:                             1
+      Number of Injection Wells:                              1
+      Well depth (or total length, if not vertical):          2.5 kilometer
+      Water loss rate:                                        0.0
+      Pump efficiency:                                       75.0
+      Injection temperature:                                 24.0 degC
+      Production Wellbore heat transmission calculated with Ramey's model
+      Average production well temperature drop:               0.0 degC
+      Flowrate per production well:                           6.0 kg/sec
+      Injection well casing ID:                               8.500 in
+      Production well casing ID:                              8.500 in
+      Number of times redrilling:                             0
+      Power plant type:                                       Supercritical ORC
+
+
+                         ***RESOURCE CHARACTERISTICS***
+
+      Maximum reservoir temperature:                        400.0 degC
+      Number of segments:                                     1
+      Geothermal gradient:                                     31.25 degC/km
+
+
+                           ***RESERVOIR PARAMETERS***
+
+The AGS models contain an intrinsic reservoir model that doesn't expose values that can be used in extensive reporting.
+
+
+                           ***RESERVOIR SIMULATION RESULTS***
+
+      Maximum Production Temperature:                        59.0 degC
+      Average Production Temperature:                        58.9 degC
+      Minimum Production Temperature:                        58.4 degC
+      Initial Production Temperature:                        59.0 degC
+The AGS models contain an intrinsic reservoir model that doesn't expose values that can be used in extensive reporting.
+
+
+                          ***CAPITAL COSTS (M$)***
+
+         Drilling and completion costs:                      37.58 MUSD
+             Drilling and completion costs per vertical production well:         3.28 MUSD
+             Drilling and completion costs per vertical injection well:          3.28 MUSD
+             Drilling and completion costs per non-vertical sections:           31.02 MUSD
+         Stimulation costs:                                   0.00 MUSD
+         Surface power plant costs:                           0.03 MUSD
+         Field gathering system costs:                        0.97 MUSD
+         Total surface equipment costs:                       0.99 MUSD
+         Exploration costs:                                   0.00 MUSD
+      Total capital costs:                                   38.57 MUSD
+
+
+                ***OPERATING AND MAINTENANCE COSTS (M$/yr)***
+
+         Wellfield maintenance costs:                         0.45 MUSD/yr
+         Power plant maintenance costs:                       0.20 MUSD/yr
+         Water costs:                                         0.00 MUSD/yr
+      Total operating and maintenance costs:                  0.65 MUSD/yr
+
+
+                           ***SURFACE EQUIPMENT SIMULATION RESULTS***
+
+      Initial geofluid availability:                          0.02 MW/(kg/s)
+      Maximum Total Electricity Generation:                   0.00 MW
+      Average Total Electricity Generation:                   0.00 MW
+      Minimum Total Electricity Generation:                   0.00 MW
+      Initial Total Electricity Generation:                   0.00 MW
+      Maximum Net Electricity Generation:                     0.00 MW
+      Average Net Electricity Generation:                     0.00 MW
+      Minimum Net Electricity Generation:                     0.00 MW
+      Initial Net Electricity Generation:                     0.00 MW
+      Average Annual Total Electricity Generation:            0.02 GWh
+      Average Annual Net Electricity Generation:              0.02 GWh
+      Average Pumping Power:                                  0.00 MW
+
+                            ************************************************************
+                            *  HEATING, COOLING AND/OR ELECTRICITY PRODUCTION PROFILE  *
+                            ************************************************************
+  YEAR       THERMAL               GEOFLUID               PUMP               NET               FIRST LAW
+             DRAWDOWN             TEMPERATURE             POWER             POWER              EFFICIENCY
+                                     (degC)               (MW)              (MW)                  (%)
+   1           1.0000                 58.97               0.0000            0.0024                0.2754
+   2           0.9908                 58.44               0.0000            0.0019                0.2282
+   3           0.9960                 58.74               0.0000            0.0022                0.2546
+   4           0.9979                 58.85               0.0000            0.0023                0.2647
+   5           0.9989                 58.91               0.0000            0.0023                0.2698
+   6           0.9994                 58.94               0.0000            0.0024                0.2725
+   7           0.9998                 58.96               0.0000            0.0024                0.2742
+   8           0.9999                 58.97               0.0000            0.0024                0.2749
+   9           1.0000                 58.97               0.0000            0.0024                0.2752
+  10           1.0000                 58.97               0.0000            0.0024                0.2754
+  11           1.0000                 58.97               0.0000            0.0024                0.2754
+  12           1.0000                 58.97               0.0000            0.0024                0.2752
+  13           0.9999                 58.97               0.0000            0.0024                0.2749
+  14           0.9998                 58.96               0.0000            0.0024                0.2743
+  15           0.9996                 58.95               0.0000            0.0024                0.2732
+  16           0.9995                 58.94               0.0000            0.0024                0.2726
+  17           0.9993                 58.94               0.0000            0.0024                0.2720
+  18           0.9992                 58.93               0.0000            0.0024                0.2713
+  19           0.9991                 58.92               0.0000            0.0023                0.2706
+  20           0.9989                 58.91               0.0000            0.0023                0.2698
+  21           0.9988                 58.90               0.0000            0.0023                0.2690
+  22           0.9986                 58.89               0.0000            0.0023                0.2681
+  23           0.9984                 58.88               0.0000            0.0023                0.2673
+  24           0.9982                 58.87               0.0000            0.0023                0.2664
+  25           0.9981                 58.86               0.0000            0.0023                0.2654
+  26           0.9979                 58.85               0.0000            0.0023                0.2645
+  27           0.9977                 58.84               0.0000            0.0023                0.2635
+  28           0.9975                 58.83               0.0000            0.0023                0.2625
+  29           0.9973                 58.82               0.0000            0.0023                0.2615
+  30           0.9971                 58.81               0.0000            0.0022                0.2606
+
+
+                              *******************************************************************
+                              *  ANNUAL HEATING, COOLING AND/OR ELECTRICITY PRODUCTION PROFILE  *
+                              *******************************************************************
+  YEAR             ELECTRICITY                   HEAT                RESERVOIR            PERCENTAGE OF
+                    PROVIDED                   EXTRACTED            HEAT CONTENT        TOTAL HEAT MINED
+                   (GWh/year)                  (GWh/year)            (10^15 J)                 (%)
+   1                   0.0                         6.8                  0.86                   2.79
+   2                   0.0                         6.8                  0.83                   5.55
+   3                   0.0                         6.8                  0.81                   8.33
+   4                   0.0                         6.8                  0.78                  11.12
+   5                   0.0                         6.8                  0.76                  13.91
+   6                   0.0                         6.8                  0.73                  16.70
+   7                   0.0                         6.8                  0.71                  19.50
+   8                   0.0                         6.8                  0.68                  22.29
+   9                   0.0                         6.8                  0.66                  25.08
+  10                   0.0                         6.8                  0.64                  27.88
+  11                   0.0                         6.8                  0.61                  30.67
+  12                   0.0                         6.8                  0.59                  33.47
+  13                   0.0                         6.8                  0.56                  36.26
+  14                   0.0                         6.8                  0.54                  39.05
+  15                   0.0                         6.8                  0.51                  41.84
+  16                   0.0                         6.8                  0.49                  44.64
+  17                   0.0                         6.8                  0.46                  47.43
+  18                   0.0                         6.8                  0.44                  50.22
+  19                   0.0                         6.8                  0.41                  53.01
+  20                   0.0                         6.8                  0.39                  55.79
+  21                   0.0                         6.8                  0.37                  58.58
+  22                   0.0                         6.8                  0.34                  61.37
+  23                   0.0                         6.8                  0.32                  64.16
+  24                   0.0                         6.8                  0.29                  66.94
+  25                   0.0                         6.8                  0.27                  69.73
+  26                   0.0                         6.8                  0.24                  72.51
+  27                   0.0                         6.8                  0.22                  75.29
+  28                   0.0                         6.8                  0.19                  78.07
+  29                   0.0                         6.8                  0.17                  80.85
+  30                   0.0                         5.1                  0.15                  82.94
+
+
+                             ********************************
+                             *  REVENUE & CASHFLOW PROFILE  *
+                             ********************************
+Year            Electricity             |            Heat                  |           Cooling                 |         Carbon                    |          Project
+Since     Price   Ann. Rev.  Cumm. Rev. |   Price   Ann. Rev.   Cumm. Rev. |  Price   Ann. Rev.   Cumm. Rev.   |   Price   Ann. Rev.   Cumm. Rev.  | OPEX    Net Rev.      Net Cashflow
+Start    (cents/kWh)(MUSD/yr) (MUSD)    |(cents/kWh) (MUSD/yr)    (MUSD)   |(cents/kWh) (MUSD/yr)    (MUSD)    |(USD/tonne) (MUSD/yr)    (MUSD)    |(MUSD/yr) (MUSD/yr)    (MUSD)
+________________________________________________________________________________________________________________________________________________________________________________________
+  1      0.00          -38.57   0.00     |    0.00     0.00         0.00    |    0.00     0.00         0.00     |    0.00     0.00         0.00     |  0.00     -38.57     -38.57
+  2     19.00          -0.64   0.00     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -39.21
+  3     19.00          -0.64   0.01     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -39.86
+  4     19.00          -0.64   0.01     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -40.50
+  5     19.00          -0.64   0.01     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -41.14
+  6     19.00          -0.64   0.02     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -41.78
+  7     19.00          -0.64   0.02     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -42.42
+  8     19.00          -0.64   0.02     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -43.07
+  9     19.00          -0.64   0.03     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -43.71
+ 10     19.00          -0.64   0.03     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -44.35
+ 11     19.00          -0.64   0.03     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -44.99
+ 12     19.00          -0.64   0.04     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -45.63
+ 13     19.00          -0.64   0.04     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -46.28
+ 14     19.00          -0.64   0.05     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -46.92
+ 15     19.00          -0.64   0.05     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -47.56
+ 16     19.00          -0.64   0.05     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -48.20
+ 17     19.00          -0.64   0.06     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -48.84
+ 18     19.00          -0.64   0.06     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -49.48
+ 19     19.00          -0.64   0.06     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -50.13
+ 20     19.00          -0.64   0.07     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -50.77
+ 21     19.00          -0.64   0.07     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -51.41
+ 22     19.00          -0.64   0.07     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -52.05
+ 23     19.00          -0.64   0.08     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -52.69
+ 24     19.00          -0.64   0.08     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -53.34
+ 25     19.00          -0.64   0.08     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -53.98
+ 26     19.00          -0.64   0.09     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -54.62
+ 27     19.00          -0.64   0.09     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -55.26
+ 28     19.00          -0.64   0.09     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -55.90
+ 29     19.00          -0.64   0.10     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -56.55
+ 30     19.00          -0.64   0.10     |    2.50     0.00         0.00    |    2.50     0.00         0.00     |    0.00     0.00         0.00     |  0.65     -0.64     -57.19
diff --git a/tests/examples/example_SBT_Lo_T.txt b/tests/examples/example_SBT_Lo_T.txt
new file mode 100644
index 00000000..2b032c6d
--- /dev/null
+++ b/tests/examples/example_SBT_Lo_T.txt
@@ -0,0 +1,50 @@
+************reservoir*************
+Reservoir Model, 8
+Reservoir Depth, 2.4 kilometer
+Gradient 1, 31.25
+Reservoir Volume Option, 4
+Reservoir Volume, 8136407202.64, pi * 1 mi * 1 mi * 1000
+Reservoir Heat Capacity, 1112
+Reservoir Density, 2663
+Reservoir Thermal Conductivity, 2.25
+
+Lateral Endpoint Depth, 2.5 kilometer
+Lateral Inclination Angle, 89
+Junction Depth, 2.4 kilometer
+Vertical Section Length, 2.4 kilometer
+Number of Multilateral Sections, 2
+SBT Accuracy Desired, 5
+Lateral Spacing, 75
+Discretization Length, 250
+SBT Initial Timestep Count, 5
+SBT Initial to Final Timestep Transition, 10000, seconds = one year
+SBT Final Timestep Count, 120, ~4 per year
+
+****************Wellbore***********
+Is AGS, True
+Well Geometry Configuration, 5
+Number of Production Wells, 1
+Number of Injection Wells, 1
+Production Well Diameter, 8.5
+Injection Well Diameter, 8.5
+Nonvertical Wellbore Diameter, 0.216
+Production Flow Rate per Well, 6
+Reservoir Impedance, 1E-4, there is very little impedance in an AGS system
+Multilaterals Cased, False
+
+****************Surface plant************
+End-Use Option, 1
+Power Plant Type, 2
+Plant Lifetime, 30
+Ambient Temperature, 3
+Surface Temperature, 3
+Injection Temperature, 24
+
+***************Economics****************
+Economic Model, 3
+Reservoir Stimulation Capital Cost, 0
+Exploration Capital Cost, 0
+Starting Electricity Sale Price, 0.19
+Ending Electricity Sale Price, 0.19
+
+xSBT Generate Wireframe Graphics, True