Incorporate into the shape files
data_path_panden = get_abs_data_path(
name="oppervlaktewater_pwn_shapes_panden",
version="latest",
location="get_from_env",
)
panden_shp = gpd.read_file(os.path.join(data_path_panden, "Panden_ICAS_IKIEF.shp"))
panden_shp["c"] = 1.0 # bodemweerstand
panden_shp["stage"] = np.nan
panden_shp.loc[panden_shp.Naam.str.contains("ICAS"), "stage"] = 2.8 # mNAP
panden_shp.loc[panden_shp.Naam.str.contains("IKIEF"), "stage"] = 5.8 # mNAP
panden_shp["rbot"] = panden_shp["stage"] - 2.0