|
2 | 2 | Reading Simple GemPy model in Subsurface |
3 | 3 | ======================================== |
4 | 4 | """ |
5 | | -# |
6 | | -# import sys |
7 | | -# |
8 | | -# sys.path.insert(0, 'subsurface/') |
9 | | -# import pooch |
10 | | -# import subsurface as ss |
11 | | -# import subsurface.reader.read_netcdf |
12 | | -# |
13 | | -# # %% md |
14 | | -# # Single surface |
15 | | -# # -------------- |
16 | | -# |
17 | | -# # Pull gempy model |
18 | | -# model_files = pooch.retrieve( |
19 | | -# url="https://github.com/cgre-aachen/gempy_data/raw/master/" |
20 | | -# "data/subsurface/example1.zip", |
21 | | -# known_hash="fb79a63eeb874cf0cdca557106c62c67eace23811db5935e57c3448fed7f8978", |
22 | | -# processor=pooch.Unzip() |
23 | | -# ) |
24 | | -# |
25 | | -# # %% |
26 | | -# fname, = [i for i in model_files if "meshes.nc" in i] |
27 | | -# dataset = ss.reader.read_netcdf.read_unstruct(fname) |
28 | | -# |
| 5 | + |
| 6 | +import sys |
| 7 | + |
| 8 | +sys.path.insert(0, 'subsurface/') |
| 9 | +import pooch |
| 10 | +import subsurface as ss |
| 11 | +import subsurface.reader.read_netcdf |
| 12 | + |
| 13 | +# %% md |
| 14 | +# Single surface |
| 15 | +# -------------- |
| 16 | + |
| 17 | +# Pull gempy model |
| 18 | +model_files = pooch.retrieve( |
| 19 | + url="https://github.com/cgre-aachen/gempy_data/raw/master/" |
| 20 | + "data/subsurface/example1.zip", |
| 21 | + known_hash="fb79a63eeb874cf0cdca557106c62c67eace23811db5935e57c3448fed7f8978", |
| 22 | + processor=pooch.Unzip() |
| 23 | +) |
| 24 | + |
29 | 25 | # # %% |
30 | | -# obj = ss.TriSurf(dataset) |
31 | | -# print(obj.mesh.points_attributes_to_dict) |
32 | | -# |
33 | | -# mesh = ss.visualization.to_pyvista_mesh(obj) |
34 | | -# ss.visualization.pv_plot([mesh]) |
35 | | -# |
| 26 | +fname, = [i for i in model_files if "meshes.nc" in i] |
| 27 | +dataset = ss.reader.read_netcdf.read_unstruct(fname) |
| 28 | + |
| 29 | +# %% |
| 30 | +obj = ss.TriSurf(dataset) |
| 31 | +print(obj.mesh.points_attributes_to_dict) |
| 32 | + |
| 33 | +mesh = ss.visualization.to_pyvista_mesh(obj) |
| 34 | +ss.visualization.pv_plot([mesh]) |
| 35 | + |
36 | 36 | # |
37 | 37 | # # %% md |
38 | 38 | # # Four Layers |
|
0 commit comments