.. DO NOT EDIT.
.. THIS FILE WAS AUTOMATICALLY GENERATED BY SPHINX-GALLERY.
.. TO MAKE CHANGES, EDIT THE SOURCE PYTHON FILE:
.. "examples/real/Claudius.py"
.. LINE NUMBERS ARE GIVEN BELOW.

.. only:: html

    .. note::
        :class: sphx-glr-download-link-note

        :ref:`Go to the end <sphx_glr_download_examples_real_Claudius.py>`
        to download the full example code

.. rst-class:: sphx-glr-example-title

.. _sphx_glr_examples_real_Claudius.py:


Claudius
========

.. GENERATED FROM PYTHON SOURCE LINES 8-18

.. code-block:: Python

    import sys, os

    # Importing gempy
    import gempy as gp
    import gempy_viewer as gpv

    # Aux imports
    import numpy as np
    import pandas as pn








.. GENERATED FROM PYTHON SOURCE LINES 19-25

Loading data from repository:
~~~~~~~~~~~~~~~~~~~~~~~~~~~~~

With pandas we can do it directly from the web and with the right args
we can directly tidy the data in gempy style:


.. GENERATED FROM PYTHON SOURCE LINES 27-58

.. code-block:: Python


    data_path = os.path.abspath('../../data/input_data/Claudius')

    reduce_data_by = 30

    dfs = []
    for letter in 'ABCD':
        dfs.append(
            pn.read_csv(
                filepath_or_buffer=f"{data_path}/{letter}Points.csv",
                sep=';',
                names=['X', 'Y', 'Z', 'surface', 'cutoff'],
                header=0
            )[::reduce_data_by]
        )

    # Add fault:
    dfs.append(
        pn.read_csv(
            filepath_or_buffer=f"{data_path}/Fault.csv",
            names=['X', 'Y', 'Z', 'surface'],
            header=0,
            sep=','
        )
    )

    surface_points = pn.concat(dfs, sort=True)
    surface_points['surface'] = surface_points['surface'].astype('str')
    surface_points.reset_index(inplace=True, drop=False)
    surface_points.tail()






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <div>
    <style scoped>
        .dataframe tbody tr th:only-of-type {
            vertical-align: middle;
        }

        .dataframe tbody tr th {
            vertical-align: top;
        }

        .dataframe thead th {
            text-align: right;
        }
    </style>
    <table border="1" class="dataframe">
      <thead>
        <tr style="text-align: right;">
          <th></th>
          <th>index</th>
          <th>X</th>
          <th>Y</th>
          <th>Z</th>
          <th>cutoff</th>
          <th>surface</th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <th>755</th>
          <td>88</td>
          <td>551099.250977</td>
          <td>7.817652e+06</td>
          <td>-10466.863281</td>
          <td>NaN</td>
          <td>Claudius_fault</td>
        </tr>
        <tr>
          <th>756</th>
          <td>89</td>
          <td>551160.807495</td>
          <td>7.817503e+06</td>
          <td>-10356.463867</td>
          <td>NaN</td>
          <td>Claudius_fault</td>
        </tr>
        <tr>
          <th>757</th>
          <td>90</td>
          <td>551131.898438</td>
          <td>7.817659e+06</td>
          <td>-10383.323242</td>
          <td>NaN</td>
          <td>Claudius_fault</td>
        </tr>
        <tr>
          <th>758</th>
          <td>91</td>
          <td>551164.412476</td>
          <td>7.817654e+06</td>
          <td>-10299.957031</td>
          <td>NaN</td>
          <td>Claudius_fault</td>
        </tr>
        <tr>
          <th>759</th>
          <td>92</td>
          <td>551197.192139</td>
          <td>7.817647e+06</td>
          <td>-10216.820312</td>
          <td>NaN</td>
          <td>Claudius_fault</td>
        </tr>
      </tbody>
    </table>
    </div>
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 59-61

.. code-block:: Python

    surface_points.dtypes





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    index        int64
    X          float64
    Y          float64
    Z          float64
    cutoff     float64
    surface     object
    dtype: object



.. GENERATED FROM PYTHON SOURCE LINES 62-64

How many points are per surface


.. GENERATED FROM PYTHON SOURCE LINES 66-68

