To add the bathymetry we will use the Gebco dataset as used in delft dashboard. We start by using the following line
bathymetryData = GetGebcoBathymetryData(xOffset, yOffset, xOffset + gridWidth, yOffset + gridHeigth, 3857)
Here we retrieve the bathymetry data from the gebco dataset using the extent of our grid and the coordinate system EPSG code. This is needed because the dataset uses a specific (different) coordinate system.
The next step is to remove the cells that are on land. We can determine this by using our bathymetry data. Add and run the following line to remove all land cells that have a depth above -2
CleanupLandCells(fmModel, bathymetryData, -2)
The grid should now look something like this:
The last step is to add the bathymetry to the model (defined as depth at the vertex) with the following lines :
values = []
for vertex in fmModel.Grid.Vertices :
values.append(GetGebcoBathymetryValueFor(vertex.CoordinateValue, 3857, bathymetryData))
fmModel.Bathymetry.SetValues(values)
First we declare a variable called values as a list. Then we continue by looping trough the vertices of the grid and calling the function "GetGebcoBathymetryValueFor" to get the depth value on the location of the vertex and adding it to the list. Finally we add the bathymetry values to the model by using the SetValues method on the bathymetry property of the model. This should result in the following :
