diff --git a/README.md b/README.md
index 4562803..b863a83 100644
--- a/README.md
+++ b/README.md
@@ -7,7 +7,8 @@ ULiege, 2018-2020
 ## Features and limitations
 
 ### What is GeoGen?
-GeoGen is a python code that generates an aerodynamic geometry which can then be meshed in gmsh and used in CFD solvers.
+GeoGen is a python code that generates an aerodynamic geometry which can then be meshed using gmsh and used in CFD solvers.  
+GeoGen v2 and up is compatible with python 3, while previous versions are compatible with python 2.
 
 ### What can GeoGen do?
 GeoGen currently supports the following configurations:
@@ -22,49 +23,43 @@ GeoGen currently supports the following configurations:
 
 GeoGen was primarly designed to be used with [SU2](https://github.com/su2code/SU2) and [waves](https://gitlab.uliege.be/am-dept/waves), but any solver interfaced with [gmsh](http://gmsh.info/) can be used!
 
-### Cite us!
-If you use this work, please acknowledge the authors:  
-```text 
-"GeoGen - a python/gmsh automated geometry generator for CFD solvers" by Adrien Crovato, University of Liege
-```  
-
 ## Usage
 The script is run through the command line:
 ```sh
 python geoGen.py path/to/config/file.py <-o path/to/output/file.geo>
 ```
-If no output file is provided, a workspace directory will be created and the geometry will be stored inside as grid.geo.
+If no output file is provided, a workspace directory will be created and the geometry will be stored inside as `grid.geo`.
 
 The geometry is generated from a python file containing a dictionary of parameters. Examples are given in [config](config/) and the main options are summurized hereunder.
 
 **Parameters**
 
 Wing definition:
- - airfPath: relative path to the airfoils directory
- - airfName: array (size: nP+1) of names of file containing airfoil (Selig formatted) coordinates
- - span: array (size: nP) of span for each planform of the wing
- - taper: array (size: nP)) of taper of each planform of the wing
- - sweep: array (size: nP) of leading edge sweep of each planform of the wing 
- - dihedral: array (size: nP) of dihedral angle of each planform of the wing
- - twist: array (size: nP+1) of twist angle of each airfoil of the wing
- - rootChord: root chord (scalar) of the wing
- - offset: array of x and z offset (size: 2) applied to the leading edge of the root section
- - coWingtip: boolean, True for cutoof wingtip, Fasle for rounded wingtip (not supported yet)
+ - `airfPath`: relative path to the airfoils directory
+ - `airfName`: array (size: nP+1) of names of file containing airfoil (Selig formatted) coordinates
+ - `span`: array (size: nP) of span for each planform of the wing
+ - `taper`: array (size: nP)) of taper of each planform of the wing
+ - `sweep`: array (size: nP) of leading edge sweep of each planform of the wing 
+ - `dihedral`: array (size: nP) of dihedral angle of each planform of the wing
+ - `twist`: array (size: nP+1) of twist angle of each airfoil of the wing
+ - `rootChord`: root chord (scalar) of the wing
+ - `offset`: array of x and z offset (size: 2) applied to the leading edge of the root section
+ - `coWingtip`: boolean, True for cutoof wingtip, Fasle for rounded wingtip (not supported yet)
 
 
 Domain definition:
- - domType: string, box for box-shaped domain or shpere for shperical-shaped domain
+ - `domType`: string, box for box-shaped domain or shpere for shperical-shaped domain
 
 
 if domain type is shpere, typical for Euler equations:
- - rSphere: radius of the sphere (scalar)
+ - `rSphere`: radius of the sphere (scalar)
 
 
 if domain type is a box (a wake will then be defined), typically for potential equations:
- - xoBox: x-coordinate (scalar) of the origin of the box 
- - xfBox: x-coordinate (scalar) of the end of the box
- - yfBox: y-coordinate (scalar) of the end of the box
- - zoBox: z-coordinate (scalar) of the origin of the box
- - zfBox: z-coordinate (scalar) of the end of the box
- - nSlope: number (scalar) of airfoil geometrical points counted from TE used to compute wake slope
+ - `xoBox`: x-coordinate (scalar) of the origin of the box 
+ - `xfBox`: x-coordinate (scalar) of the end of the box
+ - `yfBox`: y-coordinate (scalar) of the end of the box
+ - `zoBox`: z-coordinate (scalar) of the origin of the box
+ - `zfBox`: z-coordinate (scalar) of the end of the box
+ - `nSlope`: number (scalar) of airfoil geometrical points counted from TE used to compute wake slope
 
diff --git a/domain.py b/domain.py
index ae96c0c..b4fa243 100644
--- a/domain.py
+++ b/domain.py
@@ -1,4 +1,4 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 # -*- coding: utf-8 -*-
 
 ''' 
diff --git a/geoGen.py b/geoGen.py
index 2689636..bbe9751 100755
--- a/geoGen.py
+++ b/geoGen.py
@@ -1,4 +1,4 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 # -*- coding: utf-8 -*-
 
 ''' 
diff --git a/tip.py b/tip.py
index 6b5bb46..af3c1fb 100644
--- a/tip.py
+++ b/tip.py
@@ -1,4 +1,4 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 # -*- coding: utf-8 -*-
 
 ''' 
diff --git a/wake.py b/wake.py
index 7b3abaf..d90ec1c 100644
--- a/wake.py
+++ b/wake.py
@@ -1,4 +1,4 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 # -*- coding: utf-8 -*-
 
 ''' 
diff --git a/wing.py b/wing.py
index 3aacdf7..431642b 100644
--- a/wing.py
+++ b/wing.py
@@ -1,4 +1,4 @@
-#!/usr/bin/env python
+#!/usr/bin/env python3
 # -*- coding: utf-8 -*-
 
 '''