In [1]:
%matplotlib inline
import python_subgrid.wrapper
import python_subgrid.plotting
import matplotlib.cm
import matplotlib.colors
import matplotlib.pyplot as plt
import numpy as np
import logging
import shapely.ops
import scipy.spatial
import scipy.interpolate
import shapely.geometry 
import python_subgrid.plotting
import json

Initialisation


In [2]:
subgrid = python_subgrid.wrapper.SubgridWrapper(mdu="/Users/baart_f/models/kockengen/kockengen_kockengen.mdu")
#subgrid = python_subgrid.wrapper.SubgridWrapper(mdu="/home/fedor/Checkouts/models/delfland-model-voor-3di/hhdlipad.mdu")

python_subgrid.wrapper.logger.setLevel(logging.WARN)

In [3]:
# lot's of problems
subgrid.start()


WARNING:python_subgrid.wrapper:File admin/gridkockengen_kockengen.grd  does not exist. Grid will be generated.
WARNING:python_subgrid.wrapper:Could not open grid admin file 'admin/gridkockengen_kockengen.grd' for writing; iostat=              2
WARNING:python_subgrid.wrapper:Crest level of structure connection_1 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_1 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_1 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_1 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_1 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_1 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_1 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_2 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_2 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_2 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_2 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_2 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_2 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_2 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_3 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_4 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_4 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_4 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_4 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_4 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_4 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_5 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_5 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_5 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_5 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_5 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_6 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_6 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_6 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_6 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_6 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_6 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_7 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_8 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_8 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_8 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_8 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_8 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_9 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_9 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_9 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_9 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_9 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_9 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_9 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_9 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_10 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_10 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_10 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_10 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_10 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_10 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_10 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_10 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_11 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_12 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_13 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_14 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_14 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_14 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_14 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_14 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_14 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_14 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_15 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_15 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_15 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_15 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_15 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_15 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_16 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_17 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_18 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_19 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_20 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_21 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_22 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_23 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_24 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_25 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_26 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_26 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_26 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_26 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_26 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_26 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_26 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_27 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_28 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_29 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_30 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_31 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_31 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_31 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_31 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_31 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_31 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_31 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_31 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_32 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_32 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_32 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_32 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_33 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_33 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_33 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_33 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_33 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_34 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_34 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_34 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_34 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_34 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_35 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_35 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_35 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_35 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_35 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_36 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_37 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_37 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_38 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_39 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_40 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_41 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_44 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_44 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_44 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_44 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_44 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_44 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_44 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_44 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_45 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_46 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_47 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_47 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_47 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_47 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_47 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_47 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_48 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_48 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_48 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_48 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_48 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_49 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_49 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_50 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_50 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_50 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_50 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_50 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_50 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_51 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_52 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_53 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_53 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_53 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_53 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_53 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_53 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_53 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_53 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_54 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_55 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_55 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_55 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_55 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_55 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_55 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_55 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_57 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_57 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_57 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_57 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_57 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_57 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_58 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_58 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_58 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_58 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_58 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_58 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_61 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_61 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_61 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_61 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_61 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_62 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_62 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_62 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_62 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_62 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_63 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_65 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_66 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_67 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_68 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_69 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_69 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_69 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_69 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_69 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_69 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_69 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_70 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_70 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_70 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_70 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_70 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_70 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_70 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_71 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_72 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_73 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_74 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_75 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_75 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_75 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_76 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_76 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_76 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_76 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_77 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_77 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_77 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_77 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_77 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_77 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_77 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_77 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_78 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_78 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_78 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_78 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_78 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_78 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_78 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_79 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_80 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_81 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_81 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_81 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_81 is below bed level of adjacent cells, crest level is raised to bed level
WARNING:python_subgrid.wrapper:Crest level of structure connection_81 is below bed level of adjacent cells, crest level is raised to bed level

In [4]:
# Get the variables needed or the levee update
grid = {}
for var in {'FlowElem_xcc', 'FlowElem_ycc', 'nod_type', 'FlowElemContour_x', 
            'FlowElemContour_y', 'dmax', 'dps', 'x0p', 'y0p', 'x1p', 'y1p',
            'dxp', 'dyp',
            'levnd0', 'levnd1', 'xleve', 'yleve', 'dlev', 'xlevnd', 'ylevnd',
            'aii', 'bi', 'aijl'}:
    value = subgrid.get_nd(var, sliced=True)
    grid[var] = value
grid['quad_grid'] = python_subgrid.plotting.make_quad_grid(subgrid)


