TNO Intern

Commit fbb6e3a2 authored by Hen Brett's avatar Hen Brett 🐔
Browse files

Changing the expected value of hp_cop for the volumetric flow test following the heat pump refactor

parent 58119407
Loading
Loading
Loading
Loading
Loading
+0 −1
Original line number Diff line number Diff line
@@ -49,7 +49,6 @@ def optimize_well_distance(
            dist_min = well_distance
        else:
            dist_max = well_distance

    # If no convergence in 1000 iterations
    else:
        print(f"WARNING: Well distance optimization failed to converge. Final dist={well_distance}")

tests/utc/test_aquifer.py

deleted100644 → 0
+0 −19
Original line number Diff line number Diff line
import numpy as np
from pygridsio import read_grid,resample_grid
from pathlib import Path

test_files_in_path = Path(__file__).parent.parent / "resources" / "test_input" / "example_data"

def read_input_grids():
    new_cellsize=1000 # in m
    input_grids = resample_grid(read_grid(test_files_in_path / "ROSL_ROSLU__thick.zmap"), new_cellsize=new_cellsize).to_dataset(name="thickness_mean")
    input_grids["thickness_sd"] = resample_grid(read_grid(test_files_in_path / "ROSL_ROSLU__thick_sd.zmap"), new_cellsize=new_cellsize)
    input_grids["ntg"] = resample_grid(read_grid(test_files_in_path / "ROSL_ROSLU__ntg.zmap"), new_cellsize=new_cellsize)
    input_grids["porosity"] = resample_grid(read_grid(test_files_in_path / "ROSL_ROSLU__poro.zmap"), new_cellsize=new_cellsize) / 100
    input_grids["depth"] = resample_grid(read_grid(test_files_in_path / "ROSL_ROSLU__top.zmap"), new_cellsize=new_cellsize)
    input_grids["ln_permeability_mean"] = np.log(resample_grid(read_grid(test_files_in_path / "ROSL_ROSLU__perm.zmap"), new_cellsize=new_cellsize))
    input_grids["ln_permeability_sd"] = resample_grid(read_grid(test_files_in_path / "ROSL_ROSLU__ln_perm_sd.zmap"), new_cellsize=new_cellsize)
    return input_grids

def test_pydoubletcalc_on_aquifer():
    input_grids = read_input_grids()
+1 −1
Original line number Diff line number Diff line
@@ -21,7 +21,7 @@ def test_calculate_volumetric_flow():

    flow_results = calculate_volumetric_flow(props, input_data, 600000, 1550, 40)

    assert np.isclose(flow_results.hp_cop, 3.0, rtol=0.01)
    assert np.isclose(flow_results.hp_cop, 0.0, rtol=0.01)
    assert np.isclose(flow_results.hp_added_power, 0.0, rtol=0.01)
    assert np.isclose(flow_results.heat_power_per_doublet, 0.027926914290870505, rtol=0.01)
    assert np.isclose(flow_results.cop, 5.396068601571214, rtol=0.01)