diff --git a/README.md b/README.md
index 2cde4067a7f0deca6c92997386d8191ab7330f37..0f44675c3d9e0c30c396f3e8d4b5f4a60d229179 100644
--- a/README.md
+++ b/README.md
@@ -11,59 +11,64 @@ How to run it
 
 // Comment : For some collision operators (Sugama and Full Coulomb) you have to run COSOlver from B.J.Frei first in order to generate the required matrices in HeLaZ/iCa folder.
 
-# Road map
-(Current version : 2.5)
+# Changelog
 
-0. Write MOLI matlab solver in Fortran using Monli1D as starting point
+2. MPI parallel version
 
-	0.0 go from 1D space to 2D fourier and from Hermite basis to Hermite-Laguerre basis
+	>2.5 GK Sugama collision operator
 
-	0.1 implement linear Poisson equation in fourier space
+	2.4 2D cartesian parallel (along p and kr)
 
-	0.2 implement moment hierarchy linear terms
+	2.3 GK Dougherty operator
 
-	0.3 RK4 time solver
+	2.2 Allow restart with different P,J values (results are not concluents)
 
-	0.4 Benchmark with MOLI matlab results for Z-pinch (cf. kz_linear script)
+	2.1 First compilable parallel version (1D parallel along kr)
 
-	0.5 Load COSOlver matrices
+1. Implementation of the non linear Poisson brackets term
 
-	0.6 Benchmarks now include Dougherty, Lenard-Bernstein and Full Coulomb collision operators
+	1.4 Quantitative study with stationary average particle flux \Gamma_\infty
 
-1. Implementation of the non linear Poisson brackets term
+	1.3 Linear analysis showed that a certain amount of PJ are recquired to trigger mode
 
-	1.0 FFTW3 has been used to treat the convolution as a product and discrete fourier transform
+	1.2 Zonal flows are observed in a similar way to Ricci Rogers 2006 with GS2
+
+	1.1 Qualitative test : find similar turbulences as Hasegawa Wakatani system with few moments
 
 	1.1 Methods in fourier_mod.f90 have been validated by tests on Hasegawa Wakatani system
 
-	1.1 Qualitative test : find similar turbulences as Hasegawa Wakatani system with few moments
+	1.1 Methods in fourier_mod.f90 have been validated by tests on Hasegawa Wakatani system
 
-	1.2 Zonal flows are observed in a similar way to Ricci Rogers 2006 with GS2
+	1.0 FFTW3 has been used to treat the convolution as a product and discrete fourier transform
 
-	1.3 Linear analysis showed that a certain amount of PJ are recquired to trigger mode
+0. Write MOLI matlab solver in Fortran using Monli1D as starting point
 
-	1.4 Quantitative study with stationary average particle flux \Gamma_\infty
+	0.6 Benchmarks now include Dougherty, Lenard-Bernstein and Full Coulomb collision operators
 
-2. MPI parallel version
+	0.5 Load COSOlver matrices
 
-	2.1 First compilable parallel version (1D parallel along kr)
+	0.4 Benchmark with MOLI matlab results for Z-pinch (cf. kz_linear script)
 
-	2.2 Allow restart with different P,J values (results are not concluents)
+	0.3 RK4 time solver
 
-	2.3 GK Dougherty operator
+	0.2 implement moment hierarchy linear terms
 
-	2.4 2D cartesian parallel (along p and kr)
+	0.1 implement linear Poisson equation in fourier space
+
+	0.0 go from 1D space to 2D fourier and from Hermite basis to Hermite-Laguerre basis
 
->2.5 GK Sugama collision operator
+# Roadmap
+
+2. MPI parallel version
 
 	2.6 GPU accelerated version
 
 	2.7 GK Full Coulomb collision operator
 
-4. GK 3D version, kr,kz,kpar for linear device
+3. GK 3D version, kr,kz,kpar for linear device
 
-5. DK 3D version, kr,kz,kpar for linear device
+4. DK 3D version, kr,kz,kpar for linear device
 
-6. DK+GK 3D version, kr,kz,kpar for linear device
+5. DK+GK 3D version, kr,kz,kpar for linear device
 
-7. 3D version with curvature
+6. 3D version with curvature