.. code-block:: Python

    surface_points.groupby('surface').count()






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <div>
    <style scoped>
        .dataframe tbody tr th:only-of-type {
            vertical-align: middle;
        }

        .dataframe tbody tr th {
            vertical-align: top;
        }

        .dataframe thead th {
            text-align: right;
        }
    </style>
    <table border="1" class="dataframe">
      <thead>
        <tr style="text-align: right;">
          <th></th>
          <th>index</th>
          <th>X</th>
          <th>Y</th>
          <th>Z</th>
          <th>cutoff</th>
        </tr>
        <tr>
          <th>surface</th>
          <th></th>
          <th></th>
          <th></th>
          <th></th>
          <th></th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <th>0</th>
          <td>167</td>
          <td>167</td>
          <td>167</td>
          <td>167</td>
          <td>167</td>
        </tr>
        <tr>
          <th>250</th>
          <td>167</td>
          <td>167</td>
          <td>167</td>
          <td>167</td>
          <td>167</td>
        </tr>
        <tr>
          <th>330</th>
          <td>166</td>
          <td>166</td>
          <td>166</td>
          <td>166</td>
          <td>166</td>
        </tr>
        <tr>
          <th>60</th>
          <td>167</td>
          <td>167</td>
          <td>167</td>
          <td>167</td>
          <td>167</td>
        </tr>
        <tr>
          <th>Claudius_fault</th>
          <td>93</td>
          <td>93</td>
          <td>93</td>
          <td>93</td>
          <td>0</td>
        </tr>
      </tbody>
    </table>
    </div>
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 69-71

Now we do the same with the orientations:


.. GENERATED FROM PYTHON SOURCE LINES 73-92

.. code-block:: Python

    dfs = []

    for surf in ['0', '330']:
        o = pn.read_csv(
            filepath_or_buffer=f"{data_path}/Dips.csv",
            sep=';',
            names=['X', 'Y', 'Z', 'G_x', 'G_y', 'G_z', '-'],
            header=1
        )

        # Orientation needs to belong to a surface. This is mainly to categorize to which series belong and to
        # use the same color
        o['surface'] = surf
        dfs.append(o)
    orientations = pn.concat(dfs, sort=True)
    orientations.reset_index(inplace=True, drop=False)

    orientations.tail()






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <div>
    <style scoped>
        .dataframe tbody tr th:only-of-type {
            vertical-align: middle;
        }

        .dataframe tbody tr th {
            vertical-align: top;
        }

        .dataframe thead th {
            text-align: right;
        }
    </style>
    <table border="1" class="dataframe">
      <thead>
        <tr style="text-align: right;">
          <th></th>
          <th>index</th>
          <th>-</th>
          <th>G_x</th>
          <th>G_y</th>
          <th>G_z</th>
          <th>X</th>
          <th>Y</th>
          <th>Z</th>
          <th>surface</th>
        </tr>
      </thead>
      <tbody>
        <tr>
          <th>43</th>
          <td>19</td>
          <td>0.976692</td>
          <td>0.187195</td>
          <td>0.138165</td>
          <td>-0.972558</td>
          <td>550989.3105</td>
          <td>7817210.527</td>
          <td>-9782.967773</td>
          <td>330</td>
        </tr>
        <tr>
          <th>44</th>
          <td>20</td>
          <td>0.248187</td>
          <td>-0.080561</td>
          <td>-0.043063</td>
          <td>-0.995819</td>
          <td>550939.3105</td>
          <td>7821227.309</td>
          <td>-9958.425781</td>
          <td>330</td>
        </tr>
        <tr>
          <th>45</th>
          <td>21</td>
          <td>0.649480</td>
          <td>-0.161328</td>
          <td>0.075208</td>
          <td>-0.984031</td>
          <td>549276.8105</td>
          <td>7820682.980</td>
          <td>-9985.125977</td>
          <td>330</td>
        </tr>
        <tr>
          <th>46</th>
          <td>22</td>
          <td>0.050025</td>
          <td>-0.012103</td>
          <td>-0.153309</td>
          <td>-0.988104</td>
          <td>548976.8105</td>
          <td>7820345.121</td>
          <td>-9974.265625</td>
          <td>330</td>
        </tr>
        <tr>
          <th>47</th>
          <td>23</td>
          <td>0.759421</td>
          <td>0.369490</td>
          <td>-0.187053</td>
          <td>-0.910213</td>
          <td>549764.3105</td>
          <td>7820457.738</td>
          <td>-9901.208984</td>
          <td>330</td>
        </tr>
      </tbody>
    </table>
    </div>
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 93-95

.. code-block:: Python

    orientations.dtypes





