def METHOD_NAME():
    resolution = [30, 20, 50]
    extent = np.array([729550.0, 751500.0, 1913500.0, 1923650.0, -50, 800.0])
    path_interf = data_path + "/2018_interf.csv"
    path_orient = data_path + "/2018_orient_clust_n_init5_0.csv"
    path_dem = data_path + "/_cropped_DEM_coarse.tif"
    geo_model = gp.create_model('Alesmodel')
    gp.init_data(geo_model, extent=extent, resolution=resolution,
                 path_i=path_interf,
                 path_o=path_orient)
    geo_model.set_topography(source='gdal', filepath=path_dem)
    p2d = gp.plot_2d(geo_model, section_names=['topography'], show_topography=True,
                     kwargs_topography={'hillshade': False, 'fill_contour': False})
    plt.show()
    p2d = gp.plot_2d(geo_model, section_names=['topography'], show_topography=True,
                     kwargs_topography={'hillshade': False, 'fill_contour': True})
    plt.show()
    p2d = gp.plot_2d(geo_model, section_names=['topography'], show_topography=True,
                     kwargs_topography={'hillshade': True, 'fill_contour': True})
    plt.show()
    if False:
        gp.map_stack_to_surfaces(geo_model, {'fault_left'     : ('fault_left'),
                                             'fault_right'    : ('fault_right'),
                                             'fault_lr'       : ('fault_lr'),
                                             'Trias_Series'   : ('TRIAS', 'LIAS'),
                                             'Carbon_Series'  : ('CARBO'),
                                             'Basement_Series': ('basement')}, remove_unused_series=True)
        geo_model.set_is_fault(['fault_right', 'fault_left', 'fault_lr'], change_color=True)
        gp.set_interpolator(geo_model,
                            output=['geology'], compile_aesara=True,
                            aesara_optimizer='fast_run', dtype='float64',
                            verbose=[])
        gp.compute_model(geo_model, compute_mesh=True)
        geo_model._grid.regular_grid.set_topography_mask(geo_model._grid.topography)
    gpv = gp.plot.plot_3d(geo_model,
                          plotter_type='basic', off_screen=True,
                          show_topography=True,
                          show_scalar=False,
                          show_lith=True,
                          show_surfaces=False,
                          kwargs_plot_structured_grid={'opacity'   : 1,
                                                       'show_edges': False},
                          ve=10,
                          image=True,
                          kwargs_plot_topography={'scalars': 'topography'})
    #
    # gpv.p.set_scale(zscale=10)
    #
    # img = gpv.p.show(screenshot=True)
    # plt.imshow(img[1])
    # plt.show()