diff --git a/Doxygenfile b/Doxygenfile
index dbddd16f..e55b2533 100644
--- a/Doxygenfile
+++ b/Doxygenfile
@@ -869,7 +869,7 @@ FILTER_SOURCE_PATTERNS =
 # (index.html). This can be useful if you have a project on for instance GitHub
 # and want to reuse the introduction page also for the doxygen output.
 
-USE_MDFILE_AS_MAINPAGE =
+USE_MDFILE_AS_MAINPAGE = mainpage.md
 
 #---------------------------------------------------------------------------
 # Configuration options related to source browsing
diff --git a/doc/BG_Flood.bib b/doc/BG_Flood.bib
new file mode 100644
index 00000000..afa7a5d7
--- /dev/null
+++ b/doc/BG_Flood.bib
@@ -0,0 +1,18 @@
+@article{Vacondio2017,
+author = {Vacondio, Renato and Dal Pal\`{u}, Alessandro and Ferrari, Alessia and Mignosa, Paolo and Aureli, Francesca and Dazzi, Susanna},
+title = {A Non-Uniform Efficient Grid Type for GPU-Parallel Shallow Water Equations Models},
+year = {2017},
+issue_date = {February 2017},
+publisher = {Elsevier Science Publishers B. V.},
+address = {NLD},
+volume = {88},
+number = {C},
+issn = {1364-8152},
+url = {https://doi.org/10.1016/j.envsoft.2016.11.012},
+doi = {10.1016/j.envsoft.2016.11.012},
+abstract = {A GPU-parallel numerical model for the solution of the 2D Shallow Water Equations, based on a novel type of grid called Block-Uniform Quadtree (BUQ), is presented. BUQ grids are based on a data structure which allows to exploit the computational capability of GPUs with minimum overheads, while discretizing the domain with non-uniform resolution. Different cases have been simulated in order to assess the efficiency of the BUQ grids. Theoretical and laboratory tests demonstrate that speed-ups of up to one order of magnitude can be achieved in comparison with uniform Cartesian grids. In the simulation of a hypothetical flood event induced by a levee breach in a real 83\'{z}km long river reach, with maximum resolution of 5\'{z}m, a ratio of physical to computational time of about 12 was obtained, opening scenarios of quasi real-time 2D simulations in large domains, still retaining a high resolution where necessary. Shallow Water Equations parallel numerical scheme suitable for GPU.Novel Block-Uniform Quadtree (BUQ) grid-type, exploiting the computational capability of GPUs with non-uniform resolution.Very large domain can be simulated (38\'{z}km river reach) with high resolution (2\'{z}m).Speed-up of one order of magnitude is achieved, in comparison with Uniform resolution Cartesian Grids.High ratio between physical and computational time which enables real-time simulations.},
+journal = {Environ. Model. Softw.},
+month = {feb},
+pages = {119–137},
+numpages = {19}
+}
diff --git a/doc/mainpage.md b/doc/mainpage.md
index 580f7591..af99212d 100644
--- a/doc/mainpage.md
+++ b/doc/mainpage.md
@@ -9,3 +9,8 @@ Test picture:
 ![ILCL-table](ILCL-values.png)
 
 
+Test biblio:
+paper from Vacondo: \cite{Vacondio2017}
+Biblography:
+\bibliographystyle{unsrt}
+\bibliography{reference}
\ No newline at end of file