Example: 1D Shallow Water Equations Dam Break
Bottom topography
    def bottom(x):
        bump = 5.0 * np.exp(-((x - 800)**2) / 12)
        trough = -3.3 * np.exp(-((x - 900)**2) / 600)
        return bump + trough
relatedh_dam_break_filled_with_velocity_bottomtopo