.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    index        int64
    -          float64
    G_x        float64
    G_y        float64
    G_z        float64
    X          float64
    Y          float64
    Z          float64
    surface     object
    dtype: object



.. GENERATED FROM PYTHON SOURCE LINES 96-104

Data initialization:
~~~~~~~~~~~~~~~~~~~~

Suggested size of the axis-aligned modeling box: Origin: 548800 7816600
-8400 Maximum: 552500 7822000 -11010

Suggested resolution: 100m x 100m x -90m (grid size 38 x 55 x 30)


.. GENERATED FROM PYTHON SOURCE LINES 106-107

Number of voxels:

.. GENERATED FROM PYTHON SOURCE LINES 107-159

.. code-block:: Python

    np.array([38, 55, 30]).prod()

    surface_points_table: gp.data.SurfacePointsTable = gp.data.SurfacePointsTable.from_arrays(
        x=surface_points['X'].values,
        y=surface_points['Y'].values,
        z=surface_points['Z'].values,
        names=surface_points['surface'].values
    )

    orientations_table: gp.data.OrientationsTable = gp.data.OrientationsTable.from_arrays(
        x=orientations['X'].values,
        y=orientations['Y'].values,
        z=orientations['Z'].values,
        G_x=orientations['G_x'].values,
        G_y=orientations['G_y'].values,
        G_z=orientations['G_z'].values,
        names=orientations['surface'].values,
        name_id_map=surface_points_table.name_id_map  # ! Make sure that ids and names are shared
    )

    structural_frame: gp.data.StructuralFrame = gp.data.StructuralFrame.from_data_tables(
        surface_points=surface_points_table,
        orientations=orientations_table
    )

    geo_model: gp.data.GeoModel = gp.create_geomodel(
        project_name='Claudius',
        extent=[548800, 552500, 7816600, 7822000, -11010, -8400],
        resolution=[38, 55, 30],
        refinement=5,
        structural_frame=structural_frame
    )

    group_fault = gp.data.StructuralGroup(
        name='Fault1',
        elements=[geo_model.structural_frame.structural_elements.pop(-2)],
        structural_relation=gp.data.StackRelationType.FAULT,
        fault_relations=gp.data.FaultsRelationSpecialCase.OFFSET_ALL
    )

    geo_model.structural_frame.get_group_by_name("default_formation").elements.pop(-1)

    # Insert the fault group into the structural frame:
    geo_model.structural_frame.insert_group(0, group_fault)

    gp.set_is_fault(
        frame=geo_model.structural_frame,
        fault_groups=[geo_model.structural_frame.get_group_by_name('Fault1')]
    )

    print(geo_model)





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    {'_interpolation_options': InterpolationOptions(kernel_options={'range': 5, 'c_o': 10, 'uni_degree': 1, 'i_res': 4, 'gi_res': 2, 'number_dimensions': 3, 'kernel_function': <AvailableKernelFunctions.cubic: KernelFunction(base_function=<function cubic_function at 0x7f857a0a1b40>, derivative_div_r=<function cubic_function_p_div_r at 0x7f857a0a1bd0>, second_derivative=<function cubic_function_a at 0x7f857a0a1c60>, consume_sq_distance=False)>, 'kernel_solver': <Solvers.DEFAULT: 1>, 'compute_condition_number': False, 'optimizing_condition_number': False, 'condition_number': None}, _number_octree_levels=5, current_octree_level=0, block_solutions_type=BlockSolutionType.OCTREE, compute_scalar_gradient=False, mesh_extraction=True, mesh_extraction_masking_options=MeshExtractionMaskingOptions.INTERSECT, mesh_extraction_fancy=True, debug=True, debug_water_tight=False, sigmoid_slope=50000, evaluation_chunk_size=50000, _number_octree_levels_surface=4, _model_name=None),
     'grid': <gempy.core.data.grid.Grid object at 0x7f8561b42260>,
     'meta': GeoModelMeta(name='Claudius',
                          creation_date=None,
                          last_modification_date=None,
                          owner=None),
     'structural_frame': StructuralFrame(
            structural_groups=[
    StructuralGroup(
            name=Fault1,
            structural_relation=StackRelationType.FAULT,
            elements=[
    Element(
            name=Claudius_fault,
            color=#527682,
            is_active=True
    )
    ]
    ),
    StructuralGroup(
            name=default_formation,
            structural_relation=StackRelationType.ERODE,
            elements=[
    Element(
            name=0,
            color=#015482,
            is_active=True
    ),
    Element(
            name=250,
            color=#9f0052,
            is_active=True
    ),
    Element(
            name=330,
            color=#ffbe00,
            is_active=True
    ),
    Element(
            name=60,
            color=#728f02,
            is_active=True
    )
    ]
    )
    ],
            fault_relations=
    [[False,  True],
     [False, False]],
    ,
     'transform': {'_is_default_transform': False,
     'position': array([ -550658.0605    , -7819213.97425586,     9824.6748045 ]),
     'rotation': array([0., 0., 0.]),
     'scale': array([9.56397917e-05, 9.56397917e-05, 9.56397917e-05])}}




.. GENERATED FROM PYTHON SOURCE LINES 160-163

We are going to increase the smoothness (nugget) of the data to increase
the conditional number of the matrix:


.. GENERATED FROM PYTHON SOURCE LINES 165-167

.. code-block:: Python

    gp.modify_surface_points(geo_model, nugget=0.01)






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">

            <table>
              <tr><td>Structural Groups:</td><td>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Fault1</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.FAULT</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #527682;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Claudius_fault</td></tr>
        </table>
            </td></tr>
        </table>
            <br>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>default_formation</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #015482;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>0</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #9f0052;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>250</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #ffbe00;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>330</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #728f02;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>60</td></tr>
        </table>
            </td></tr>
        </table>
            </td></tr>
              <tr><td>Fault Relations:</td><td><table style="border-collapse: collapse; table-layout: fixed;"><th></th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Fault1</th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">default_fo...</th><tr><th>Fault1</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></td></tr><tr><th>default_formation</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr></table></td></tr>
              <tr><td></td><td>
            <table>
              <tr>
                <td><div style="display: inline-block; background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></div> True</td>
                <td><div style="display: inline-block; background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></div> False</td>
              </tr>
            </table>
            </td></tr>
            </table>
        
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 168-171

Also the original poles are pointing downwards. We can change the
direction by calling the following:


.. GENERATED FROM PYTHON SOURCE LINES 173-175

.. code-block:: Python

    gp.modify_orientations(geo_model, polarity=-1)






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">

            <table>
              <tr><td>Structural Groups:</td><td>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Fault1</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.FAULT</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #527682;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Claudius_fault</td></tr>
        </table>
            </td></tr>
        </table>
            <br>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>default_formation</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #015482;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>0</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #9f0052;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>250</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #ffbe00;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>330</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #728f02;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>60</td></tr>
        </table>
            </td></tr>
        </table>
            </td></tr>
              <tr><td>Fault Relations:</td><td><table style="border-collapse: collapse; table-layout: fixed;"><th></th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Fault1</th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">default_fo...</th><tr><th>Fault1</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></td></tr><tr><th>default_formation</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr></table></td></tr>
              <tr><td></td><td>
            <table>
              <tr>
                <td><div style="display: inline-block; background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></div> True</td>
                <td><div style="display: inline-block; background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></div> False</td>
              </tr>
            </table>
            </td></tr>
            </table>
        
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 176-180

We need an orientation per series/fault. The faults does not have
orientation so the easiest is to create an orientation from the surface
points availablle:


.. GENERATED FROM PYTHON SOURCE LINES 180-194

.. code-block:: Python


    element = geo_model.structural_frame.get_element_by_name("Claudius_fault")
    new_orientations: gp.data.OrientationsTable = gp.create_orientations_from_surface_points_coords(
        xyz_coords=element.surface_points.xyz
    )
    gp.add_orientations(
        geo_model=geo_model,
        x=new_orientations.data['X'],
        y=new_orientations.data['Y'],
        z=new_orientations.data['Z'],
        pole_vector=new_orientations.grads,
        elements_names="Claudius_fault"
    )






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">

            <table>
              <tr><td>Structural Groups:</td><td>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Fault1</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.FAULT</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #527682;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Claudius_fault</td></tr>
        </table>
            </td></tr>
        </table>
            <br>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>default_formation</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #015482;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>0</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #9f0052;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>250</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #ffbe00;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>330</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #728f02;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>60</td></tr>
        </table>
            </td></tr>
        </table>
            </td></tr>
              <tr><td>Fault Relations:</td><td><table style="border-collapse: collapse; table-layout: fixed;"><th></th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Fault1</th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">default_fo...</th><tr><th>Fault1</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></td></tr><tr><th>default_formation</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr></table></td></tr>
              <tr><td></td><td>
            <table>
              <tr>
                <td><div style="display: inline-block; background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></div> True</td>
                <td><div style="display: inline-block; background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></div> False</td>
              </tr>
            </table>
            </td></tr>
            </table>
        
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 195-197

.. code-block:: Python

    gpv.plot_2d(geo_model, direction='y')




.. image-sg:: /examples/real/images/sphx_glr_Claudius_001.png
   :alt: Cell Number: mid Direction: y
   :srcset: /examples/real/images/sphx_glr_Claudius_001.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f8517306080>



.. GENERATED FROM PYTHON SOURCE LINES 198-200

We will need to separate with surface belong to each series:


.. GENERATED FROM PYTHON SOURCE LINES 200-209

.. code-block:: Python


    gp.map_stack_to_surfaces(
        gempy_model=geo_model,
        mapping_object={
            'Default series': ('0', '60', '250'),
            'Fault': 'Claudius_fault',
            'Uncomformity': '330',
        }
    )





.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">

            <table>
              <tr><td>Structural Groups:</td><td>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Default series</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #015482;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>0</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #728f02;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>60</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #9f0052;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>250</td></tr>
        </table>
            </td></tr>
        </table>
            <br>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Fault</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #527682;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Claudius_fault</td></tr>
        </table>
            </td></tr>
        </table>
            <br>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Uncomformity</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #ffbe00;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>330</td></tr>
        </table>
            </td></tr>
        </table>
            </td></tr>
              <tr><td>Fault Relations:</td><td><table style="border-collapse: collapse; table-layout: fixed;"><th></th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Default se...</th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Fault</th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Uncomformi...</th><tr><th>Default series</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr><tr><th>Fault</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr><tr><th>Uncomformity</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr></table></td></tr>
              <tr><td></td><td>
            <table>
              <tr>
                <td><div style="display: inline-block; background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></div> True</td>
                <td><div style="display: inline-block; background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></div> False</td>
              </tr>
            </table>
            </td></tr>
            </table>
        
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 210-212

So far we did not specify which series/faults are actula faults:


.. GENERATED FROM PYTHON SOURCE LINES 214-221

.. code-block:: Python

    gp.set_is_fault(
        frame=geo_model.structural_frame,
        fault_groups=[geo_model.structural_frame.get_group_by_name('Fault')]
    )

    geo_model.structural_frame






.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">

            <table>
              <tr><td>Structural Groups:</td><td>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Default series</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #015482;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>0</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #728f02;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>60</td></tr>
        </table>
            <br>
        <table width="50%" style="border-left:15px solid #9f0052;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>250</td></tr>
        </table>
            </td></tr>
        </table>
            <br>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Fault</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.FAULT</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #527682;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>Claudius_fault</td></tr>
        </table>
            </td></tr>
        </table>
            <br>
        <table style="border-left:1.2px solid black;>
          <tr><th colspan="2"><b>StructuralGroup:</b></th></tr>
          <tr><td>Name:</td><td>Uncomformity</td></tr>
          <tr><td>Structural Relation:</td><td>StackRelationType.ERODE</td></tr>
          <tr><td>Elements:</td><td>
        <table width="50%" style="border-left:15px solid #ffbe00;">
          <tr><th colspan="2"><b>StructuralElement:</b></th></tr>
          <tr><td>Name:</td><td>330</td></tr>
        </table>
            </td></tr>
        </table>
            </td></tr>
              <tr><td>Fault Relations:</td><td><table style="border-collapse: collapse; table-layout: fixed;"><th></th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Default se...</th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Fault</th><th style="transform: rotate(-35deg); height:150px; vertical-align: bottom; text-align: center;">Uncomformi...</th><tr><th>Default series</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr><tr><th>Fault</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></td></tr><tr><th>Uncomformity</th><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td><td style="background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></td></tr></table></td></tr>
              <tr><td></td><td>
            <table>
              <tr>
                <td><div style="display: inline-block; background-color: #527682; width: 20px; height: 20px; border: 1px solid black;"></div> True</td>
                <td><div style="display: inline-block; background-color: #FFB6C1; width: 20px; height: 20px; border: 1px solid black;"></div> False</td>
              </tr>
            </table>
            </td></tr>
            </table>
        
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 222-232

.. code-block:: Python

    geo_model.interpolation_options.kernel_options.range = 1
    gp.compute_model(
        geo_model,
        gp.data.GemPyEngineConfig(
            backend=gp.data.AvailableBackends.numpy,
            use_gpu=False,
            dtype='float64'
        )
    )





.. rst-class:: sphx-glr-script-out

 .. code-block:: none

    Setting Backend To: AvailableBackends.numpy


.. raw:: html

    <div class="output_subarea output_html rendered_html output_result">
    <b>Solutions:</b> 5 Octree Levels, 5 DualContouringMeshes
    </div>
    <br />
    <br />

.. GENERATED FROM PYTHON SOURCE LINES 233-237

.. code-block:: Python

    sect = ['mid']

    gpv.plot_2d(geo_model, cell_number=sect, series_n=1, show_scalar=True, direction='x')




.. image-sg:: /examples/real/images/sphx_glr_Claudius_002.png
   :alt: Cell Number: mid Direction: x
   :srcset: /examples/real/images/sphx_glr_Claudius_002.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f851715caf0>



.. GENERATED FROM PYTHON SOURCE LINES 238-240

.. code-block:: Python

    gpv.plot_2d(geo_model, cell_number=sect, show_data=True, direction='x')




.. image-sg:: /examples/real/images/sphx_glr_Claudius_003.png
   :alt: Cell Number: mid Direction: x
   :srcset: /examples/real/images/sphx_glr_Claudius_003.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f851715ef50>



.. GENERATED FROM PYTHON SOURCE LINES 241-245

.. code-block:: Python

    gpv.plot_2d(geo_model, cell_number=[28], series_n=0, direction='y', show_scalar=True)
    gpv.plot_2d(geo_model, cell_number=[28], series_n=1, direction='y', show_scalar=True)
    gpv.plot_2d(geo_model, cell_number=[28], series_n=2, direction='y', show_scalar=True)




.. rst-class:: sphx-glr-horizontal


    *

      .. image-sg:: /examples/real/images/sphx_glr_Claudius_004.png
         :alt: Cell Number: 28 Direction: y
         :srcset: /examples/real/images/sphx_glr_Claudius_004.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/real/images/sphx_glr_Claudius_005.png
         :alt: Cell Number: 28 Direction: y
         :srcset: /examples/real/images/sphx_glr_Claudius_005.png
         :class: sphx-glr-multi-img

    *

      .. image-sg:: /examples/real/images/sphx_glr_Claudius_006.png
         :alt: Cell Number: 28 Direction: y
         :srcset: /examples/real/images/sphx_glr_Claudius_006.png
         :class: sphx-glr-multi-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f850ccf29e0>



.. GENERATED FROM PYTHON SOURCE LINES 246-248

.. code-block:: Python

    gpv.plot_2d(geo_model, cell_number=[28], show_data=True, direction='y')




.. image-sg:: /examples/real/images/sphx_glr_Claudius_007.png
   :alt: Cell Number: 28 Direction: y
   :srcset: /examples/real/images/sphx_glr_Claudius_007.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gempy_viewer.modules.plot_2d.visualization_2d.Plot2D object at 0x7f850cd782e0>



.. GENERATED FROM PYTHON SOURCE LINES 249-252

.. code-block:: Python


    # sphinx_gallery_thumbnail_number = 8
    gpv.plot_3d(geo_model, show_lith=True, show_data=True, show_boundaries=True)



.. image-sg:: /examples/real/images/sphx_glr_Claudius_008.png
   :alt: Claudius
   :srcset: /examples/real/images/sphx_glr_Claudius_008.png
   :class: sphx-glr-single-img


.. rst-class:: sphx-glr-script-out

 .. code-block:: none


    <gempy_viewer.modules.plot_3d.vista.GemPyToVista object at 0x7f850c35cac0>




.. rst-class:: sphx-glr-timing

   **Total running time of the script:** (1 minutes 2.331 seconds)


.. _sphx_glr_download_examples_real_Claudius.py:

.. only:: html

  .. container:: sphx-glr-footer sphx-glr-footer-example

    .. container:: sphx-glr-download sphx-glr-download-jupyter

      :download:`Download Jupyter notebook: Claudius.ipynb <Claudius.ipynb>`

    .. container:: sphx-glr-download sphx-glr-download-python

      :download:`Download Python source code: Claudius.py <Claudius.py>`


.. only:: html

 .. rst-class:: sphx-glr-signature

    `Gallery generated by Sphinx-Gallery <https://sphinx-gallery.github.io>`_