grid['xu'] = subgrid.get_nd('FlowLink_xu', sliced=False)
grid['yu'] = subgrid.get_nd('FlowLink_yu', sliced=False)
grid['line'] = subgrid.get_nd('FlowLink', sliced=False)
grid['link_type'] = subgrid.get_nd('link_type', sliced=False)
# fill the boolean indices
grid['liutot'] = subgrid.get_nd('liutot')
grid['livtot'] = subgrid.get_nd('livtot')


grid['is_1d'] = grid['link_type'] == 1
grid['line'] = grid['line'][:,grid['is_1d']]
grid['xu'] = grid['xu'][grid['is_1d']]
grid['yu'] = grid['yu'][grid['is_1d']]
grid['is_u'] = np.zeros_like(grid['xu'], dtype='bool') 
grid['is_v'] = np.zeros_like(grid['xu'], dtype='bool') 

grid['is_u'][:grid['liutot']] = True
grid['is_v'][(grid['liutot']):(grid['liutot']+grid['livtot'])] = True
assert np.all(grid['is_u'] == ~grid['is_v'])

Convergence check


In [5]:
# convergence problem is not traced yet. 
subgrid.update(-1)
fig, (ax) = plt.subplots()

for name, color in  [('aii', 'k'), ('bi', 'b'), ('aijl', 'g')]:
    var = subgrid.get_nd(name)
    ax.plot(var, color + '.', alpha=0.3, label=name)
ax.legend()
ax.set_yscale('symlog')
subgrid.get_nd('aijl')


WARNING:python_subgrid.wrapper:File admin/tablekockengen_kockengen.tbl  does not exist. Tables will be generated.
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000010
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
Out[5]:
array([       0.        ,    -2632.71901632,        0.        , ...,
       -1444842.29625433, -1444842.29625433,  -722421.14812716])

Checks


In [6]:
for i in range(1):
    subgrid.update(-1)
vars = {}
for var in {'s1', 'q', 'vol1'}:
    vars[var] = subgrid.get_nd(var, sliced=True).copy()
vars['wl'] = vars['s1'] - -grid['dmax']
colors = matplotlib.cm.Blues(matplotlib.colors.Normalize(vmin=vars['wl'].min(), vmax=vars['wl'].max())(vars['wl']))


WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000

In [7]:
import networkx as nx

# in case you want to do topological things, quickly determine if a levee ring is breached
# we need something like this
# for now ignore
if 'xlevnd' in grid: 
    # we want to have separate levees
    # lookup the split points
    start = grid['levnd0'] - 1
    end = grid['levnd1'] - 1
    xlev = grid['xlevnd']
    ylev = grid['ylevnd']

    edges = np.c_[start, end]
    G = nx.Graph()
    G.add_edges_from(edges)
    print nx.number_connected_components(G)
    lines = list(nx.DiGraph(g.edges()) for g in nx.connected_component_subgraphs(G))
    for line in lines:
        for edge in line.edges():
            plt.plot([xlev[edge[0]], xlev[edge[1]]], [ylev[edge[0]],ylev[edge[1]]], 'k-')


90

In [8]:
# we treat levees as a collections of unconnected independent closed line segments
# this is just plotting
if 'xlevnd' in grid: 
    # we want to have separate levees
    # lookup the split points
    start = grid['levnd0'] - 1
    end = grid['levnd1'] - 1
    xlev = grid['xlevnd']
    ylev = grid['ylevnd']
for edge in zip(start, end):
    plt.plot([xlev[edge[0]], xlev[edge[1]]], [ylev[edge[0]],ylev[edge[1]]], 'k-')


Experiment


In [9]:
# we want to test out this example
manhole =     {
        "itype": 1, 
        "timestep_start": 16, 
        "hash": "0bba855b78e8042b9ae9d8ba1c9a1695d3fe34c6", 
        "timestep_end": None, 
        "wgs84_y": "52.138430275373516", 
        "wgs84_x": "4.979805351667136", 
        "amount": 2.0, 
        "y": 461215.50900446286, 
        "x": 127109.67984852854, 
        "type": "manhole", 
        "unique_id": "00f81085-144a-48f8-aaf6-9d0b50521142"
    }


polygon = shapely.geometry.shape({
    "type": "Polygon", 
    "coordinates": [[
                [126966.098,461322.395],
                [126926.560,461235.412],
                [127057.299,461170.570],
                [127095.782,461247.009]
                ]] }
)
# go 5 meters deep
breach = 5

