|
1 | 1 | # Integrating diagenetic equations using Python
|
2 | 2 |
|
3 |
| -This repo was created as an attempt to reproduce the plots shown at the kickoff of the AstroTOM ("Turing or Milankovitch") project by Niklas Hohmann, from his Matlab scripts (available at [github.com/MindTheGap-ERC/LMA-Matlab](https://github.com/MindTheGap-ERC/LMA-Matlab)). |
| 3 | +This repo was forked off [Integrating-diagenetic-equations-using-Python](https://github.com/astro-turing/Integrating-diagenetic-equations-using-Python), which is part of the [Astro-turing organisation](https://github.com/astro-turing), created as part of the AstroToM ("Turing or Milankovitch") project, an OpenSSI 2021b project from the Netherlands eScience Center and Utrecht University (UU). |
4 | 4 |
|
5 |
| -AstroTOM is an OpenSSI 2021b project from the Netherlands eScience Center and Utrecht University (UU). |
| 5 | +The diagenetic modelling efforts within AstoToM can be regarded as a precursor to this repo and to [rhythmite](https://github.com/MindTheGap-ERC/rhythmite), while [LMA_utils](https://github.com/MindTheGap-ERC/LMA_utils), [LHeureuxEqs](https://github.com/MindTheGap-ERC/LHeureuxEqs) and [Cross-comparison](https://github.com/MindTheGap-ERC/Cross-comparison) are auxiliary repos. |
6 | 6 |
|
7 |
| -Dr. Emilia Jarochowska (UU) is the lead applicant of this project. |
| 7 | +[LMA-Matlab](https://github.com/MindTheGap-ERC/LMA-Matlab) was the first attempt to reproduce the results from [L'Heureux (2018)](https://onlinelibrary.wiley.com/doi/10.1155/2018/4968315). That repo is coded in MATLAB, while the [original diagenetic model from L'Heureux](https://github.com/astro-turing/Diagenetic_model_LHeureux_2018) was written in FORTRAN. |
| 8 | +[Integrating-diagenetic-equations-using-Python](https://github.com/astro-turing/Integrating-diagenetic-equations-using-Python) was inspired by [LMA-Matlab](https://github.com/MindTheGap-ERC/LMA-Matlab). |
8 | 9 |
|
9 |
| -After replacing central differencing for the gradients in the five diagenetic equations 40-43 from [L'Heureux (2018)](https://www.hindawi.com/journals/geofluids/2018/4968315/) by forward and backward differencing depending on the sign of U and W as a first step and a Fiadeiro-Veronis spatial difference scheme as a second step, it turns out that these equations can be integrated for more than 13.190 years (the full T*) with an implicit or explicit (in time) solver, but not with a simple Eulerian scheme. A Runge-Kutta solver, with an adaptive timestep will, however, suffice. |
10 |
| -After correcting the value of b (5-->5e-4) it turned out that a stable integration is also possible without a Fiadeiro-Veronis scheme. The `main` branch makes use of a constant porosity diffusion coefficient. |
| 10 | +MindTheGap is led by dr. Emilia Jarochowska (UU). |
11 | 11 |
|
12 |
| -Implicit (in time) solvers with use of Jacobians (in functional forms, so without numerical approximations) are available in the `Use_solve_ivp_without_py-pde_wrapper` branch. |
| 12 | +Wide use is made of the [py-pde](https://py-pde.readthedocs.io/en/latest/) package. |
13 | 13 |
|
14 |
| -Wide use is made of the [py-pde](https://py-pde.readthedocs.io/en/latest/) package, especially in the `main` branch. |
| 14 | +The porosity diffusion coefficient is held constant. |
15 | 15 |
|
16 | 16 | ## Installing and using
|
17 | 17 | To run this code, you need `git` and `conda` or `pip` to install .
|
18 | 18 | ```
|
19 |
| -git clone [email protected]:astro-turing/Integrating-diagenetic-equations-using-Python.git |
| 19 | +git clone [email protected]:MindTheGap-ERC/reactive-transport-model-for-limestone-marl-sequences.git |
20 | 20 | ```
|
21 | 21 | or
|
22 | 22 | ```
|
23 |
| -git clone https://github.com/astro-turing/Integrating-diagenetic-equations-using-Python.git |
| 23 | +git clone https://github.com/MindTheGap-ERC/reactive-transport-model-for-limestone-marl-sequences.git |
24 | 24 | ```
|
25 | 25 | Next,
|
26 | 26 | ```
|
27 |
| -cd Integrating-diagenetic-equations-using-Python |
| 27 | +cd reactive-transport-model-for-limestone-marl-sequences |
28 | 28 | git switch release_v1.0.0
|
29 | 29 | pipenv install
|
30 | 30 | ```
|
|
0 commit comments