From 4ba4a2fbd8c0ddc9d55dcb069650457f838386ba Mon Sep 17 00:00:00 2001 From: RuudHurkmans <59835874+RuudHurkmans@users.noreply.github.com> Date: Fri, 29 Jan 2021 14:42:37 +0100 Subject: [PATCH] Bug fix: land use type bares soil was (incorrectly) not taken into account. --- delft3dfmpy/__init__.py | 2 +- delft3dfmpy/converters/hydamo_to_dflowrr.py | 2 +- 2 files changed, 2 insertions(+), 2 deletions(-) diff --git a/delft3dfmpy/__init__.py b/delft3dfmpy/__init__.py index 3ab9285..7c851d1 100644 --- a/delft3dfmpy/__init__.py +++ b/delft3dfmpy/__init__.py @@ -4,7 +4,7 @@ __author__ = """Guus Rongen""" __email__ = 'rongen@hkv.nl' -__version__ = '1.1.0' +__version__ = '1.1.2' from delft3dfmpy.core.dfm import DFlowFMModel from delft3dfmpy.core.drr import DFlowRRModel diff --git a/delft3dfmpy/converters/hydamo_to_dflowrr.py b/delft3dfmpy/converters/hydamo_to_dflowrr.py index e9c251e..ba701e9 100644 --- a/delft3dfmpy/converters/hydamo_to_dflowrr.py +++ b/delft3dfmpy/converters/hydamo_to_dflowrr.py @@ -86,7 +86,7 @@ def generate_unpaved(catchments, landuse, surface_level, soiltype, surface_stor #ms = [ms for ms in meteo_areas.itertuples() if ms.geometry.contains(cat.geometry.centroid)] #ms = ms[0] if ms != [] else meteo_areas.iloc[0,:][0] mapping = np.zeros(16, dtype=int) - for i in range(1,12): + for i in range(1,13): if i in lu_counts[num]: mapping[sobek_indices[i-1]-1] = lu_counts[num][i]*px_area lu_map = ' '.join(map(str,mapping)) elev = mean_elev[num]['median']