A legacy GM-SYS 2D model file reader coded in Python
This Python file is used to read GM-SYS legacy model files:
- The
sur
file contains model bodies' geometries. - The
blk
file contains model bodies' physical properties (density, magnetic susceptibility, and remanent magnetisation). - The
ecs
file contains model origin, azimuth, and coordinate system information. - The
wel
file contains wells used in the model. - The
grv
file contains the gravity data used to make the model. - The
mag
file contains magnetic observation data.
ecs
file reader is optimised for New Zealand use (it reads NZTM and NZMG coordinate systems). You may want to tweak the gmsys_reader.py
to suit your local coordinate system.
Download the gmsys_reader.py
into a folder containing the legacy GM-SYS files, and import it onto your project by typing from gmsys_reader import *
. BLK and SUR files will be returned as a Pandas DataFrame, while the other files will be read as a mix of NumPy array and Pandas DataFrame. Please see sample_rangitoto_plot.ipynb
for example.
The model used in the example comes from Fig. 5C of my publication Luthfian et al. (2023). The article and figures are licensed under CC BY 4.0 DEED Creative Commons license. You are free to share, copy, redistribute, adapt, and build upon the material for any commercial or non-commercial purpose, with attribution and no additional restrictions applicable.