To add the bathymetry, we will employ the Gebco data set, as used in Delft dashboard. The following code retrieves the bathymetry data from the Gebco data set using the extent of our grid. The coordinate system EPSG code is required because the data set uses a specific (different) coordinate system.
bathymetryData = GetGebcoBathymetryData(xOffset, yOffset, xOffset + gridWidth, yOffset + gridHeigth, EPSGCode)
Next, the cells that are located on land will be removed. This can be determined from the bathymetry data. The CleanupLandCells method removes all land cells that have a z coordinate larger than a certain value. In our case, above 0:
CleanupLandCells(fmModel, bathymetryData, 0)
The grid should now look 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 :
Note: The icon 
indicates functions that you might find interesting to see how they have been built. You can do this by opening the corresponding library (another normal script python file on its own) in the toolbox and checking the code which defines the function or method you are interested in.
