Update README.md

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