diff --git a/nature/fig/U.pdf b/nature/fig/U.pdf
index 132a588..33ba638 100644
Binary files a/nature/fig/U.pdf and b/nature/fig/U.pdf differ
diff --git a/nature/fig/out_orbitals.pdf b/nature/fig/out_orbitals.pdf
index 2ef9811..8c32d65 100644
Binary files a/nature/fig/out_orbitals.pdf and b/nature/fig/out_orbitals.pdf differ
diff --git a/nature/fig/ts.pdf b/nature/fig/ts.pdf
index 953f5ff..2334f1e 100644
Binary files a/nature/fig/ts.pdf and b/nature/fig/ts.pdf differ
diff --git a/nature/library.bib b/nature/library.bib
index ee1285c..d44e247 100644
--- a/nature/library.bib
+++ b/nature/library.bib
@@ -84,15 +84,10 @@
 	journal={PhD. Universidade de Cantabria},
 	year={2015}
 }
-@article{poncet2022,
-	title = {In-situ measurements of energetic depth-limited wave loading},
-	journal = {Applied Ocean Research},
-	volume = {125},
-	pages = {103216},
-	year = {2022},
-	issn = {0141-1187},
-	doi = {https://doi.org/10.1016/j.apor.2022.103216},
-	url = {https://www.sciencedirect.com/science/article/pii/S0141118722001572},
-	author = {P.A. Poncet and B. Liquet and B. Larroque and D. D’Amico and D. Sous and S. Abadie},
-	keywords = {Wave impact, Breaking wave, Loading, Breakwater, Field measurement, Pressure impulse, Multiple linear regression, Wind, Water level},
+@phdthesis{poncet2021,
+	title={Characterization of wave impact loading on structures at full scale: field experiment, statistical analysis and 3D advanced numerical modeling},
+	author={Poncet, Pierre-Antoine},
+	year={2021},
+	school={Université de Pau et des Pays de l'Adour},
+	chapter={4},
 }
diff --git a/nature/main.tex b/nature/main.tex
index 3b2746b..de39e93 100644
--- a/nature/main.tex
+++ b/nature/main.tex
@@ -65,7 +65,7 @@ for the study of multiphase incompressible flows.
 In this paper, we first use a one-dimensionnal depth-averaged non-linear non-hydrostatic model to verify that the
 signal measured by the wave buoy can be used as an incident wave input for the determination of hydrodynamic conditions
 near the breakwater. For this model, we use a SWASH model \parencite{zijlema2011} already calibrated by
-\textcite{poncet2022} on a domain reaching 1450m offshore of the breakwater.
+\textcite{poncet2021} on a domain reaching 1450m offshore of the breakwater.
 
 Then, we use a nested VOF model in two vertical dimensions that uses the output from the larger scale SWASH model as
 initial and boundary conditions to obtain the hydrodynamic conditions on the breakwater. The models uses olaFlow
@@ -129,12 +129,16 @@ the crest increases, with a zone reaching 400m long in front of the wave where t
 \subsection{Hydrodynamic conditions on the breakwater}
 
 The two-dimensionnal olaFlow model near the breakwater allowed to compute the flow velocity near and on the breakwater
-during the passage of the identified wave.
+during the passage of the identified wave. The results displayed in Figure~\ref{fig:U} show that the flow velocity
+reaches a maximum of 14.5m/s towards the breakwater during the identified extreme wave. The maximum reached velocity is
+similar to earlier shorter waves (at t=100s and t=120s), but the flow velocity remains high for twice as long as during
+those earlier waves. The tail of the identified wave also exhibits a water level over 5m for over 40s.
 
 \begin{figure*}
 	\centering
 	\includegraphics{fig/U.pdf}
-	\caption{Horizontal velocity computed with the olaFlow model at x=-20m on the breakwater armor}\label{fig:U}
+	\caption{Horizontal flow velocity computed with the olaFlow model at x=-20m on the breakwater armor. The identified
+	wave reaches this point around t=175s.}\label{fig:U}
 \end{figure*}
 
 \section{Discussion}