Implement climate projections

Snakefile

@@ -57,6 +57,8 @@ wildcard_constraints:
     ll="(v|c)([0-9\.]+|opt|all)|all",
     opts="[-+a-zA-Z0-9\.]*",
     unc="[-+a-zA-Z0-9\.]*",
+    # "^proj-" is reserved for projections, negative lookbehind assertion used in the regexp
+    cutout=".+(?<!proj-)",
 
 
 if config["custom_rules"] is not []:
@@ -363,6 +365,33 @@ if config["enable"].get("build_cutout", False):
             "scripts/build_cutout.py"
 
 
+if config["enable"].get("climate_projection_cutout", False):
+
+    rule build_climate_projections:
+        params:
+            snapshots=config["snapshots"],
+            climate_scenario=config["projection"]["climate_scenario"],
+            base_year=config["projection"]["base_year"],
+            future_year=config["projection"]["future_year"],
+            years_window=config["projection"]["years_window"],
+            cmip6_nn_fl=config["projection"]["cmip6_nn_fl"],
+            cmip6_xx_fl=config["projection"]["cmip6_xx_fl"],
+        input:
+            cutout="cutouts/" + CDIR + "{cutout}.nc",
+            cmip6_avr="data/cmip6/t_CMIP6_ssp245_mon_201501-210012.nc",
+        output:
+            "cutouts/" + CDIR + "proj-{cutout}.nc",
+        log:
+            "logs/" + RDIR + "build_climate_projections/proj-{cutout}.log",
+        benchmark:
+            "benchmarks/" + RDIR + "proj-{cutout}"
+        threads: ATLITE_NPROCESSES
+        resources:
+            mem_mb=ATLITE_NPROCESSES * 1000,
+        script:
+            "scripts/build_climate_projections.py"
+
+
 if config["enable"].get("build_natura_raster", False):
 
     rule build_natura_raster:

config.default.yaml

@@ -21,6 +21,7 @@ enable:
   # If "build_cutout" : true, then environmental data is extracted according to `snapshots` date range and `countries`
   # requires cds API key https://cds.climate.copernicus.eu/api-how-to
   # More information https://atlite.readthedocs.io/en/latest/introduction.html#datasets
+  climate_projection_cutout: false
 
 custom_rules: []  # Default empty [] or link to custom rule file e.g. ["my_folder/my_rules.smk"] that add rules to Snakefile
 
@@ -203,6 +204,16 @@ atlite:
       # y: [33., 72]    # manual set cutout range
 
 
+projection:
+  climate_scenario: ssp245
+  base_year: 2020
+  future_year: 2050
+  years_window: 5 # a half of width currently
+  gcm_selection: false # false to use all the available global climate models
+  cmip6_nn_fl: data/cmip6/tnn_CMIP6_ssp245_mon_201501-210012.nc # false or names of min-extremal files of climate inputs
+  cmip6_xx_fl: data/cmip6/txx_CMIP6_ssp245_mon_201501-210012.nc # false or names of max-extremal files of climate inputs
+
+
 renewable:
   onwind:
     cutout: cutout-2013-era5

config.tutorial.yaml

@@ -30,6 +30,7 @@ enable:
   # If "build_cutout" : true # requires cds API key https://cds.climate.copernicus.eu/api-how-to
   # More information https://atlite.readthedocs.io/en/latest/introduction.html#datasets
   build_cutout: false
+  climate_projection_cutout: false
   build_natura_raster: true  # If True, then build_natura_raster can be run
 
 custom_rules: []  # Default empty [] or link to custom rule file e.g. ["my_folder/my_rules.smk"] that add rules to Snakefile
@@ -217,6 +218,17 @@ atlite:
       # x: [-12., 35.]  # set cutout range manual, instead of automatic by boundaries of country
       # y: [33., 72]    # manual set cutout range
 
+
+projection:
+  climate_scenario: ssp245
+  base_year: 2020
+  future_year: 2050
+  years_window: 5 # a half of width currently
+  gcm_selection: false # false to use all the available global climate models
+  cmip6_nn_fl: data/cmip6/tnn_CMIP6_ssp245_mon_201501-210012.nc # false or names of min-extremal files of climate inputs
+  cmip6_xx_fl: data/cmip6/txx_CMIP6_ssp245_mon_201501-210012.nc # false or names of max-extremal files of climate inputs
+
+
 renewable:
   onwind:
     cutout: cutout-2013-era5-tutorial