def METHOD_NAME():
    METHOD_NAME = RasterModelGrid((10, 10), xy_spacing=10.0)
    METHOD_NAME.set_closed_boundaries_at_grid_edges(True, True, True, True)
    z = METHOD_NAME.add_zeros("node", "topographic__elevation")
    METHOD_NAME.add_zeros("node", "soil__depth")
    z += METHOD_NAME.x_of_node.copy() + METHOD_NAME.y_of_node.copy()
    z[25] -= 40
    z[35] -= 40
    z[26] -= 40
    z[36] -= 40
    z[24] -= 40
    z[34] -= 40
    return METHOD_NAME