PoroStrip_v1.0 – A code to compute anisotropic porosity fields for urban flood modelling

This repository contains the code "PoroStrip_v1.0". It allows obtaining the spatially-distributed porosity fields for anisotropic porosity models in differential form, that can be used for simulating large-scale urban floods. The method and its implementation are described in the paper: Ferrari, A., and Viero, D.P., 2020. Floodwater pathways in urban areas: a method to compute porosity fields for anisotropic subgrid models in differential form, Journal of Hydrology 589, 125193, https://doi.org/10.1016/j.jhydrol.2020.125193. The repository also contains the input (and some output) files for the three case studies described in the above paper. This repository contains the following folders: ⇒"Code_Porostrip_v1.0" with the code "PoroStrip_v1.0" available as both compiled and source code. The folder contains: •"Executable" with: -the Win86 executable file "PoroStrip_v1.0.exe"; -the "PoroStrip_v1.0_cfg.txt" configuration file, which determines the output files to be created; •"Source" with the source code (ASCII files), where some routines and intrinsic routine calls are specific for Intel Fortan; •the Intel Fortran project files. ⇒"CaseStudy1_Toce", "CaseStudy2_Spinea", and "CaseStudy3_Palmanova" folders, referring to the case study described in the companion paper. Each case study folder contains: •"input" with the files that have to be provided to the code as inputs: -the .bln file with the coordinates of the obstacle footprints; -the .asc files with the characteristics of the chosen (coarse) grid, structured as the header of ESRI ASCII grid files (xLLcenter, yLLcenter, cellsize, ncols, nrows, NODATA_VALUE); •"ouput" with some of the output files produced by the code for each tested grid size and strip width; •"visualization": the .dxf files produced by the code have been modified with AutoCAD Map3D to produce a .pdf file helpfull for the overall visualization of the results (e.g. location of the obstacles and resulting roseplots of conveyance porosity). The alternative .dxf version generated with AutoCAD is denoted by the "_prnt" suffix ("_prnt.dxf"). The .pdf files have been structured with the following switchable layers: -CELLS showing the (coarse) grid cells (pink-lined squares); -OBSTACLES representing the building footprint (gray areas); -PSI_L showing the L principal direction of maximum conveyance (blue segments); -PSI_T showing the T principal direction othogonal to L (red segments); -ROSEPLOTS representing the Psi roseplots (green lines). ⇒the "Layers_PorosityMaps.las" that can be used to assign proper styles to the AutoDesk AutoCAD layers of the dxf files produced by the "Porostrip_v1.0.exe" code. See also the "README.txt" file contained in this repository.

本仓库包含代码"PoroStrip_v1.0"，该工具可获取微分形式各向异性孔隙率模型的空间分布孔隙率场，可用于模拟大规模城市洪水。该方法及其实现详见下述论文：Ferrari, A. 与 Viero, D.P. 于2020年发表于《Journal of Hydrology》第589卷的《Floodwater pathways in urban areas: a method to compute porosity fields for anisotropic subgrid models in differential form》，DOI链接为https://doi.org/10.1016/j.jhydrol.2020.125193。

本仓库还包含上述论文中提及的三个案例研究的输入（及部分输出）文件。本仓库包含以下文件夹：
1. "Code_Porostrip_v1.0"：存放可执行版与源码版的"PoroStrip_v1.0"代码，该文件夹内含：
- "Executable"子文件夹：包含Win86架构可执行文件"PoroStrip_v1.0.exe"，以及用于指定待生成输出文件的配置文件"PoroStrip_v1.0_cfg.txt"；
- "Source"子文件夹：内含源代码（ASCII格式文件），其中部分子程序与内置函数调用专为Intel Fortran设计；
- Intel Fortran项目文件。
2. "CaseStudy1_Toce"、"CaseStudy2_Spinea"与"CaseStudy3_Palmanova"文件夹：对应上述论文描述的三个案例研究。每个案例研究文件夹包含：
- "input"子文件夹：存放代码所需的输入文件：
- .bln格式文件：存储障碍物轮廓坐标；
- .asc格式文件：包含所选（粗分辨率）网格的特征信息，格式遵循ESRI ASCII网格文件（ESRI ASCII grid files）头规范，包含xLLcenter、yLLcenter、cellsize、ncols、nrows、NODATA_VALUE字段；
- "output"子文件夹：存放代码针对每种测试网格尺寸与条带宽度生成的部分输出文件；
- "visualization"子文件夹：代码生成的.dxf文件已通过AutoCAD Map3D进行修改，以生成便于整体可视化结果的.pdf文件（例如障碍物位置与输导孔隙率的玫瑰图）。通过AutoCAD生成的替代.dxf版本带有"_prnt"后缀（即"_prnt.dxf"）。该.pdf文件包含以下可切换图层：
- CELLS：显示（粗分辨率）网格单元（粉色边框正方形）；
- OBSTACLES：表示建筑轮廓（灰色区域）；
- PSI_L：显示最大输导率的L主方向（蓝色线段）；
- PSI_T：显示与L正交的T主方向（红色线段）；
- ROSEPLOTS：表示Psi玫瑰图（绿色线条）。
3. "Layers_PorosityMaps.las"文件：可用于为"Porostrip_v1.0.exe"代码生成的DXF文件配置欧特克AutoCAD（AutoDesk AutoCAD）图层的合适样式。

请参阅本仓库内含的"README.txt"文件。