def METHOD_NAME():
    cwd = os.path.dirname(__file__)
    data_path = cwd + '/../../../examples/'
    geo_model = gp.load_model(r'Tutorial_ch1-8_Onlap_relations',
                              path=data_path + 'data/gempy_models/Tutorial_ch1-8_Onlap_relations',
                              recompile=True)
    geo_model.set_regular_grid([-200, 1000, -500, 500, -1000, 0], [50, 50, 50])
   # geo_model.set_topography(d_z=np.array([-600, -100]))
    s = gp.compute_model(geo_model, compute_mesh=True, debug=False)
    gp.plot.plot_2d(geo_model, cell_number=2)
    gp.plot_2d(geo_model, cell_number=[2],
               regular_grid=geo_model.solutions.mask_matrix_pad[0],
               show_data=True, kwargs_regular_grid={'cmap': 'gray', 'norm': None})
    gp.plot_2d(geo_model, cell_number=[2],
               regular_grid=geo_model.solutions.mask_matrix_pad[1],
               show_data=True, kwargs_regular_grid={'cmap': 'gray', 'norm': None})
    gp.plot_2d(geo_model, cell_number=[2],
               regular_grid=geo_model.solutions.mask_matrix_pad[2],
               show_data=True, kwargs_regular_grid={'cmap': 'gray', 'norm': None})
    gp.plot_2d(geo_model, cell_number=[2],
               regular_grid=geo_model.solutions.mask_matrix_pad[3],
               show_data=True, kwargs_regular_grid={'cmap': 'gray', 'norm': None})
    p3d = gp.plot_3d(geo_model, show_surfaces=True, show_data=True,
                     image=True,
                     kwargs_plot_structured_grid={'opacity': .2})