def test_set_orientations():
# Importing the data from CSV-files and setting extent and resolution
geo_data = gp.create_data(
extent=[0, 2000, 0, 2000, 0, 2000],
resolution=[50, 50, 50],
path_o=input_path +
'/input_data/tut_chapter1/simple_fault_model_orientations.csv',
path_i=input_path +
'/input_data/tut_chapter1/simple_fault_model_points.csv')
gp.get_data(geo_data)
# Assigning series to formations as well as their order (timewise)
gp.map_stack_to_surfaces(geo_data, {
"Fault_Series": 'Main_Fault',
"Strat_Series": ('Sandstone_2', 'Siltstone')
})
geo_data._orientations.create_orientation_from_surface_points(
geo_data.surface_points, [0, 1, 2])
gp.set_orientation_from_surface_points(geo_data, [0, 1, 2])
# Also the original poles are pointing downwards. We can change the # direction by calling the following: # # %% geo_model.modify_orientations(geo_model.orientations.df.index, polarity=-1).df.tail() # %% # 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: # # %% fault_idx = geo_model.surface_points.df.index[geo_model.surface_points.df['surface'] == 'Claudius_fault'] gp.set_orientation_from_surface_points(geo_model, fault_idx).df.tail() # %% # Now we can see how the data looks so far: # # %% geo_model.surfaces # %% gp.plot_2d(geo_model, direction='y') # %% # By default all surfaces belong to one unique series. #
