diff --git a/r-cormas-cowat-example.R b/r-cormas-cowat-example.R
index d50134ca4af3f708ad93ab2bc98b0e0a43756dc3..86fa79d67ba0e15c56d95f3afbea847d4824fb28 100644
--- a/r-cormas-cowat-example.R
+++ b/r-cormas-cowat-example.R
@@ -73,5 +73,8 @@ for (day in 1:31){
 
 riverFlow <- getNumericProbe("riverFlow","COWAT")
 naturalRiverFlow <- getNumericProbe("naturalRiverFlow","COWAT")
-plot(0:240, riverLevel)
+plot(0:744, riverFlow, type ="l")
+lines(0:744, naturalRiverFlow, type ="l", col="red", lwd=2)
+points(days * 24, riverFlows, col="blue")
+ plot(river)