diff --git a/Scenario_wolfgpu/Q10-5-5/__bath_assembly.vrt b/Scenario_wolfgpu/Q10-5-5/__bath_assembly.vrt
new file mode 100644
index 0000000000000000000000000000000000000000..08279a1c3c9190fb659aacb5f48c35082545ce83
--- /dev/null
+++ b/Scenario_wolfgpu/Q10-5-5/__bath_assembly.vrt
@@ -0,0 +1,16 @@
+<VRTDataset rasterXSize="480" rasterYSize="1160">
+  <SRS dataAxisToSRSAxisMapping="1,2">PROJCS["BD72 / Belgian Lambert 72",GEOGCS["BD72",DATUM["Reseau_National_Belge_1972",SPHEROID["International 1924",6378388,297,AUTHORITY["EPSG","7022"]],AUTHORITY["EPSG","6313"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4313"]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["latitude_of_origin",90],PARAMETER["central_meridian",4.36748666666667],PARAMETER["standard_parallel_1",51.1666672333333],PARAMETER["standard_parallel_2",49.8333339],PARAMETER["false_easting",150000.013],PARAMETER["false_northing",5400088.438],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","31370"]]</SRS>
+  <GeoTransform>  2.5100000000000000e+05,  5.0000000000000000e+00,  0.0000000000000000e+00,  1.4130000000000000e+05,  0.0000000000000000e+00, -5.0000000000000000e+00</GeoTransform>
+  <VRTRasterBand dataType="Float32" band="1">
+    <NoDataValue>99999</NoDataValue>
+    <ColorInterp>Gray</ColorInterp>
+    <ComplexSource resampling="bilinear">
+      <SourceFilename relativeToVRT="1">..\bathymetry.tif</SourceFilename>
+      <SourceBand>1</SourceBand>
+      <SourceProperties RasterXSize="480" RasterYSize="1160" DataType="Float32" BlockXSize="480" BlockYSize="4" />
+      <SrcRect xOff="0" yOff="0" xSize="480" ySize="1160" />
+      <DstRect xOff="0" yOff="0" xSize="480" ySize="1160" />
+      <NODATA>99999</NODATA>
+    </ComplexSource>
+  </VRTRasterBand>
+</VRTDataset>
diff --git a/Scenario_wolfgpu/Q10-5-5/__bathymetry.tif b/Scenario_wolfgpu/Q10-5-5/__bathymetry.tif
new file mode 100644
index 0000000000000000000000000000000000000000..80698e8f53a1ed1edde5d7821cecbbddb9113569
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/__bathymetry.tif differ
diff --git a/Scenario_wolfgpu/Q10-5-5/__bathymetry_after_scripts.tif b/Scenario_wolfgpu/Q10-5-5/__bathymetry_after_scripts.tif
new file mode 100644
index 0000000000000000000000000000000000000000..71a1c41ca9b2147ed7d9d36e5f64e62205880b7f
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/__bathymetry_after_scripts.tif differ
diff --git a/Scenario_wolfgpu/Q10-5-5/__infiltration_after_scripts.tif b/Scenario_wolfgpu/Q10-5-5/__infiltration_after_scripts.tif
new file mode 100644
index 0000000000000000000000000000000000000000..af886f11e3cd339b8d47bb61d11d83c601cd27b2
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/__infiltration_after_scripts.tif differ
diff --git a/Scenario_wolfgpu/Q10-5-5/__mann_assembly.vrt b/Scenario_wolfgpu/Q10-5-5/__mann_assembly.vrt
new file mode 100644
index 0000000000000000000000000000000000000000..74a8e50ea9faa86ecd6b61143df48cf54de7c165
--- /dev/null
+++ b/Scenario_wolfgpu/Q10-5-5/__mann_assembly.vrt
@@ -0,0 +1,16 @@
+<VRTDataset rasterXSize="480" rasterYSize="1160">
+  <SRS dataAxisToSRSAxisMapping="1,2">PROJCS["BD72 / Belgian Lambert 72",GEOGCS["BD72",DATUM["Reseau_National_Belge_1972",SPHEROID["International 1924",6378388,297,AUTHORITY["EPSG","7022"]],AUTHORITY["EPSG","6313"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4313"]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["latitude_of_origin",90],PARAMETER["central_meridian",4.36748666666667],PARAMETER["standard_parallel_1",51.1666672333333],PARAMETER["standard_parallel_2",49.8333339],PARAMETER["false_easting",150000.013],PARAMETER["false_northing",5400088.438],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","31370"]]</SRS>
+  <GeoTransform>  2.5100000000000000e+05,  5.0000000000000000e+00,  0.0000000000000000e+00,  1.4130000000000000e+05,  0.0000000000000000e+00, -5.0000000000000000e+00</GeoTransform>
+  <VRTRasterBand dataType="Float32" band="1">
+    <NoDataValue>99999</NoDataValue>
+    <ColorInterp>Gray</ColorInterp>
+    <ComplexSource resampling="bilinear">
+      <SourceFilename relativeToVRT="1">..\manning.tif</SourceFilename>
+      <SourceBand>1</SourceBand>
+      <SourceProperties RasterXSize="480" RasterYSize="1160" DataType="Float32" BlockXSize="480" BlockYSize="4" />
+      <SrcRect xOff="0" yOff="0" xSize="480" ySize="1160" />
+      <DstRect xOff="0" yOff="0" xSize="480" ySize="1160" />
+      <NODATA>99999</NODATA>
+    </ComplexSource>
+  </VRTRasterBand>
+</VRTDataset>
diff --git a/Scenario_wolfgpu/Q10-5-5/__manning.tif b/Scenario_wolfgpu/Q10-5-5/__manning.tif
new file mode 100644
index 0000000000000000000000000000000000000000..13a22b2c8dfbaf39f6081a2e098d080f6a0646f4
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/__manning.tif differ
diff --git a/Scenario_wolfgpu/Q10-5-5/__manning_after_scripts.tif b/Scenario_wolfgpu/Q10-5-5/__manning_after_scripts.tif
new file mode 100644
index 0000000000000000000000000000000000000000..b72534d15572482a1441fc969b3b7fca70db0882
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/__manning_after_scripts.tif differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/NAP.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/NAP.npy
new file mode 100644
index 0000000000000000000000000000000000000000..4d7f4d8355c6692d860b9aa69abdc03adac17531
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/NAP.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/bathymetry.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/bathymetry.npy
new file mode 100644
index 0000000000000000000000000000000000000000..125bb3bdc40899c3fe6fef1965ad724a6b89f712
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/bathymetry.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/h.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/h.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/h.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/infiltration_zones.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/infiltration_zones.npy
new file mode 100644
index 0000000000000000000000000000000000000000..9e40e78cb4a2fcb263ce31074a12d9fcb99dd446
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/infiltration_zones.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/manning.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/manning.npy
new file mode 100644
index 0000000000000000000000000000000000000000..ede955942e91c3c995ff5fbbca1f42b9c3a544ec
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/manning.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/parameters.json b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/parameters.json
new file mode 100644
index 0000000000000000000000000000000000000000..d2633f4e441c7b1f16527dfa2d9dc49062b79812
--- /dev/null
+++ b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/parameters.json
@@ -0,0 +1,130 @@
+{
+  "spec_version": "0.1",
+  "comment": null,
+  "creation_date": "21/12/2023 14:31:45",
+  "parameters": {
+    "courant": 0.4,
+    "froude_max": 3,
+    "nx": 480,
+    "ny": 1160,
+    "dx": 5.0,
+    "dy": 5.0,
+    "timestep_strategy": "OPTIMIZED_TIME_STEP",
+    "timestep": 0.0,
+    "infiltration_interp": "NONE",
+    "duration": {
+      "type": "SECONDS",
+      "duration": 100000.0
+    },
+    "report_period": {
+      "type": "SECONDS",
+      "duration": 3600.0
+    },
+    "runge_kutta": 0.5,
+    "base_coord_x": 251000.0,
+    "base_coord_y": 135500.0
+  },
+  "boundary_conditions": [
+    {
+      "i": 3,
+      "j": 997,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 998,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 999,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1000,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1001,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1002,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    }
+  ],
+  "infiltrations": [
+    {
+      "start_time": 0.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 3600.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 7200.0,
+      "q": [
+        20.0,
+        10.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 43200.0,
+      "q": [
+        100.0,
+        50.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 86400.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 100000.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    }
+  ],
+  "maps": {
+    "h": "h.npy",
+    "qx": "qx.npy",
+    "qy": "qy.npy",
+    "bathymetry": "bathymetry.npy",
+    "manning": "manning.npy",
+    "infiltration_zones": "infiltration_zones.npy",
+    "NAP": "NAP.npy"
+  }
+}
\ No newline at end of file
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/quickrun.py b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/quickrun.py
new file mode 100644
index 0000000000000000000000000000000000000000..2b6352abd0f623e3f8422b1a08d0784afb17ff02
--- /dev/null
+++ b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/quickrun.py
@@ -0,0 +1,11 @@
+from pathlib import Path
+from wolfgpu.simple_simulation import SimpleSimulation
+from wolfgpu.SimulationRunner import SimulationRunner
+from wolfgpu.glsimulation import ResultsStore
+
+def main():
+	sim = SimpleSimulation.load(Path(__file__).parent)
+	result_store:ResultsStore = SimulationRunner.quick_run(sim, Path(__file__).parent)
+
+if __name__ == '__main__':
+	main()
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/qx.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/qx.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/qx.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/qy.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/qy.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_Mono25/qy.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/NAP.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/NAP.npy
new file mode 100644
index 0000000000000000000000000000000000000000..4d7f4d8355c6692d860b9aa69abdc03adac17531
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/NAP.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/bathymetry.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/bathymetry.npy
new file mode 100644
index 0000000000000000000000000000000000000000..125bb3bdc40899c3fe6fef1965ad724a6b89f712
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/bathymetry.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/h.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/h.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/h.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/infiltration_zones.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/infiltration_zones.npy
new file mode 100644
index 0000000000000000000000000000000000000000..9e40e78cb4a2fcb263ce31074a12d9fcb99dd446
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/infiltration_zones.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/manning.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/manning.npy
new file mode 100644
index 0000000000000000000000000000000000000000..ede955942e91c3c995ff5fbbca1f42b9c3a544ec
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/manning.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/parameters.json b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/parameters.json
new file mode 100644
index 0000000000000000000000000000000000000000..f8aca177253ac7271aa147eaa0ac308cc0ce8f1b
--- /dev/null
+++ b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/parameters.json
@@ -0,0 +1,98 @@
+{
+  "spec_version": "0.1",
+  "comment": null,
+  "creation_date": "21/12/2023 14:31:45",
+  "parameters": {
+    "courant": 0.4,
+    "froude_max": 3,
+    "nx": 480,
+    "ny": 1160,
+    "dx": 5.0,
+    "dy": 5.0,
+    "timestep_strategy": "OPTIMIZED_TIME_STEP",
+    "timestep": 0.0,
+    "infiltration_interp": "NONE",
+    "duration": {
+      "type": "SECONDS",
+      "duration": 86400.0
+    },
+    "report_period": {
+      "type": "SECONDS",
+      "duration": 3600.0
+    },
+    "runge_kutta": 0.5,
+    "base_coord_x": 251000.0,
+    "base_coord_y": 135500.0
+  },
+  "boundary_conditions": [
+    {
+      "i": 3,
+      "j": 997,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 998,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 999,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1000,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1001,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1002,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    }
+  ],
+  "infiltrations": [
+    {
+      "start_time": 0.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 86400.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    }
+  ],
+  "maps": {
+    "h": "h.npy",
+    "qx": "qx.npy",
+    "qy": "qy.npy",
+    "bathymetry": "bathymetry.npy",
+    "manning": "manning.npy",
+    "infiltration_zones": "infiltration_zones.npy",
+    "NAP": "NAP.npy"
+  }
+}
\ No newline at end of file
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/quickrun.py b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/quickrun.py
new file mode 100644
index 0000000000000000000000000000000000000000..2b6352abd0f623e3f8422b1a08d0784afb17ff02
--- /dev/null
+++ b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/quickrun.py
@@ -0,0 +1,11 @@
+from pathlib import Path
+from wolfgpu.simple_simulation import SimpleSimulation
+from wolfgpu.SimulationRunner import SimulationRunner
+from wolfgpu.glsimulation import ResultsStore
+
+def main():
+	sim = SimpleSimulation.load(Path(__file__).parent)
+	result_store:ResultsStore = SimulationRunner.quick_run(sim, Path(__file__).parent)
+
+if __name__ == '__main__':
+	main()
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/qx.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/qx.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/qx.npy differ
diff --git a/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/qy.npy b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/qy.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/Q10-5-5/simulations/sim_T25/qy.npy differ
diff --git a/Scenario_wolfgpu/__pycache__/bc_aval_pepinster.cpython-310.pyc b/Scenario_wolfgpu/__pycache__/bc_aval_pepinster.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..e523bc2c039cac9b6428f937bf1971aae9e67234
Binary files /dev/null and b/Scenario_wolfgpu/__pycache__/bc_aval_pepinster.cpython-310.pyc differ
diff --git a/Scenario_wolfgpu/__pycache__/infiltration.cpython-310.pyc b/Scenario_wolfgpu/__pycache__/infiltration.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..1d8a2925ae48b7f0cb3dbf7b8d7974af345bd9c3
Binary files /dev/null and b/Scenario_wolfgpu/__pycache__/infiltration.cpython-310.pyc differ
diff --git a/Scenario_wolfgpu/__pycache__/other_script.cpython-310.pyc b/Scenario_wolfgpu/__pycache__/other_script.cpython-310.pyc
new file mode 100644
index 0000000000000000000000000000000000000000..dba1a3efc82631465550c063a550af238ff375c9
Binary files /dev/null and b/Scenario_wolfgpu/__pycache__/other_script.cpython-310.pyc differ
diff --git a/Scenario_wolfgpu/bathymetry.tif b/Scenario_wolfgpu/bathymetry.tif
new file mode 100644
index 0000000000000000000000000000000000000000..0e092a073567fd2bc11b3cece0c3d0465d04693f
Binary files /dev/null and b/Scenario_wolfgpu/bathymetry.tif differ
diff --git a/Scenario_wolfgpu/bc_aval_pepinster.py b/Scenario_wolfgpu/bc_aval_pepinster.py
new file mode 100644
index 0000000000000000000000000000000000000000..a70d8208e2770cb89edf3ea8fe85ea180f68e7c0
--- /dev/null
+++ b/Scenario_wolfgpu/bc_aval_pepinster.py
@@ -0,0 +1,25 @@
+from wolfhece.scenario.imposebc_void import Impose_Boundary_Conditions
+from wolfhece.wolf_array import WolfArray
+from wolfgpu.simple_simulation import SimpleSimulation, BoundaryConditionsTypes, Direction
+from wolfhece.PyVertexvectors import Zones, zone, vector, wolfvertex
+
+class Impose_BC_Scenario(Impose_Boundary_Conditions):
+
+    def impose_bc(self, simul: SimpleSimulation):
+        """
+        FR
+        --
+        Cette fonction doit être particularisée afin d'appliquer des conditions aux limites à la simulation.
+
+        La simulation est passée en argument et ne doit pas être créee/modifiée.
+
+        EN
+        --
+        This function must be specialized in order to apply boundary conditions to the simulation.
+
+        The simulation is passed as an argument and must not be created/modified.
+        """
+        # Imposition des conditions aux limites - 1-based
+        for j in range(997,1003):
+            simul.add_boundary_condition(i=3, j=j, bc_type=BoundaryConditionsTypes.FROUDE_NORMAL, bc_value=.3, border=Direction.LEFT)
+        pass
diff --git a/Scenario_wolfgpu/deck.tif b/Scenario_wolfgpu/deck.tif
new file mode 100644
index 0000000000000000000000000000000000000000..fa935a78b6fe3debac7f88551c6c70815e9ce90c
Binary files /dev/null and b/Scenario_wolfgpu/deck.tif differ
diff --git a/Scenario_wolfgpu/discharges/Mono25.txt b/Scenario_wolfgpu/discharges/Mono25.txt
new file mode 100644
index 0000000000000000000000000000000000000000..54a1a7022e8c8cfe4ca5d33a4aa7543d8415bec5
--- /dev/null
+++ b/Scenario_wolfgpu/discharges/Mono25.txt
@@ -0,0 +1,7 @@
+Temps	Vesdre	Hoëgne	Wayai
+0	10	5	5
+3600	10	5	5
+7200	20	10	5
+43200	100	50	5
+86400	10	5	5
+100000	10	5	5
\ No newline at end of file
diff --git a/Scenario_wolfgpu/discharges/T10.txt b/Scenario_wolfgpu/discharges/T10.txt
new file mode 100644
index 0000000000000000000000000000000000000000..20f3b9efef0d6188228a4d7e55ffef9d3f5dfbc1
--- /dev/null
+++ b/Scenario_wolfgpu/discharges/T10.txt
@@ -0,0 +1,3 @@
+Temps	Vesdre	Hoëgne	Wayai
+0	10	5	5
+999999	10	5	5
\ No newline at end of file
diff --git a/Scenario_wolfgpu/discharges/T100.txt b/Scenario_wolfgpu/discharges/T100.txt
new file mode 100644
index 0000000000000000000000000000000000000000..20f3b9efef0d6188228a4d7e55ffef9d3f5dfbc1
--- /dev/null
+++ b/Scenario_wolfgpu/discharges/T100.txt
@@ -0,0 +1,3 @@
+Temps	Vesdre	Hoëgne	Wayai
+0	10	5	5
+999999	10	5	5
\ No newline at end of file
diff --git a/Scenario_wolfgpu/discharges/T1000.txt b/Scenario_wolfgpu/discharges/T1000.txt
new file mode 100644
index 0000000000000000000000000000000000000000..751127980f3eb54cda2318c6223882affa18ffe9
--- /dev/null
+++ b/Scenario_wolfgpu/discharges/T1000.txt
@@ -0,0 +1,3 @@
+Temps	Vesdre	Hoëgne	Wayai
+0	10	5	5
+999999	10	5	5
diff --git a/Scenario_wolfgpu/discharges/T25.txt b/Scenario_wolfgpu/discharges/T25.txt
new file mode 100644
index 0000000000000000000000000000000000000000..1f55dad70595a1387f0b9d60f32c7a787b32999e
--- /dev/null
+++ b/Scenario_wolfgpu/discharges/T25.txt
@@ -0,0 +1,2 @@
+Temps	Vesdre	Hoëgne	Wayai
+0	10	5	5
\ No newline at end of file
diff --git a/Scenario_wolfgpu/discharges/T5.txt b/Scenario_wolfgpu/discharges/T5.txt
new file mode 100644
index 0000000000000000000000000000000000000000..20f3b9efef0d6188228a4d7e55ffef9d3f5dfbc1
--- /dev/null
+++ b/Scenario_wolfgpu/discharges/T5.txt
@@ -0,0 +1,3 @@
+Temps	Vesdre	Hoëgne	Wayai
+0	10	5	5
+999999	10	5	5
\ No newline at end of file
diff --git a/Scenario_wolfgpu/discharges/T50.txt b/Scenario_wolfgpu/discharges/T50.txt
new file mode 100644
index 0000000000000000000000000000000000000000..1f55dad70595a1387f0b9d60f32c7a787b32999e
--- /dev/null
+++ b/Scenario_wolfgpu/discharges/T50.txt
@@ -0,0 +1,2 @@
+Temps	Vesdre	Hoëgne	Wayai
+0	10	5	5
\ No newline at end of file
diff --git a/Scenario_wolfgpu/infiltration.py b/Scenario_wolfgpu/infiltration.py
new file mode 100644
index 0000000000000000000000000000000000000000..4f0e45cb72de90b5a4e162a7694f224239d131eb
--- /dev/null
+++ b/Scenario_wolfgpu/infiltration.py
@@ -0,0 +1,72 @@
+from wolfhece.scenario.update_void import Update_Sim
+from wolfhece.wolf_array import WolfArray
+from wolfhece.PyVertexvectors import Zones, zone, vector, wolfvertex
+
+class Update_Sim_Scenario(Update_Sim):
+
+    def update_infiltration(self, infiltration_zones: WolfArray):
+      infiltration_zones.array[477,1087:1092] = 1
+      infiltration_zones.array[477,121] = 2
+      infiltration_zones.array[384,2] = 3
+
+    def update_manning(self, manning: WolfArray):
+        """
+        FR
+        --
+        Cette fonction doit être particularisée afin d'appliquer des modifications
+        à la distribution du coefficient de Manning.
+
+        L'information de base se trouve dans le paramètre 'manning' de classe 'WolfArray'.
+
+        La matrice Numpy, en <Float32>, (avec masque - cf np.ma ) est accessible via :
+          - manning.array[:,:]
+          - manning.array.mask[:,:]
+          - manning.array.data[:,:]
+
+        Il n'est pas permis de remplacer la matrice (création d'une nouvelle matrice et/ou pointage d'una autre matrice). Toutes les operations doivent se faire dans l'esace alloué.
+
+        EN
+        --
+        This function must be customized to apply modifications to the Manning roughness parameter.
+
+        The basic information is located in the 'manning' parameter ('WolfArray' class).
+
+        The Numpy array, in <Float32>, (with a mask - see np.ma), can be accessed via:
+          - manning.array[:,:]
+          - manning.array.mask[:,:]
+          - manning.array.data[:,:]
+
+        It is not allowed to replace the array (creating a new array and/or pointing an other array). All operations must be performed within the allocated space.
+        """
+        pass
+
+    def update_topobathy(self, topobahty: WolfArray):
+        """
+        FR
+        --
+        Cette fonction doit être particularisée afin d'appliquer des modifications
+        à la topographie du lit majeur ou à la bathymétrie du lit mineur.
+
+        L'information de base se trouve dans le paramètre  'topobahty' de classe 'WolfArray'.
+
+        La matrice Numpy, en <Float32>, (avec masque - cf np.ma ) est accessible via :
+          - topobahty.array[:,:]
+          - topobahty.array.mask[:,:]
+          - topobahty.array.data[:,:]
+
+        Il n'est pas permis de remplacer la matrice (création d'une nouvelle matrice et/ou pointage d'una autre matrice). Toutes les operations doivent se faire dans l'esace alloué.
+
+        EN
+        --
+        This function must be customized to apply modifications to the topography of the main bed or the bathymetry of the minor bed.
+
+        The basic information is located in the 'topobahty' parameter ('WolfArray' class).
+
+        The Numpy array, in <Float32>, (with a mask - see np.ma), can be accessed via:
+          - topobahty.array[:,:]
+          - topobahty.array.mask[:,:]
+          - topobahty.array.data[:,:]
+
+        It is not allowed to replace the array (creating a new array and/or pointing an other array). All operations must be performed within the allocated space.
+        """
+        pass
diff --git a/Scenario_wolfgpu/manning.tif b/Scenario_wolfgpu/manning.tif
new file mode 100644
index 0000000000000000000000000000000000000000..fbaa356ef151c923476659485c08e58f2a4eb490
Binary files /dev/null and b/Scenario_wolfgpu/manning.tif differ
diff --git a/Scenario_wolfgpu/other_script.py b/Scenario_wolfgpu/other_script.py
new file mode 100644
index 0000000000000000000000000000000000000000..73f60835b0f5bd3907de9a62821887df4d0784aa
--- /dev/null
+++ b/Scenario_wolfgpu/other_script.py
@@ -0,0 +1,3 @@
+import numpy as np
+
+print('hello world!')
\ No newline at end of file
diff --git a/Scenario_wolfgpu/readme.txt b/Scenario_wolfgpu/readme.txt
new file mode 100644
index 0000000000000000000000000000000000000000..3688223b741bd56af855c7baf809e515124a835d
--- /dev/null
+++ b/Scenario_wolfgpu/readme.txt
@@ -0,0 +1 @@
+Ceci est un exemple de scénario pour le code wolfgpu
\ No newline at end of file
diff --git a/Scenario_wolfgpu/roof.tif b/Scenario_wolfgpu/roof.tif
new file mode 100644
index 0000000000000000000000000000000000000000..fa935a78b6fe3debac7f88551c6c70815e9ce90c
Binary files /dev/null and b/Scenario_wolfgpu/roof.tif differ
diff --git a/Scenario_wolfgpu/var/__bath_assembly.vrt b/Scenario_wolfgpu/var/__bath_assembly.vrt
new file mode 100644
index 0000000000000000000000000000000000000000..08279a1c3c9190fb659aacb5f48c35082545ce83
--- /dev/null
+++ b/Scenario_wolfgpu/var/__bath_assembly.vrt
@@ -0,0 +1,16 @@
+<VRTDataset rasterXSize="480" rasterYSize="1160">
+  <SRS dataAxisToSRSAxisMapping="1,2">PROJCS["BD72 / Belgian Lambert 72",GEOGCS["BD72",DATUM["Reseau_National_Belge_1972",SPHEROID["International 1924",6378388,297,AUTHORITY["EPSG","7022"]],AUTHORITY["EPSG","6313"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4313"]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["latitude_of_origin",90],PARAMETER["central_meridian",4.36748666666667],PARAMETER["standard_parallel_1",51.1666672333333],PARAMETER["standard_parallel_2",49.8333339],PARAMETER["false_easting",150000.013],PARAMETER["false_northing",5400088.438],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","31370"]]</SRS>
+  <GeoTransform>  2.5100000000000000e+05,  5.0000000000000000e+00,  0.0000000000000000e+00,  1.4130000000000000e+05,  0.0000000000000000e+00, -5.0000000000000000e+00</GeoTransform>
+  <VRTRasterBand dataType="Float32" band="1">
+    <NoDataValue>99999</NoDataValue>
+    <ColorInterp>Gray</ColorInterp>
+    <ComplexSource resampling="bilinear">
+      <SourceFilename relativeToVRT="1">..\bathymetry.tif</SourceFilename>
+      <SourceBand>1</SourceBand>
+      <SourceProperties RasterXSize="480" RasterYSize="1160" DataType="Float32" BlockXSize="480" BlockYSize="4" />
+      <SrcRect xOff="0" yOff="0" xSize="480" ySize="1160" />
+      <DstRect xOff="0" yOff="0" xSize="480" ySize="1160" />
+      <NODATA>99999</NODATA>
+    </ComplexSource>
+  </VRTRasterBand>
+</VRTDataset>
diff --git a/Scenario_wolfgpu/var/__bathymetry.tif b/Scenario_wolfgpu/var/__bathymetry.tif
new file mode 100644
index 0000000000000000000000000000000000000000..80698e8f53a1ed1edde5d7821cecbbddb9113569
Binary files /dev/null and b/Scenario_wolfgpu/var/__bathymetry.tif differ
diff --git a/Scenario_wolfgpu/var/__bathymetry_after_scripts.tif b/Scenario_wolfgpu/var/__bathymetry_after_scripts.tif
new file mode 100644
index 0000000000000000000000000000000000000000..71a1c41ca9b2147ed7d9d36e5f64e62205880b7f
Binary files /dev/null and b/Scenario_wolfgpu/var/__bathymetry_after_scripts.tif differ
diff --git a/Scenario_wolfgpu/var/__infiltration_after_scripts.tif b/Scenario_wolfgpu/var/__infiltration_after_scripts.tif
new file mode 100644
index 0000000000000000000000000000000000000000..af886f11e3cd339b8d47bb61d11d83c601cd27b2
Binary files /dev/null and b/Scenario_wolfgpu/var/__infiltration_after_scripts.tif differ
diff --git a/Scenario_wolfgpu/var/__mann_assembly.vrt b/Scenario_wolfgpu/var/__mann_assembly.vrt
new file mode 100644
index 0000000000000000000000000000000000000000..74a8e50ea9faa86ecd6b61143df48cf54de7c165
--- /dev/null
+++ b/Scenario_wolfgpu/var/__mann_assembly.vrt
@@ -0,0 +1,16 @@
+<VRTDataset rasterXSize="480" rasterYSize="1160">
+  <SRS dataAxisToSRSAxisMapping="1,2">PROJCS["BD72 / Belgian Lambert 72",GEOGCS["BD72",DATUM["Reseau_National_Belge_1972",SPHEROID["International 1924",6378388,297,AUTHORITY["EPSG","7022"]],AUTHORITY["EPSG","6313"]],PRIMEM["Greenwich",0,AUTHORITY["EPSG","8901"]],UNIT["degree",0.0174532925199433,AUTHORITY["EPSG","9122"]],AUTHORITY["EPSG","4313"]],PROJECTION["Lambert_Conformal_Conic_2SP"],PARAMETER["latitude_of_origin",90],PARAMETER["central_meridian",4.36748666666667],PARAMETER["standard_parallel_1",51.1666672333333],PARAMETER["standard_parallel_2",49.8333339],PARAMETER["false_easting",150000.013],PARAMETER["false_northing",5400088.438],UNIT["metre",1,AUTHORITY["EPSG","9001"]],AXIS["Easting",EAST],AXIS["Northing",NORTH],AUTHORITY["EPSG","31370"]]</SRS>
+  <GeoTransform>  2.5100000000000000e+05,  5.0000000000000000e+00,  0.0000000000000000e+00,  1.4130000000000000e+05,  0.0000000000000000e+00, -5.0000000000000000e+00</GeoTransform>
+  <VRTRasterBand dataType="Float32" band="1">
+    <NoDataValue>99999</NoDataValue>
+    <ColorInterp>Gray</ColorInterp>
+    <ComplexSource resampling="bilinear">
+      <SourceFilename relativeToVRT="1">..\manning.tif</SourceFilename>
+      <SourceBand>1</SourceBand>
+      <SourceProperties RasterXSize="480" RasterYSize="1160" DataType="Float32" BlockXSize="480" BlockYSize="4" />
+      <SrcRect xOff="0" yOff="0" xSize="480" ySize="1160" />
+      <DstRect xOff="0" yOff="0" xSize="480" ySize="1160" />
+      <NODATA>99999</NODATA>
+    </ComplexSource>
+  </VRTRasterBand>
+</VRTDataset>
diff --git a/Scenario_wolfgpu/var/__manning.tif b/Scenario_wolfgpu/var/__manning.tif
new file mode 100644
index 0000000000000000000000000000000000000000..13a22b2c8dfbaf39f6081a2e098d080f6a0646f4
Binary files /dev/null and b/Scenario_wolfgpu/var/__manning.tif differ
diff --git a/Scenario_wolfgpu/var/__manning_after_scripts.tif b/Scenario_wolfgpu/var/__manning_after_scripts.tif
new file mode 100644
index 0000000000000000000000000000000000000000..b72534d15572482a1441fc969b3b7fca70db0882
Binary files /dev/null and b/Scenario_wolfgpu/var/__manning_after_scripts.tif differ
diff --git a/Scenario_wolfgpu/var/quickruns.bat b/Scenario_wolfgpu/var/quickruns.bat
new file mode 100644
index 0000000000000000000000000000000000000000..194e39ba8a159ab56bc51400eaf0f52fdd3a9a34
--- /dev/null
+++ b/Scenario_wolfgpu/var/quickruns.bat
@@ -0,0 +1,5 @@
+D:
+cd D:\ProgrammationGitLab\HECEPython\tests\data\GPU_Scen\var\simulations\sim_Mono25
+python quickrun.py
+cd D:\ProgrammationGitLab\HECEPython\tests\data\GPU_Scen\var\simulations\sim_T25
+python quickrun.py
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/NAP.npy b/Scenario_wolfgpu/var/simulations/sim_Mono25/NAP.npy
new file mode 100644
index 0000000000000000000000000000000000000000..4d7f4d8355c6692d860b9aa69abdc03adac17531
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_Mono25/NAP.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/bathymetry.npy b/Scenario_wolfgpu/var/simulations/sim_Mono25/bathymetry.npy
new file mode 100644
index 0000000000000000000000000000000000000000..125bb3bdc40899c3fe6fef1965ad724a6b89f712
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_Mono25/bathymetry.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/h.npy b/Scenario_wolfgpu/var/simulations/sim_Mono25/h.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_Mono25/h.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/infiltration_zones.npy b/Scenario_wolfgpu/var/simulations/sim_Mono25/infiltration_zones.npy
new file mode 100644
index 0000000000000000000000000000000000000000..9e40e78cb4a2fcb263ce31074a12d9fcb99dd446
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_Mono25/infiltration_zones.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/manning.npy b/Scenario_wolfgpu/var/simulations/sim_Mono25/manning.npy
new file mode 100644
index 0000000000000000000000000000000000000000..ede955942e91c3c995ff5fbbca1f42b9c3a544ec
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_Mono25/manning.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/parameters.json b/Scenario_wolfgpu/var/simulations/sim_Mono25/parameters.json
new file mode 100644
index 0000000000000000000000000000000000000000..ef1a4c3d53273be7d8bf61a0a51aff154fbd0d38
--- /dev/null
+++ b/Scenario_wolfgpu/var/simulations/sim_Mono25/parameters.json
@@ -0,0 +1,130 @@
+{
+  "spec_version": "0.1",
+  "comment": null,
+  "creation_date": "21/12/2023 14:48:23",
+  "parameters": {
+    "courant": 0.4,
+    "froude_max": 3,
+    "nx": 480,
+    "ny": 1160,
+    "dx": 5.0,
+    "dy": 5.0,
+    "timestep_strategy": "OPTIMIZED_TIME_STEP",
+    "timestep": 0.0,
+    "infiltration_interp": "NONE",
+    "duration": {
+      "type": "SECONDS",
+      "duration": 100000.0
+    },
+    "report_period": {
+      "type": "SECONDS",
+      "duration": 3600.0
+    },
+    "runge_kutta": 0.5,
+    "base_coord_x": 251000.0,
+    "base_coord_y": 135500.0
+  },
+  "boundary_conditions": [
+    {
+      "i": 3,
+      "j": 997,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 998,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 999,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1000,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1001,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1002,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    }
+  ],
+  "infiltrations": [
+    {
+      "start_time": 0.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 3600.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 7200.0,
+      "q": [
+        20.0,
+        10.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 43200.0,
+      "q": [
+        100.0,
+        50.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 86400.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 100000.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    }
+  ],
+  "maps": {
+    "h": "h.npy",
+    "qx": "qx.npy",
+    "qy": "qy.npy",
+    "bathymetry": "bathymetry.npy",
+    "manning": "manning.npy",
+    "infiltration_zones": "infiltration_zones.npy",
+    "NAP": "NAP.npy"
+  }
+}
\ No newline at end of file
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/quickrun.py b/Scenario_wolfgpu/var/simulations/sim_Mono25/quickrun.py
new file mode 100644
index 0000000000000000000000000000000000000000..2b6352abd0f623e3f8422b1a08d0784afb17ff02
--- /dev/null
+++ b/Scenario_wolfgpu/var/simulations/sim_Mono25/quickrun.py
@@ -0,0 +1,11 @@
+from pathlib import Path
+from wolfgpu.simple_simulation import SimpleSimulation
+from wolfgpu.SimulationRunner import SimulationRunner
+from wolfgpu.glsimulation import ResultsStore
+
+def main():
+	sim = SimpleSimulation.load(Path(__file__).parent)
+	result_store:ResultsStore = SimulationRunner.quick_run(sim, Path(__file__).parent)
+
+if __name__ == '__main__':
+	main()
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/qx.npy b/Scenario_wolfgpu/var/simulations/sim_Mono25/qx.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_Mono25/qx.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_Mono25/qy.npy b/Scenario_wolfgpu/var/simulations/sim_Mono25/qy.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_Mono25/qy.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/NAP.npy b/Scenario_wolfgpu/var/simulations/sim_T25/NAP.npy
new file mode 100644
index 0000000000000000000000000000000000000000..4d7f4d8355c6692d860b9aa69abdc03adac17531
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_T25/NAP.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/bathymetry.npy b/Scenario_wolfgpu/var/simulations/sim_T25/bathymetry.npy
new file mode 100644
index 0000000000000000000000000000000000000000..125bb3bdc40899c3fe6fef1965ad724a6b89f712
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_T25/bathymetry.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/h.npy b/Scenario_wolfgpu/var/simulations/sim_T25/h.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_T25/h.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/infiltration_zones.npy b/Scenario_wolfgpu/var/simulations/sim_T25/infiltration_zones.npy
new file mode 100644
index 0000000000000000000000000000000000000000..9e40e78cb4a2fcb263ce31074a12d9fcb99dd446
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_T25/infiltration_zones.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/manning.npy b/Scenario_wolfgpu/var/simulations/sim_T25/manning.npy
new file mode 100644
index 0000000000000000000000000000000000000000..ede955942e91c3c995ff5fbbca1f42b9c3a544ec
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_T25/manning.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/parameters.json b/Scenario_wolfgpu/var/simulations/sim_T25/parameters.json
new file mode 100644
index 0000000000000000000000000000000000000000..2986c4133a8030fa5705661181186ef3ecc268b3
--- /dev/null
+++ b/Scenario_wolfgpu/var/simulations/sim_T25/parameters.json
@@ -0,0 +1,98 @@
+{
+  "spec_version": "0.1",
+  "comment": null,
+  "creation_date": "21/12/2023 14:48:29",
+  "parameters": {
+    "courant": 0.4,
+    "froude_max": 3,
+    "nx": 480,
+    "ny": 1160,
+    "dx": 5.0,
+    "dy": 5.0,
+    "timestep_strategy": "OPTIMIZED_TIME_STEP",
+    "timestep": 0.0,
+    "infiltration_interp": "NONE",
+    "duration": {
+      "type": "SECONDS",
+      "duration": 86400.0
+    },
+    "report_period": {
+      "type": "SECONDS",
+      "duration": 3600.0
+    },
+    "runge_kutta": 0.5,
+    "base_coord_x": 251000.0,
+    "base_coord_y": 135500.0
+  },
+  "boundary_conditions": [
+    {
+      "i": 3,
+      "j": 997,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 998,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 999,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1000,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1001,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    },
+    {
+      "i": 3,
+      "j": 1002,
+      "type": "FROUDE_NORMAL",
+      "val": 0.3,
+      "direction": "LEFT"
+    }
+  ],
+  "infiltrations": [
+    {
+      "start_time": 0.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    },
+    {
+      "start_time": 86400.0,
+      "q": [
+        10.0,
+        5.0,
+        5.0
+      ]
+    }
+  ],
+  "maps": {
+    "h": "h.npy",
+    "qx": "qx.npy",
+    "qy": "qy.npy",
+    "bathymetry": "bathymetry.npy",
+    "manning": "manning.npy",
+    "infiltration_zones": "infiltration_zones.npy",
+    "NAP": "NAP.npy"
+  }
+}
\ No newline at end of file
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/quickrun.py b/Scenario_wolfgpu/var/simulations/sim_T25/quickrun.py
new file mode 100644
index 0000000000000000000000000000000000000000..2b6352abd0f623e3f8422b1a08d0784afb17ff02
--- /dev/null
+++ b/Scenario_wolfgpu/var/simulations/sim_T25/quickrun.py
@@ -0,0 +1,11 @@
+from pathlib import Path
+from wolfgpu.simple_simulation import SimpleSimulation
+from wolfgpu.SimulationRunner import SimulationRunner
+from wolfgpu.glsimulation import ResultsStore
+
+def main():
+	sim = SimpleSimulation.load(Path(__file__).parent)
+	result_store:ResultsStore = SimulationRunner.quick_run(sim, Path(__file__).parent)
+
+if __name__ == '__main__':
+	main()
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/qx.npy b/Scenario_wolfgpu/var/simulations/sim_T25/qx.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_T25/qx.npy differ
diff --git a/Scenario_wolfgpu/var/simulations/sim_T25/qy.npy b/Scenario_wolfgpu/var/simulations/sim_T25/qy.npy
new file mode 100644
index 0000000000000000000000000000000000000000..91e808caa68be2d2c7e78b0dc806e01bf741cc71
Binary files /dev/null and b/Scenario_wolfgpu/var/simulations/sim_T25/qy.npy differ