In [10]:
# let's make a multilinestring object that we can use for spatial indexing
if 'xlevnd' in grid: 
    # we want to have separate levees
    
    x = grid['FlowElem_xcc'][grid['nod_type']==1]
    y = grid['FlowElem_ycc'][grid['nod_type']==1]
    xlev = grid['xlevnd']
    ylev = grid['ylevnd']
    points = np.r_[
                 np.c_[x,y], 
                 np.c_[xlev, ylev]
                 ]
    # split up the levees, we now have a list of arrays with start end points
    start = grid['levnd0'] - 1
    end = grid['levnd1'] - 1
    
    # reconstruct the separate lines
    lines = []
    # match up begin and end
    for (seg_begin, seg_end) in zip(start, end):
        lines.append(
            [
                (xlev[seg_begin], ylev[seg_begin]), 
                (xlev[seg_end], ylev[seg_end])
            ]
        )
    levees = shapely.geometry.MultiLineString(lines)
levees.to_wkt()[:100] + '...'


Out[10]:
'MULTILINESTRING ((126845.9613950000057230 458871.7821289999992587, 126857.6895150000054855 458850.35...'

In [11]:
# define a plot function
def render(subgrid, grid, patch, labels=False):
    fig, ax = plt.subplots(figsize=(20,13))
    origin = (grid['x0p'], grid['y0p'])
    extent = (grid['x0p'], grid['x1p'], grid['y0p'], grid['y1p'])
    quad_grid = grid['quad_grid']
    wl = subgrid.get_nd('s1', sliced=True)[quad_grid.filled(-1)]- -grid['dps']
    vol_mask = (subgrid.get_nd('vol1', sliced=True)==0)[quad_grid.filled(-1)]
    wl_masked = np.ma.masked_array(wl, vol_mask)
    vol = subgrid.get_nd('vol1', sliced=True)[quad_grid.filled(-1)]
    im = ax.imshow(wl_masked, cmap='Blues', vmin=0.0, vmax=5, origin=origin, extent=extent)
    #im = ax.imshow(subgrid.get_nd('s1', sliced=True)[quad_grid.filled(-1)], cmap='Blues', vmin=-5, vmax=5, origin=origin, extent=extent)
    xc = grid['FlowElemContour_x']
    yc = grid['FlowElemContour_y']
    verts = [np.c_[xc_i, yc_i] for xc_i, yc_i in zip(xc, yc)]

    poly_collection = matplotlib.collections.PolyCollection(verts, closed=True, linewidth=0.1, transOffset=ax.transData, facecolor='none')
    if 'xlevnd' in grid:
        for i, levee in enumerate(levees):
            plt.plot(levee.xy[0], levee.xy[1], 'k-', linewidth=0.1)
            if labels:
                plt.text(np.mean(levee.xy[0]), np.mean(levee.xy[1]), str(i), fontsize=6)
        for line in levees.intersection(polygon):
            plt.plot(line.xy[0], line.xy[1], 'r-', linewidth=1)
        sc = ax.scatter(grid['xleve'].mean(0), grid['yleve'].mean(0), c=grid['dlev'], edgecolor='none', vmin=-3,vmax=3)
        plt.colorbar(sc, ax=ax)
    plt.colorbar(im, ax=ax)
    
    # split up u and v
    u1 = subgrid.get_nd('u1', sliced=False)[:]
    u1 = u1[grid['is_1d']]
    q = subgrid.get_nd('q', sliced=False)[:]
    q = q[grid['is_1d']]

    u = grid['is_u']*u1
    v = grid['is_v']*u1
    
    q_u = grid['is_u']*q
    q_v = grid['is_v']*q
        
    ax.quiver(grid['xu'], grid['yu'], q_u, q_v, 
               alpha=0.7,
               scale=1000,
               color='red')
    #trans = fig.dpi_scale_trans + matplotlib.transforms.Affine2D().scale(1.0/72.0)
    #poly_collection.set_transform(trans)  # the points to pixels transform
    ax.add_collection(poly_collection)
    ax.add_patch(patch)
    ax.autoscale_view()

In [12]:
# toss in some water
subgrid.discharge(manhole['x'], manhole['y'], manhole['unique_id'], manhole['itype'], manhole['amount']*100)
manhole


Out[12]:
{'amount': 2.0,
 'hash': '0bba855b78e8042b9ae9d8ba1c9a1695d3fe34c6',
 'itype': 1,
 'timestep_end': None,
 'timestep_start': 16,
 'type': 'manhole',
 'unique_id': '00f81085-144a-48f8-aaf6-9d0b50521142',
 'wgs84_x': '4.979805351667136',
 'wgs84_y': '52.138430275373516',
 'x': 127109.67984852854,
 'y': 461215.50900446286}

In [13]:
# let the simulation run for a bit
for i in range(20):
    subgrid.update(-1)


WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000004
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000

In [14]:
# create an area where we want to update the levees
# assume that if  you lower the rectangle, that the area north of it will flood
xy = polygon.exterior.coords
patch = matplotlib.patches.Polygon(xy, facecolor='none')
# plot the model
render(subgrid, grid, patch)
plt.plot(manhole['x'], manhole['y'], 'rx')
# zoom in
plt.xlim(126000, 128000)
plt.ylim(461000, 462000)


Out[14]:
(461000, 462000)

In [14]:


In [15]:
# intersect the levees with the polygon
shapely.geometry.mapping(levees.intersection(polygon))


Out[15]:
{'coordinates': (((126997.679734959, 461200.1390833171),
   (126998.338164, 461202.086183)),
  ((126998.338164, 461202.086183), (127014.072488, 461248.616472)),
  ((127014.072488, 461248.616472), (127027.0665881455, 461286.95363494393))),
 'type': 'MultiLineString'}

In [16]:
import shapely.prepared
# create a prepared polygon for performance
polygon_prepared = shapely.prepared.prep(polygon)

In [17]:
# lookup the levees that are at least partially in the polygon

start_points = np.c_[grid['xleve'][0], grid['yleve'][0]]
end_points = np.c_[grid['xleve'][1], grid['yleve'][1]]

lines = []
for start, end in zip(start_points, end_points):
    line = shapely.geometry.LineString([start, end])
    lines.append(line)
# Use contains if you only want the levees that are fully contained
in_levee = [polygon_prepared.intersects(line) for line in lines]
levee_idx, = np.where(in_levee)
levee_idx


Out[17]:
array([1064, 1065, 1066, 1067, 1068, 1069, 1070, 1071, 1072, 1073, 1074,
       1075, 1076, 1077, 1078])

In [18]:
# lower the bathymetry
geojson = json.dumps(shapely.geometry.mapping(polygon))
grid['dps']
x0 = grid['x0p']
x1 = grid['x1p']
y0 = grid['y0p']
y1 = grid['y1p']
rr, cc = python_subgrid.plotting.draw_shape_on_raster(geojson, grid['dps'], breach, extent=(x0, y0, x1, y1))

quad_cells = set(grid['quad_grid'][row, col] + 1 for row, col in zip(rr, cc))
subgrid.update_tables('dps', list(quad_cells) )
plt.subplots(figsize=(10,8))
grid['dps'] = subgrid.get_nd('dps', sliced=True)
plt.imshow(-np.ma.masked_less(grid['dps'][::3,::3], -100), vmin=-breach, vmax=breach, origin='bottom', cmap='gist_earth')
plt.colorbar()


DEBUG:root:drawing using value 5 and geometry {"type": "Polygon", "coordinates": [[[126966.098, 461322.395], [126926.56, 461235.412], [127057.299, 461170.57], [127095.782, 461247.009], [126966.098, 461322.395]]]} within extent (array(124434.5), array(458773.0), array(129571.5), array(464846.0))
DEBUG:root:drawing in [4796 4796 4796 ..., 5097 5098 5098], [5244 5245 5246 ..., 5066 5063 5064], old values:
 [ 0.95999998  0.97000003  0.95999998 ...,  2.3499999   2.3499999   2.3499999 ]
WARNING:python_subgrid.wrapper:Something broke the friction table, rereading....
Out[18]:
<matplotlib.colorbar.Colorbar instance at 0x163f9acf8>

In [19]:
# lower the levees
grid['dlev'][levee_idx] = -breach
grid['dlev']
# update tables, nodes are ignored for now
subgrid.update_tables('dlev', levee_idx + 1)
grid['dlev'] = subgrid.get_nd('dlev')

In [20]:
# let the simulation run for a bit
for i in range(20):
    subgrid.update(-1)
# plot the model


WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000004
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000006
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000005
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000004
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000005
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000005
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000001
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000008
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000011
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000004
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000004
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000002
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000004
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000005
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000008
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000003
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000005
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000004
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000004
WARNING:python_subgrid.wrapper:!! no convergence for CG simple method. noiter, eps=           1001      0.000000

In [21]:
# create an area where we want to update the levees
# assume that if  you lower the rectangle, that the area south of it will flood
xy = polygon.exterior.coords
patch = matplotlib.patches.Polygon(xy, facecolor='none')
# plot the model
render(subgrid, grid, patch)
plt.plot(manhole['x'], manhole['y'], 'rx')
# zoom in
plt.xlim(126000, 128000)
plt.ylim(461000, 462000)


Out[21]:
(461000, 462000)