Update notebooks to produce gridcerf data v2.1

[kmongird]

This pull request:

• Updates the build_compiled_suitability.ipynb notebook: incorporates several improvements to code execution, output handling, and file management. The main changes include enabling compression for raster outputs, and simplifying directory creation logic.
• Updates the build_geothermal_resource_potential.ipynb notebook: incorporates updated geothermal potential data to build new geothermal class rasters
• Creates the build_offshore_bathymetry_suitability.ipynb notebook: Adds a new constraint layer for offshore wind related to bathymetry

[Copilot]

🟡 Changes recommended

The new offshore bathymetry notebook contains masking/metadata logic issues that can yield incorrect suitability outputs and/or inconsistent raster write metadata.

Once you've addressed the issues Copilot identified, you can request another Copilot review.

Pull request overview

This PR updates the notebook-based data build pipeline to generate GRIDCERF data v2.1 outputs, including adding an offshore bathymetry constraint layer and improving compiled raster generation.

Changes:

• Adds a new build_offshore_bathymetry_suitability.ipynb notebook to generate an offshore bathymetry exclusion raster from GEBCO data.
• Updates build_compiled_suitability.ipynb to improve scenario-layer naming, simplify output directory creation, and enable DEFLATE compression when writing compiled rasters.

File summaries

File	Description
notebooks/build_offshore_bathymetry_suitability.ipynb	New notebook to warp GEBCO bathymetry to GRIDCERF CRS/extent and produce an offshore depth suitability raster.
notebooks/build_compiled_suitability.ipynb	Improves compilation pipeline logic (layer renaming + directory creation) and writes compressed compiled rasters.

Review details

Suppressed comments (1)

notebooks/build_offshore_bathymetry_suitability.ipynb:231

• The land-mask combination logic can incorrectly flip some land pixels to "suitable" when arr is already excluded (e.g., if a land cell is NaN/previously masked): arr=1 + land marker (2) becomes 3, then arr == 2 maps it to 0. Land exclusion should be an OR with the existing exclusion mask, not an equality check on a summed value.

    "            land_array = np.where(land_array==0, 1, 2)\n",
    "            \n",
    "            # apply land area mask\n",
    "            arr += land_array\n",
    "            arr = np.where(arr ==2, 1, 0)\n",

• Files reviewed: 2/3 changed files
• Comments generated: 3
• Review effort level: Lite

---

💡 Add a code-review agent skill or configure MCP servers for context-aware, tailored reviews. Learn more in the docs.