Fix: Invalid aquifer shapes in build_ates_potentials (PyPSA#1724)

* fix: ensure aquifer geometries are valid in build_ates_potentials.py

* fix: add release note for fixing invalid aquifer shape geometries in build_ates_potentials.py

Author: amos-schledorn

doc/release_notes.rst

@@ -11,6 +11,8 @@ Release Notes
 .. Upcoming Release
 .. ================
 
+* Fix: Invalid aquifer shape geometries are now fixed in `build_ates_potentials.py` (fixing https://github.com/PyPSA/pypsa-eur/issues/1696)
+
 * Fix: Sanitize columns in `add_brownfield` as it's done for `add_exisiting_baseyear` (https://github.com/PyPSA/pypsa-eur/pull/1676).
 
 * (Breaking) Consolidate gap-filling strategies options under a new configuration section `load:fill_gaps` and add a switch (https://github.com/PyPSA/pypsa-eur/pull/1677). The options `load:interpolate_limit` and `load:time_shift_for_large_gaps` are now located under `load:fill_gaps` as `load:fill_gaps:interpolate_limit` and `load:fill_gaps:time_shift_for_large_gaps`.

scripts/build_ates_potentials.py

@@ -391,6 +391,8 @@ def check_aquifer_coverage(
     aquifer_shapes = gpd.read_file(snakemake.input.aquifer_shapes_shp).to_crs(
         regions_onshore.crs
     )
+    # fix any invalid geometries
+    aquifer_shapes.geometry = aquifer_shapes.geometry.make_valid()
 
     dh_areas = gpd.read_file(snakemake.input.dh_areas).to_crs(regions_onshore.crs)