Data and code for: From roughness uncertainty to model-configuration effects in 2D flood inundation models: a Monte Carlo comparison of SFINCS and HEC-RAS

This dataset contains the model setups, raw simulation outputs, tabular data and analysis code required to reproduce all results reported in: Navas, S., Galán, Á., & del Jesus, M. (2026). From roughness uncertainty to model-configuration effects in 2D flood inundation models: a Monte Carlo comparison of SFINCS and HEC-RAS. Environmental Modelling & Software (submitted July 2026). One thousand Monte Carlo hydraulic simulations were performed with each model — 1,000 with SFINCS and 1,000 with HEC-RAS — using stochastic Manning roughness coefficient combinations for nine land-use classes in the Besaya River floodplain at Corrales de Buelna, Cantabria, Spain (T = 500-year flood event). Five simulations (0.5%) affected by numerical non-convergence were removed, leaving 995 valid paired simulations for analysis. Version notes (this version): corrects the Manning combination matrix (data/monte_carlo_combinations.csv), which in earlier versions did not match the matrix that actually generated the production simulations, and updates every file derived from it (the HEC-RAS/SFINCS sensitivity tables, the paired comparison table, the sampling notebook, and several figures). This version also brings figures.zip fully in sync with the current manuscript: adds the twelfth main-text figure (the T=500 inflow hydrograph, showing the HEC-RAS 43-hour truncation, added in response to a reviewer request), replaces the superseded combined location figure with the two figures actually used, and adds five supplementary figures that were missing from earlier versions. The raw hydraulic simulation outputs themselves, and their depth and flooded-area metrics, are unchanged. See README.md for the full account. The dataset includes: Tabular inputs: Manning roughness reference values, fitted distributions, 1,000 × 9 Manning combination matrix, hourly reporting series of the T=500 inflow hydrograph (native 2-minute model input retained in the model files).Tabular outputs: per-simulation SFINCS and HEC-RAS statistics, matched 995-pair ensemble, anomalous simulation diagnostics.Raw maximum-depth GeoTIFF rasters for all 1,000 simulations from each model (~290 MB total, 5 m resolution, EPSG:25830).Model setup files for HEC-RAS 6.6 (~1.4 GB) and SFINCS v2 (~7 MB).Seven Jupyter notebooks reproducing the full analysis pipeline.Publication-quality figures (PDF) for all twelve main-text figures and all nine supplementary figures. The tabular data in data/ are sufficient to reproduce all tables and figures without rerunning the hydraulic models. Raw simulation rasters are included for spatial reproducibility. Analysis code uses the pyhydra Python package: https://github.com/navass11/pyhydra. NotesSubmitted to Environmental Modelling & Software (under review, July 2026).