from IPython.core.display import HTML
css_file = '../style/style.css'
HTML(open(css_file, "r").read())
Define the function topo_withDownstreamWall that receives two arrays $x$ and $y$ and returns another array with the same topography as in the example of Part 1 in this workshop, but with the addition of a wall in the downstream channel defined by a 15m height in the region where $(x,y)\in[50,55]\times[-75,25]$.
Define a function that writes the topography function in 1. into a file 'channel_withDownstreamWall.topotype2' with $topo\_type=2$ and same discretization as in Part 1.
Similarly define a function qinit_withDownstreamWall that represents the initial condition, as in the Part 1, and write it to a file called 'qinit_withDownstreamWall.topotype2' with $topo\_type=2$.
From a new cell in this notebook write a setrun.py file that can read the new data. It must satisfy the following conditions:
Run the simulation and show a video and punctual results in the same points as in Part 1. Cmpare with Part 1. Comment.