Merge pull request #66 from PyPSA/fneum-patch-1

base_network: return zero underwater fraction if no link geometry given

scripts/base_network.py

@@ -411,10 +411,13 @@ def findforeignbus(G, i):
 
 def _set_links_underwater_fraction(n):
     if n.links.empty: return
-
-    offshore_shape = gpd.read_file(snakemake.input.offshore_shapes).unary_union
-    links = gpd.GeoSeries(n.links.geometry.dropna().map(shapely.wkt.loads))
-    n.links['underwater_fraction'] = links.intersection(offshore_shape).length / links.length
+
+    if not hasattr(n.links, 'geometry'):
+        n.links['underwater_fraction'] = 0.
+    else:
+        offshore_shape = gpd.read_file(snakemake.input.offshore_shapes).unary_union
+        links = gpd.GeoSeries(n.links.geometry.dropna().map(shapely.wkt.loads))
+        n.links['underwater_fraction'] = links.intersection(offshore_shape).length / links.length
 
 def _adjust_capacities_of_under_construction_branches(n):
     lines_mode = snakemake.config['lines'].get('under_construction', 'undef')