[BUG] Computing a model only with sections mess up sorting surfaces with solutions

[Leguark]

Describe the bug Computing a model only with sections mess up sorting surfaces with solutions

For some reason scalar fields at surface point are not in order

To Reproduce Provide detailed steps to reproduce the behavior:

@pytest.fixture(scope='module')
def section_model(one_fault_model_topo_solution):
    geo_model = one_fault_model_topo_solution
    section_dict = {'section_SW-NE': ([250, 250], [1750, 1750], [100, 100]),
                    'section_NW-SE': ([250, 1750], [1750, 250], [100, 100])}
    geo_model.set_section_grid(section_dict)

    geo_model.set_active_grid('sections', reset=False)

    one_fault_model_topo_solution.update_additional_data()
    one_fault_model_topo_solution.update_to_interpolator()
    gp.compute_model(geo_model, sort_surfaces=True)
    return geo_model

Expected behavior After one compute the color of the layers should converge

[Wh20-fly]

Hi, @Leguark .Now can GemPy obtain the 2D section in the NW-SE direction? A few days ago, an author said that GemPy cannot obtain the 2D section in this direction at present, as shown in the picture below. What does this BUG of yours show? Thank you very much.

[Leguark]

No this bug is completely unrelated. Can you make a pull request with a test showing that behavior? If you doing it like this, it is much easier for us to check what is the problem and submit a patch.

To do so:

1. Pull a fresh branch from master
2. Create in gempy/test//test_plotting/test_plot_section.py
3. Add the function:

def test_plot_East-west():
     # Your code replicating the problem here

4. git commit
5. git push to your repo
6. pull request

If you have any doubts on the way just contact one of us.

[Wh20-fly]

Hi, @Leguark .In your comment,the third point "Add the function:",how can I do that,could you explain more detailedly.And the second point means what,how can I make it,and whether I finish the second one or the third one will solve my problem,thank you.

[Wh20-fly]

Hi, @Leguark .Can GemPy plot any direction's 2D section,for example NW-SE,I am in a hurry to get that function,thank you, and please see the upper comment of mine.

[AlexanderJuestel]

Hello,

please see ch1-7a-CrossSections for creating custom section.

I believe you actually worked with those before.