Skip to content

Commit

Permalink
Output the Gradient results in the preferred units
Browse files Browse the repository at this point in the history
  • Loading branch information
@natakokota
natakokota committed Jun 9, 2024
1 parent becec79 commit cff4408
Showing 1 changed file with 33 additions and 19 deletions.
52 changes: 33 additions & 19 deletions src/geophires_x/Outputs.py
View file
Original file line number Diff line number Diff line change
Expand Up @@ -729,6 +729,27 @@ def read_parameters(self, model:Model) -> None:
model.logger.info(f'Complete {__class__!s}: {__name__}')

def PrintOutputs(self, model: Model):
def gradient_output(j):
if model.reserv.gradient.PreferredUnits.value != 'degC/m':
return '{0:10.4f}'.format(model.reserv.gradient.value[j] * 1000.0), model.reserv.gradient.PreferredUnits.value
else:
return '{0:10.4f}'.format(model.reserv.gradient.value[j]), model.reserv.gradient.CurrentUnits.value
def gradient_output2():
if model.reserv.gradient.PreferredUnits.value != 'degC/m':
return f' Geothermal gradient: {model.reserv.gradient.value[0] * 1000.0:10.4f} ' + model.reserv.gradient.PreferredUnits.value + NL
else:
return f' Geothermal gradient: {model.reserv.gradient.value[0]:10.4f} ' + model.reserv.gradient.CurrentUnits.value + NL
def gradient_output3():
if model.reserv.gradient.PreferredUnits.value != 'degC/m':
return f' Segment {str(i):s} Geothermal gradient: {model.reserv.gradient.value[i-1] * 1000.0:10.4f} ' + model.reserv.gradient.PreferredUnits.value +NL
else:
return f' Segment {str(i):s} Geothermal gradient: {model.reserv.gradient.value[i-1]:10.4f} ' + model.reserv.gradient.CurrentUnits.value +NL
def gradient_output4():
if model.reserv.gradient.PreferredUnits.value != 'degC/m':
return f' Segment {str(i+1):s} Geothermal gradient: {model.reserv.gradient.value[i] * 1000.0:10.4f} ' + model.reserv.gradient.PreferredUnits.value + NL
else:
return f' Segment {str(i+1):s} Geothermal gradient: {model.reserv.gradient.value[i]:10.4f} ' + model.reserv.gradient.CurrentUnits.value + NL

