Result comparison

[visr]

The integration tests compare the output states of the HWS model. For #2136 I needed a quick way to eyeball the differences on the map. Sharing my throw-away script here before throwing it away such that I don't have to rewrite it later.

import ribasim
import pandas as pd
import geopandas as gpd

new_path = r"hws_2024_7_0\results\basin_state.arrow"
old_path = r"hws_2024_7_0\benchmark\basin_state.arrow"
gpkg = r"hws_2024_7_0\database.gpkg"

# compute the differences
new = pd.read_feather(new).set_index("node_id")
old = pd.read_feather(old).set_index("node_id")
dif = new - old

# load the spatial Node layer and add the differences
node = gpd.read_file(gpkg, layer="Node", fid_as_index=True)
basin = node[node.node_type == "Basin"]
basin["dif"] = dif
basin["new"] = new
basin["old"] = old

basin.to_file(
    r"hws_2024_7_0\results\basin_dif.gpkg",
    driver="GPKG",
    index=True,
)