gempy

name	gempy
description	3D structural geological modeling using implicit methods. Create geological models with faults, folds, and unconformities from surface points and orientations. Use when Claude needs to: (1) Build 3D geological models from surface contacts and orientations, (2) Model faults, unconformities, or intrusions, (3) Compute and visualize subsurface geology, (4) Export models to VTK or numpy arrays, (5) Generate gravity forward models, (6) Create cross-sections or 3D visualizations.
version	1.0.0
author	Geoscience Skills
license	MIT
tags	["Geological Modelling","3D","Implicit Surfaces","Faults","Structural Geology"]
dependencies	["gempy>=2.3.0","numpy","pyvista"]
complements	["gemgis","loopstructural","pyvista"]
workflow_role	modelling

GemPy - 3D Geological Modelling

Quick Reference

import gempy as gp

# Create model
geo_model = gp.create_geomodel(
    project_name='Model',
    extent=[0, 1000, 0, 1000, 0, 500],  # [xmin, xmax, ymin, ymax, zmin, zmax]
    resolution=[50, 50, 25]
)

# Add surface points and orientations
gp.add_surface_points(geo_model, x=[100, 500, 900], y=[500, 500, 500],
                      z=[400, 350, 400], surface='TopFormation')
gp.add_orientations(geo_model, x=[500], y=[500], z=[375],
                    pole_vector=[0, 0, 1], surface='TopFormation')

# Compute and visualize
gp.set_interpolator(geo_model)
sol = gp.compute_model(geo_model)
gp.plot_2d(geo_model, cell_number=[25], direction='y')

Key Classes

Class	Purpose
GeoModel	Main container - holds all model data
StructuralFrame	Manages geological relationships
Grid	Computation mesh for interpolation
Solutions	Model results (lithology, scalar fields)

Essential Operations