Code in Jupyter Notebook

# Importing Libraries 
import gempy as gp 
import gempy_viewer as gpv 

data_path = 'https://raw.githubusercontent.com/cgre-aachen/gempy_data/master/'
# Create a GeoModel instance
data = gp.create_geomodel(project_name='horizontal', extent=[0, 1000, 0, 1000, 0, 1000], refinement=6, importer_helper=gp.data.ImporterHelper(
        path_to_orientations=data_path + "/data/input_data/jan_models/model1_orientations.csv",
        path_to_surface_points=data_path + "/data/input_data/jan_models/model1_surface_points.csv"
    )
)
# Map geological series to surfaces
gp.map_stack_to_surfaces(gempy_model=data, mapping_object={"Strat_Series": ('rock2', 'rock1')})
# This is generating 3 layers individually(see the result in FIGURE 1 listed in this repository)

# Compute the geological model
gp.compute_model(data)
geo_data = data

# Display results
# Plot the result of the model in the x 
gpv.plot_2d(geo_data, direction=['x'], show_data=True, show_boundaries=False)
# Generated Data can be viewed in FIGURE 2 listed in this repository

# Plot the result of the model in the y 
gpv.plot_2d(geo_data, direction=['y'], show_data=True, show_boundaries=False)
# Generated Data can be viewed in FIGURE 3 listed in this repository