diff --git a/notebooks/entsoe_data.ipynb.license b/notebooks/entsoe_data.ipynb.license
new file mode 100644
index 00000000..861897bd
--- /dev/null
+++ b/notebooks/entsoe_data.ipynb.license
@@ -0,0 +1,3 @@
+SPDX-FileCopyrightText: 2021 - 2023 The PyPSA-EUR Authors
+
+SPDX-License-Identifier: MIT
diff --git a/scripts/build_natura_raster.py b/scripts/build_natura_raster.py
index 3cd62fd9..8fdb4ea3 100644
--- a/scripts/build_natura_raster.py
+++ b/scripts/build_natura_raster.py
@@ -54,6 +54,23 @@ logger = logging.getLogger(__name__)
 
 
 def determine_cutout_xXyY(cutout_name):
+    """
+    Determine the full extent of a cutout.
+
+    Since the coordinates of the cutout data are given as the
+    center of the grid cells, the extent of the cutout is
+    calculated by adding/subtracting half of the grid cell size.
+
+
+    Parameters
+    ----------
+    cutout_name : str
+        Path to the cutout.
+
+    Returns
+    -------
+    A list of extent coordinates in the order [x, X, y, Y].
+    """
     cutout = atlite.Cutout(cutout_name)
     assert cutout.crs.to_epsg() == 4326
     x, X, y, Y = cutout.extent