depth:float=typer.Option(500,help="The depth of the aquifer, +ive downwards, units: [m]"),
temperature:float=typer.Option(None,help="The temperature of the aquifer, if not provided this will be calculated using a Temperature gradient of 8°C + 31°C/km with the depth parameter,"),
thickness:float=typer.Option(500,help="The thickness of the aquifer, units: [m]"),
ntg:float=typer.Option(0.5,help="The net-to-gross of the aquifer, units: [0-1]"),
porosity:float=typer.Option(0.5,help="The porosity of the aquifer, between 0-1 (1 = 100%), units: [0-1]"),
permeability:float=typer.Option(0.5,help="The permeability of the aquifer, units: [mD]"),
output_file:Path=typer.Option(None,help="A file to output the data to, accepted filetypes are: .nc, .csv, .xlsx, .json"),
verbose:bool=typer.Option(False,help="print the input data and output data to terminal")
)->None:
"""Simulate a Geothermal Doublet, with no stochastics (only a P50 simulation)"""
# instantiate the input_data dataset from the input
Perform a ThermoGIS Doublet performance simulation. This will occur across all dimensions of the input_data (ie. input data can have a single value for each required variable, or it can be 1Dimensional or a 2Dimensional grid)
@@ -10,6 +10,7 @@ def generate_thickness_permeability_transmissivity_for_pvalues(thickness_mean: f
-transmissivity
Note: Transmissivity is a function of ln(permeability) * thickness, and so no analytical solution exists to combine these two probability distributions and so the transmissivity distribution has to be sampled.
UNLESS pvalue = P50; then the analytical solution for Transmissivity = (np.log(mean_thickness) + ln_permeability_mean)
:param thickness_mean:
:param thickness_sd:
@@ -22,18 +23,26 @@ def generate_thickness_permeability_transmissivity_for_pvalues(thickness_mean: f
