Pymesh2d v2

methods/dynamical_downscaling/SfincsCompoundFlooding/01_Wrapper_Setup_Run.ipynb

@@ -0,0 +1,225 @@
+{
+ "cells": [
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "97b930aa",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "import geopandas as gpd\n",
+    "import numpy as np\n",
+    "import pandas as pd\n",
+    "import xarray as xr"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 7,
+   "id": "0e0f6af4",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Define grid resolution\n",
+    "resolution = 50\n",
+    "\n",
+    "# Setup grid parameters with spatial origin and grid size\n",
+    "x0, y0 = 585000, 2350004\n",
+    "x1, y1 = 632500, 2397504\n",
+    "\n",
+    "# Boundary point where ocean-driver forcing is setup\n",
+    "x, y = 601900.8911782742, 2356315.456879479\n",
+    "gdf_boundary_points = gpd.GeoDataFrame(geometry=gpd.points_from_xy([x], [y]), crs=6634)\n",
+    "\n",
+    "# Grid configuration dictionary\n",
+    "fixed_parameters = {\n",
+    "    \"path_to_dem_tif\": \"inputs/dem.tif\",\n",
+    "    \"dem_tif_var_name\": \"elevtn\",\n",
+    "    \"path_to_rgh_tif\": \"inputs/sfincs.man\",\n",
+    "    \"path_to_outflow_shp\": \"inputs/outflow_boundaries.shp\",\n",
+    "    \"path_to_waterlevel_shp\": \"inputs/inflow_boundaries.shp\",\n",
+    "    \"x0\": x0,  # Origin location of the cell edge (x0)\n",
+    "    \"y0\": y0,  # Origin location of the cell edge (y0)\n",
+    "    \"dx\": resolution,  # Grid size in x-direction\n",
+    "    \"dy\": resolution,  # Grid size in y-direction\n",
+    "    \"mmax\": (x1 - x0) // resolution,  # Number of grid cells in x-direction\n",
+    "    \"nmax\": (y1 - y0) // resolution,  # Number of grid cells in y-direction\n",
+    "    \"rotation\": 0,  # Rotation angle in degrees (anti-clockwise from east)\n",
+    "    \"epsg\": 6634,  # Coordinate reference system (CRS) EPSG code\n",
+    "    \"path_to_mask\": \"inputs/mask.shp\",\n",
+    "    \"nr_subgrid_pixels\": 1,\n",
+    "    \"gdf_boundary_points\": gdf_boundary_points,\n",
+    "}"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "id": "4c323def",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Creating time chunks as individual cases\n",
+    "K = 5"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "d5523cb9",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Waterlevels\n",
+    "waterlevels_forcing = pd.read_pickle(\"inputs/waterlevels.pkl\")\n",
+    "\n",
+    "time_chunks = np.array_split(waterlevels_forcing.index.sort_values().unique(), K)\n",
+    "time_bounds = [(edge[0], edge[-1]) for edge in time_chunks]\n",
+    "\n",
+    "waterlevels_forcing = np.array_split(waterlevels_forcing, K)\n",
+    "\n",
+    "slowly_waterlevel = [\n",
+    "    pd.DataFrame(i[\"Msetup\"]).set_axis([0], axis=1) for i in waterlevels_forcing\n",
+    "]\n",
+    "quickly_waterlevel = [\n",
+    "    pd.DataFrame(i[\"Hig\"]).set_axis([0], axis=1) for i in waterlevels_forcing\n",
+    "]"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "592f477d",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "# Precipitation\n",
+    "precipitation_forcing = xr.open_dataset(\"inputs/precipitation.nc\")\n",
+    "precipitation_forcing = precipitation_forcing.rio.write_crs(\"epsg:4326\")\n",
+    "precipitation_forcing = precipitation_forcing.sortby(\"time\")\n",
+    "\n",
+    "precipitation_forcing = [\n",
+    "    precipitation_forcing.sel(time=slice(str(start), str(end)))\n",
+    "    for start, end in time_bounds\n",
+    "]"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "0855f806",
+   "metadata": {},
+   "outputs": [],
+   "source": [
+    "metamodel_parameters = {\n",
+    "    \"precipitation_forcing\": precipitation_forcing,\n",
+    "    \"slowly_waterlevel_forcing\": slowly_waterlevel,\n",
+    "    \"quickly_waterlevel_forcing\": quickly_waterlevel,\n",
+    "}"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "12c34008",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stderr",
+     "output_type": "stream",
+     "text": [
+      "2025-10-02 03:40:51,291 - SfincsAlbaModelWrapper - WARNING - Parameter precipitation_forcing is not in the default_parameters\n",
+      "2025-10-02 03:40:51,293 - SfincsAlbaModelWrapper - WARNING - Parameter slowly_waterlevel_forcing is not in the default_parameters\n",
+      "2025-10-02 03:40:51,339 - SfincsAlbaModelWrapper - WARNING - Parameter quickly_waterlevel_forcing is not in the default_parameters\n"
+     ]
+    }
+   ],
+   "source": [
+    "from bluemath_tk.wrappers.sfincs.sfincs_wrapper import SfincsModelWrapper\n",
+    "\n",
+    "sfincs_model = SfincsModelWrapper(\n",
+    "    templates_dir=\"templates\",\n",
+    "    metamodel_parameters=metamodel_parameters,\n",
+    "    fixed_parameters=fixed_parameters,\n",
+    "    output_dir=\"CASES\",\n",
+    ")"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": 8,
+   "id": "b83fb3a2",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stderr",
+     "output_type": "stream",
+     "text": [
+      "Model dir already exists and files might be overwritten: /home/ibbisin0/globus/ricondal/Bluemath_Github/BlueMath/methods/dynamical_downscaling/SfincsExample/Alba_CASES/0000/gis.\n"
+     ]
+    },
+    {
+     "name": "stderr",
+     "output_type": "stream",
+     "text": [
+      "Model dir already exists and files might be overwritten: /home/ibbisin0/globus/ricondal/Bluemath_Github/BlueMath/methods/dynamical_downscaling/SfincsExample/Alba_CASES/0001/gis.\n",
+      "Model dir already exists and files might be overwritten: /home/ibbisin0/globus/ricondal/Bluemath_Github/BlueMath/methods/dynamical_downscaling/SfincsExample/Alba_CASES/0002/gis.\n",
+      "Model dir already exists and files might be overwritten: /home/ibbisin0/globus/ricondal/Bluemath_Github/BlueMath/methods/dynamical_downscaling/SfincsExample/Alba_CASES/0003/gis.\n",
+      "Model dir already exists and files might be overwritten: /home/ibbisin0/globus/ricondal/Bluemath_Github/BlueMath/methods/dynamical_downscaling/SfincsExample/Alba_CASES/0004/gis.\n"
+     ]
+    }
+   ],
+   "source": [
+    "sfincs_model.build_cases()"
+   ]
+  },
+  {
+   "cell_type": "code",
+   "execution_count": null,
+   "id": "f08271ab",
+   "metadata": {},
+   "outputs": [
+    {
+     "name": "stderr",
+     "output_type": "stream",
+     "text": [
+      "2025-10-02 03:18:10,417 - SfincsAlbaModelWrapper - ERROR - Error running command: docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu\n",
+      "2025-10-02 03:18:10,418 - SfincsAlbaModelWrapper - ERROR - Error: Command 'docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu' returned non-zero exit status 2.\n",
+      "2025-10-02 03:18:13,839 - SfincsAlbaModelWrapper - ERROR - Error running command: docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu\n",
+      "2025-10-02 03:18:13,840 - SfincsAlbaModelWrapper - ERROR - Error: Command 'docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu' returned non-zero exit status 2.\n",
+      "2025-10-02 03:18:16,976 - SfincsAlbaModelWrapper - ERROR - Error running command: docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu\n",
+      "2025-10-02 03:18:16,976 - SfincsAlbaModelWrapper - ERROR - Error: Command 'docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu' returned non-zero exit status 2.\n",
+      "2025-10-02 03:18:20,365 - SfincsAlbaModelWrapper - ERROR - Error running command: docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu\n",
+      "2025-10-02 03:18:20,365 - SfincsAlbaModelWrapper - ERROR - Error: Command 'docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu' returned non-zero exit status 2.\n",
+      "2025-10-02 03:18:23,651 - SfincsAlbaModelWrapper - ERROR - Error running command: docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu\n",
+      "2025-10-02 03:18:23,651 - SfincsAlbaModelWrapper - ERROR - Error: Command 'docker run --rm -v .:/case_dir -w /case_dir deltares/sfincs-cpu' returned non-zero exit status 2.\n"
+     ]
+    }
+   ],
+   "source": [
+    "sfincs_model.run_cases(launcher=\"docker\")"
+   ]
+  }
+ ],
+ "metadata": {
+  "kernelspec": {
+   "display_name": "bluemath-dev",
+   "language": "python",
+   "name": "python3"
+  },
+  "language_info": {
+   "codemirror_mode": {
+    "name": "ipython",
+    "version": 3
+   },
+   "file_extension": ".py",
+   "mimetype": "text/x-python",
+   "name": "python",
+   "nbconvert_exporter": "python",
+   "pygments_lexer": "ipython3",
+   "version": "3.11.11"
+  }
+ },
+ "nbformat": 4,
+ "nbformat_minor": 5
+}