Consider shipping routes for offshore potential calculation

Snakefile

@@ -185,6 +185,22 @@ if config['enable'].get('retrieve_natura_raster', True):
         run: move(input[0], output[0])
 
 
+if config['enable'].get('build_ship_raster', False):
+    rule build_ship_raster:
+        input:
+            ship_density="data/bundle/shipdensity/shipdensity_global.zip",
+            cutouts=expand("cutouts/{cutouts}.nc", **config['atlite'])
+        output: "resources/europe_shipdensity_raster.nc"
+        log: "logs/build_ship_raster.log"
+        script: "scripts/build_ship_raster.py"
+
+
+if config['enable'].get('retrieve_ship_raster', True):
+    rule retrieve_ship_raster:
+        input: HTTP.remote("path-to-file", keep_local=True, static=True)
+        output: "resources/europe_shipdensity_raster.nc"
+        run: move(input[0], output[0])
+
 rule build_renewable_profiles:
     input:
         base_network="networks/base.nc",
@@ -195,6 +211,9 @@ rule build_renewable_profiles:
         gebco=lambda w: ("data/bundle/GEBCO_2014_2D.nc"
                          if "max_depth" in config["renewable"][w.technology].keys()
                          else []),
+        ship_density= lambda w: ("resources/europe_shipdensity_raster.nc"
+                         if "ship_threshold" in config["renewable"][w.technology].keys()
+                         else []),
         country_shapes='resources/country_shapes.geojson',
         offshore_shapes='resources/offshore_shapes.geojson',
         regions=lambda w: ("resources/regions_onshore.geojson"

config.default.yaml

@@ -30,6 +30,8 @@ enable:
   retrieve_cutout: true
   build_natura_raster: false
   retrieve_natura_raster: true
+  build_ship_raster: false
+  retrieve_ship_raster: true
   custom_busmap: false
 
 electricity:
@@ -134,6 +136,7 @@ renewable:
     # until done more rigorously in #153
     corine: [44, 255]
     natura: true
+    ship_threshold: 400
     max_depth: 50
     max_shore_distance: 30000
     potential: simple # or conservative
@@ -151,6 +154,7 @@ renewable:
     # until done more rigorously in #153
     corine: [44, 255]
     natura: true
+    ship_threshold: 400
     max_depth: 50
     min_shore_distance: 30000
     potential: simple # or conservative

config.tutorial.yaml

@@ -31,6 +31,8 @@ enable:
   retrieve_cutout: true
   build_natura_raster: false
   retrieve_natura_raster: true
+  build_ship_raster: false
+  retrieve_ship_raster: true
   custom_busmap: false
 
 electricity:
@@ -87,6 +89,7 @@ renewable:
     # correction_factor: 0.93
     corine: [44, 255]
     natura: true
+    ship_threshold: 400
     max_shore_distance: 30000
     potential: simple # or conservative
     clip_p_max_pu: 1.e-2
@@ -100,6 +103,7 @@ renewable:
     # correction_factor: 0.93
     corine: [44, 255]
     natura: true
+    ship_threshold: 400
     min_shore_distance: 30000
     potential: simple # or conservative
     clip_p_max_pu: 1.e-2