OGSTools 0.4.0 Release Notes
Overview
Supports Python: 3.10, 3.11, 3.12.
First release after integration of ogs6py. The state of ogs6py corresponds to the functionality of v.403 of ogs6py (only light API changes, see Features/Project file)
OGS
OpenGeoSys version: 6.5.3
Migration strategy from ogs6py to ogstools
The recommended strategy is:
- Back up your environment with
pip freeze > yourfile.txt
- Upgrade ogs6py to 0.403 first and resolve issues (see breaking changes).
- Uninstall ogs6py and install ogstools
- Remove
import ogs6py
. Addimport ogstools as ogs
. The formerOGS
becomesProject
and function parameters ofOGS.__init__()
are now with lower case names. See example in Features
API changes
plot
- renamed meshplotlib to plot
- renamed plot function to contourf
- renamed setup.p_min -> setup.vmin, setup.p_max -> setup.vmax
- added several setup options available as kwargs
variables
- renamed propertylib to variables (to align with OpenGeoSys vocabulary)
- presets are now directly accessible as variables
- renamed presets to properties
- renamed Property class to Variable
meshlib
- renamed aggregate function to time_aggregate
- renamed MeshSeries.read to MeshSeries.mesh
Features
ogs6py
- new version with significant parts of the codebase rewritten. Changes visible to the user:
example
before:
import ogstools as ogs
prj = ogs.Project(input_file="mechanics.prj", output_file="old_parameter_add.prj")
prj.add_block(
blocktag="parameter",
parent_xpath="./parameters",
taglist=["name", "type", "value"],
textlist=["density", "Constant", "1900"],
)
prj.write_input()
now:
import ogstools as ogs
prj = ogs.Project(input_file="mechanics.prj", output_file="new_parameter_add.prj")
prj.parameters.add_parameter(name="density", type="Constant", value="1900")
prj.write_input()
- support for combining replacement- and build-method commands
- breaking changes: some methods renamed closer to project file tags:
add_process_variable()
split-up intoadd_secondary_variable()
andadd_process_variable()
geo.add_geom()
->geometry.add_geometry()
,timeloop
->time_loop
etc.
- support for more MPL properties
- several bugfixes and tiny new features
-
integration into ogstools
plot
- Time slice plots can be created easily: filled contour plots of transient data along a sampling line (example)
- Variables now have default symbols (e.g. used for labels in plots)
Mesh and MeshSeries
- Mesh can be created from a shapefile
mesh = ogs.Mesh.read(test_shapefile)
- Mesh can be remeshed with triangle
- MeshSeries has a performant algorithm for integration point tessellation called
MeshSeries.ip_tesselated()
- example for analyzing integration point data - MeshSeries allows multidimensional indexing on ndarrays https://numpy.org/doc/stable/user/basics.indexing.html
import ogstools as ogs
ms = ogs.MeshSeries("filepath/filename_pvd_or_xdmf")
ms.data("darcy_velocity")[-2:, 1:4, :] # shape is(2, 3, 2)
result_mesh = ms[-1]
for mesh in ms:
print(mesh)
- Added function argument to Meshseries.animate functions has more flexible parameter that allow
- transformation of the mesh and
- customization of the plot
- Meshseries has domain_aggregate function (e.g. min/max of a variable per time step)
Project files
- ogs6py added (this version corresponds to https://github.com/joergbuchwald/ogs6py/releases/tag/v.403)
- OGS class is renamed to Project
- Function parameters are with now lower case letters
project_file
is nowoutput_file
Documentation, examples and more
- Documentation got new structure of examples (now all organized under ogstools/examples/)
- pip requirements.txt files (of tested environment) are added for stable test environment reproduction (To be used by developers. Do not use it in your projects!)
Bugfixes
- Several small fixes in plotting functions (visual)
- MeshSeries closes file handle to h5 file after reading is finished
- Dependency compatibility (e.g., remove restriction to matplotlib and relaxing requirements)