diff --git a/README.md b/README.md
index d7c19882be4c03dc3729d5162b703c6df05293f6..bf1a8aed8f712f168f1260ac6320adb5fabe2638 100644
--- a/README.md
+++ b/README.md
@@ -23,6 +23,22 @@ marker function `domain_marker_meshrez*.pvd` or `interface_marker_global_meshrez
 ### Curent starte
 Currently the code can do the following:
 
+- run a flexible combination of models, i.e. each subdomain can assume Richards/Richardson's equation or the full two-phase flow equations. 
+
+- run examples with exact solutions and compare the numerical solution to it. 
+
+- run mesh studies against for these exact solution examples. 
+
+### TODO
+
+- Parallelisation: Unfortunately atm the code does not run in parallel, although the different mesh sizes during mesh study calculations do run in parallel. The code is quite capsuled already, but the method `domainPatch.calc_gli_term()` in `./LDDsimulation/domainPatch.py` needs to be adapted if one were to parallelise the code. Mainly, one needs to figure out how to share pointers to the storages `domainPatch.pressure[phase]` and 
+`domainPatch.pressure_prev_timestep[phase]` between processes. Hopefully this is possible without copying the data that is used internally in `Fenics`
+into `numpy arrays` each time.
+
+- more modells would be nice, maybe more effects. 
+
+- more examples without exact solutions. Some Code needs to be written to compare a numerical solution to one with a finer meshrez.
+
 
 ## License
 `LDDsimulation` is licensed under the terms and conditions of the GNU General