"""
PrintOutputs writes the standard outputs to the output file.
:param model: The container class of the application, giving access to everything else, including the logger
Expand Down Expand Up @@ -847,19 +868,16 @@ def PrintOutputs(self, model: Model):
model.reserv.depth.CurrentUnits.value))

if model.reserv.numseg.value == 1:
summary.append(OutputTableItem('Geothermal gradient', '{0:10.4f}'.format(model.reserv.gradient.value[0]),
model.reserv.gradient.CurrentUnits.value))
summary.append(OutputTableItem('Geothermal gradient', gradient_output(0)))
else:
for i in range(1, model.reserv.numseg.value):
summary.append(OutputTableItem(f'Segment {str(i)} Geothermal gradient',
'{0:10.4f}'.format(model.reserv.gradient.value[i - 1]),
model.reserv.gradient.CurrentUnits.value))
gradient_output(i-1)))
summary.append(OutputTableItem(f'Segment {str(i)} Thickness',
'{0:10.0f}'.format(model.reserv.layerthickness.value[i - 1]),
model.reserv.layerthickness.CurrentUnits.value))
summary.append(OutputTableItem(f'Segment {str(i + 1)} Geothermal gradient',
'{0:10.4f}'.format(model.reserv.gradient.value[i]),
model.reserv.gradient.CurrentUnits.value))
gradient_output(i)))
if model.economics.DoCarbonCalculations.value:
summary.append(OutputTableItem('Total Avoided Carbon Emissions', '{0:10.2f}'.format(
model.economics.CarbonThatWouldHaveBeenProducedTotal.value * 0.000453592), 'metric tonnes'))
Expand Down Expand Up @@ -946,19 +964,15 @@ def PrintOutputs(self, model: Model):
model.reserv.Tmax.CurrentUnits.value))
resource_characteristics.append(OutputTableItem('Number of segments', '{0:10.0f}'.format(model.reserv.numseg.value)))
if model.reserv.numseg.value == 1:
resource_characteristics.append(OutputTableItem('Geothermal gradient', '{0:10.4f}'.format(model.reserv.gradient.value[0]),
model.reserv.gradient.CurrentUnits.value))
resource_characteristics.append(OutputTableItem('Geothermal gradient', gradient_output(0)))
else:
for i in range(1, model.reserv.numseg.value):
resource_characteristics.append(OutputTableItem(f'Segment {str(i)} Geothermal gradient',
'{0:10.4f}'.format(
model.reserv.gradient.value[i - 1]),
model.reserv.gradient.CurrentUnits.value))
gradient_output(i-1)))
resource_characteristics.append(OutputTableItem(f'Segment {str(i)} Thickness', '{0:10.0f}'.format(
model.reserv.layerthickness.value[i - 1]), model.reserv.layerthickness.CurrentUnits.value))
resource_characteristics.append(OutputTableItem(f'Segment {str(i + 1)} Geothermal gradient',
'{0:10.4f}'.format(model.reserv.gradient.value[i]),
model.reserv.gradient.CurrentUnits.value))
gradient_output(i)))
if model.wellbores.IsAGS.value:
reservoir_parameters.append(OutputTableItem('The AGS models contain an intrinsic reservoir model that doesn\'t expose values that can be used in extensive reporting.'))
else:
Expand Down Expand Up @@ -1596,12 +1610,12 @@ def PrintOutputs(self, model: Model):
f.write(f' Well depth (or total length, if not vertical): {model.reserv.depth.value:10.1f} ' +model.reserv.depth.CurrentUnits.value + NL)

if model.reserv.numseg.value == 1:
f.write(f' Geothermal gradient: {model.reserv.gradient.value[0]:10.4f} ' + model.reserv.gradient.CurrentUnits.value + NL)
f.write(gradient_output2())
else:
for i in range(1, model.reserv.numseg.value):
f.write(f' Segment {str(i):s} Geothermal gradient: {model.reserv.gradient.value[i-1]:10.4f} ' + model.reserv.gradient.CurrentUnits.value +NL)
f.write(gradient_output3())
f.write(f' Segment {str(i):s} Thickness: {model.reserv.layerthickness.value[i-1]:10.0f} ' + model.reserv.layerthickness.CurrentUnits.value + NL)
f.write(f' Segment {str(i+1):s} Geothermal gradient: {model.reserv.gradient.value[i]:10.4f} ' + model.reserv.gradient.CurrentUnits.value + NL)
f.write(gradient_output4())
if model.economics.DoCarbonCalculations.value:
f.write(f' Total Avoided Carbon Emissions: {model.economics.CarbonThatWouldHaveBeenProducedTotal.value*0.000453592:10.2f} metric tonnes' + NL)

Expand Down Expand Up @@ -1666,12 +1680,12 @@ def PrintOutputs(self, model: Model):
f.write(f' Maximum reservoir temperature: {model.reserv.Tmax.value:10.1f} ' + model.reserv.Tmax.CurrentUnits.value + NL)
f.write(f' Number of segments: {model.reserv.numseg.value:10.0f} ' + NL)
if model.reserv.numseg.value == 1:
f.write(f' Geothermal gradient: {model.reserv.gradient.value[0]:10.4f} ' + model.reserv.gradient.CurrentUnits.value + NL)
f.write(gradient_output2())
else:
for i in range(1, model.reserv.numseg.value):
f.write(f' Segment {str(i):s} Geothermal gradient: {model.reserv.gradient.value[i-1]:10.4f} ' + model.reserv.gradient.CurrentUnits.value +NL)
f.write(gradient_output3())
f.write(f' Segment {str(i):s} Thickness: {model.reserv.layerthickness.value[i-1]:10.0f} ' + model.reserv.layerthickness.CurrentUnits.value + NL)
f.write(f' Segment {str(i+1):s} Geothermal gradient: {model.reserv.gradient.value[i]:10.4f} ' + model.reserv.gradient.CurrentUnits.value + NL)
f.write(gradient_output4())

f.write(NL)
f.write(NL)
Expand Down

0 comments on commit cff4408

Please sign in to comment.