diff --git a/README.md b/README.md
index 6346799..e220267 100644
--- a/README.md
+++ b/README.md
@@ -1,6 +1,6 @@
 ## PhaceXL 2D
 
-PhaceXL_2D is a script for the generation of finite-thickness interface elements for 2D polycrystal microstructure modelling. It parses a 2D polycrystal microstructure mesh file generated using the [Neper](http://neper.sourceforge.net/) software package and modified using the [Phon](http://dx.doi.org/10.5281/zenodo.16711) cohesive element generator and shrinks the existing grains by a predefined factor to create finite-thickness interface elements at the grain boundaries. If the original 2D polycrystal microstructure mesh is periodic, PhaceXL_2D retains mesh periodicity at the boundary of the modelling domain. The resulting mesh is exported in a format readable by Abaqus™ software and can be used for finite-element analyses.
+PhaceXL_2D is a script for the generation of finite-thickness interface elements for 2D polycrystal microstructure modelling. It parses a 2D polycrystal microstructure mesh file generated using the [Neper](http://neper.sourceforge.net/) software package and modified using the [Phon](http://github.com/KristofferC/Phon) cohesive element generator and shrinks the existing grains by a predefined factor to create finite-thickness interface elements at the grain boundaries. If the original 2D polycrystal microstructure mesh is periodic, PhaceXL_2D retains mesh periodicity at the boundary of the modelling domain. The resulting mesh is exported in a format readable by Abaqus™ software and can be used for finite-element analyses.
 
 ![Exemplary 2D polycrystal mesh segment](docs/example.jpg "Exemplary 2D polycrystal mesh segment with finite-thickness interface elements at grain boundaries")
 
@@ -8,7 +8,7 @@ To generate a 2D polycrystal microstructure mesh with finite-thickness interface
 
 1. Generate a 2D polycrystal microstructure mesh using [Neper](http://neper.sourceforge.net/).
 2. Export the polycrystal mesh as well as the coordinates of the grain centroids.
-3. Use [Phon](http://dx.doi.org/10.5281/zenodo.16711) to insert zero-thickness cohesive elements at the grain boundaries.
+3. Use [Phon](http://github.com/KristofferC/Phon) to insert zero-thickness cohesive elements at the grain boundaries.
 4. Parse both files in PhaceXL_2D to generate finite-thickness interface elements at the grain boundaries.
 5. Run finite-element analysis.
 
@@ -16,7 +16,7 @@ To generate a 2D polycrystal microstructure mesh with finite-thickness interface
 
 Run [PhaceXL_2D](main/phaceXL_2D.py) using   '*phaceXL_2D.py -i INPUTFILENAME -c CENTROIDFILENAME -p MPER -s SFACTOR*'   from the command line. The following input parameters have to be specified:
 
- *INPUTFILENAME*:   name of the file containing the polycrystal mesh and the cohesive elements generated by [Neper](http://neper.sourceforge.net/) and [Phon](http://dx.doi.org/10.5281/zenodo.16711),  
+ *INPUTFILENAME*:   name of the file containing the polycrystal mesh and the cohesive elements generated by [Neper](http://neper.sourceforge.net/) and [Phon](http://github.com/KristofferC/Phon),  
  *CENTROIDFILENAME*:   name of the file containing the coordinates of the grain centroids,  
  *MPER*:               indication whether the polycrystal microstructure mesh is periodic ('y') or (by default) non-periodic ('n'),  
  *SFACTOR*:            factor controlling the amount of shrinkage for each grain.