diff --git a/README.md b/README.md
index 386bc1a8c169b03920aafbff890bd4e6973aa36f..88bfcd0889dd149a5dce237bf36f289a7e0c6590 100644
--- a/README.md
+++ b/README.md
@@ -28,77 +28,44 @@ along with this program.  If not, see <https://www.gnu.org/licenses/>.
 # Changelog
 
 4. GYACOMO
-
   4.1 Miller geometry is added and benchmarked for CBC adiabatic electrons
-
   4.0 new naming and opening the code with GNU GPLv3 license
 
 3. HeLaZ 3D
-
   3.9 HeLaZ can now evolve electromagnetic fluctuations by solving Ampere equations (benchmarked linearly)
-
 	3.8 HeLaZ has been benchmarked for CBC with GENE for various gradients values (see Dimits_fig3.m)
-
 	3.7 The frequency plane has been transposed from positive kx to positive ky for easier implementation of shear. Also added 3D zpinch geometry
-
 	3.6 HeLaZ is now parallelized in p, kx and z and benchmarked for each parallel options with gbms (new molix) for linear fluxtube shearless.
-
 	3.5 Staggered grid for parallel odd/even coupling
-
 	3.4 HeLaZ can run with adiabatic electrons now!
-
 	3.3 HeLaZ 3D has been benchmarked in fluxtube salphaB geometry linear run with molix (B.J.Frei) code and works now for shear = 0 with periodic z BC
-
 	3.2 Stopping file procedure like in GBS is added
-
 	3.1 Implementation of mirror force
-
-	3.0 HeLaZ is now 3D and works like HeLaZ 2D if Nz = 1, the axis were renamed (r,z) -> (x,y,z) and now the parallel direction is ez. All arrays have been extended, diagnostics and analysis too. The linear coefficients are now precomputed with lin_coeff_and_geometry routines.
+  3.0 HeLaZ is now 3D and works like HeLaZ 2D if Nz = 1, the axis were renamed (r,z) -> (x,y,z) and now the parallel direction is ez. All arrays have been extended, diagnostics and analysis too. The linear coefficients are now precomputed with lin_coeff_and_geometry routines.
 
 2. MPI parallel version
-
 	2.7 Versatile interpolation of kperp for the cosolver matrices and corrections done on DGGK
-
 	2.6 Change of collisionality normalisation (from nu_ei to nu_ii), implementation of FCGK
-
 	2.5 GK cosolver collision implementation
-
 	2.4 2D cartesian parallel (along p and kr)
-
 	2.3 GK Dougherty operator
-
 	2.2 Allow restart with different P,J values (results are not concluents)
-
 	2.1 First compilable parallel version (1D parallel along kr)
 
 1. Implementation of the non linear Poisson brackets term
-
 	1.4 Quantitative study with stationary average particle flux \Gamma_\infty
-
 	1.3 Linear analysis showed that a certain amount of PJ are recquired to trigger mode
-
 	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 Methods in fourier_mod.f90 have been validated by tests on Hasegawa Wakatani system
-
 	1.0 FFTW3 has been used to treat the convolution as a product and discrete fourier transform
 
 0. Write MOLI matlab solver in Fortran using Monli1D as starting point
-
 	0.6 Benchmarks now include Dougherty, Lenard-Bernstein and Full Coulomb collision operators
-
 	0.5 Load COSOlver matrices
-
 	0.4 Benchmark with MOLI matlab results for Z-pinch (cf. kz_linear script)
-
 	0.3 RK4 time solver
-
 	0.2 implement moment hierarchy linear terms
-
 	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