Update README

README.md

@@ -4,15 +4,15 @@
 [![DOI](https://zenodo.org/badge/299357364.svg)](https://zenodo.org/badge/latestdoi/299357364)
 
 
-This repository contains Pseudo-Transient (PT) routines resolving chemical reactions coupled to fluid flow in viscously defroming solid pourous matrix, so-called Hydro-Mechanical-Chemical (HMC) coupling. Example of such multi-physical processes to understand, e.g, the formation of olivine veins by dehydration of ductile serpentinite ([Schmalholz et al., 2023 - submitted](), [Schmalholz et al., 2022](https://doi.org/10.1002/essoar.10512291.2)), or the brucite-periclase reactions [(Schmalholz et al., 2020)](https://doi.org/10.1029/2020GC009351).
+This repository contains Pseudo-Transient (PT) routines resolving chemical reactions coupled to fluid flow in viscously defroming solid pourous matrix, so-called Hydro-Mechanical-Chemical (HMC) coupling. Example of such multi-physical processes to understand, e.g, the formation of olivine veins by dehydration of ductile serpentinite ([Schmalholz et al., 2023](https://doi.org/10.1029/2023JB026985)), or the brucite-periclase reactions ([Schmalholz et al., 2020](https://doi.org/10.1029/2020GC009351)).
 
 ![Serpentinite dehydration and olivine vein formation during ductile shearing](docs/viz_dehy.png)
 > Serpentinite dehydration and olivine vein formation during ductile shearing
 
 Pseudo-Transient approach relies in using physics-motivated transient terms within differential equations in order to iteratively converge to an accurate solution. The PT HMC routines are written using the [Julia programming language](https://julialang.org) and build upon the high-performance [ParallelStencil.jl](https://github.com/omlins/ParallelStencil.jl) package to enable for optimal execution on graphics processing units (GPUs) and multi-threaded CPUs.
 
 ## Script list
-The Julia scripts are located in the [scripts](scripts) folder which contains the routines for the [(Schmalholz et al., _submitted_)]() study.
+The Julia scripts are located in the [scripts](scripts) folder which contains the routines for the [(Schmalholz et al., 2023)](https://doi.org/10.1029/2023JB026985) study.
 - The main script is [`DeHy.jl`](scripts/DeHy.jl);
 - The visualisation script is [`vizme_DeHy.jl`](scripts/vizme_DeHy.jl).
 
@@ -57,9 +57,7 @@ julia --project DeHy.jl
 ## References
 Schmalholz, S. M., Moulas, E., Plümper, O., Myasnikov, A. V., & Podladchikov, Y. Y. (2020). **2D hydro‐mechanical‐chemical modeling of (De)hydration reactions in deforming heterogeneous rock: The periclase‐brucite model reaction**. Geochemistry, Geophysics, Geosystems, 21, 2020GC009351. [https://doi.org/10.1029/2020GC009351](https://doi.org/10.1029/2020GC009351)
 
-Schmalholz, S. M., Moulas, E., Räss, L., & Müntener, O. (2022).  **Shear-driven formation of olivine veins by dehydration of ductile serpentinite: a numerical study with implications for transient weakening.** ESS Open Archive. [https://doi.org/10.1002/essoar.10512291.2](https://doi.org/10.1002/essoar.10512291.2)
-
-Schmalholz, S. M., Moulas, E., Räss, L., & Müntener, O. (submitted). **Shear-driven formation of olivine veins by dehydration of ductile serpentinite: a numerical study with implications for porosity production and transient weakening.** [_Submitted to JGR_]()
+Schmalholz, S. M., Moulas, E., Räss, L., & Müntener, O. (2023). **Shear-driven formation of olivine veins by dehydration of ductile serpentinite: a numerical study with implications for porosity production and transient weakening.** [https://doi.org/10.1029/2023JB026985](https://doi.org/10.1029/2023JB026985)
 
 [CUDA.jl]: https://github.com/JuliaGPU/CUDA.jl
 [Makie.jl]: https://docs.makie.org/stable/