diff --git a/vignettes/v00_teaching_hydrology.rmd b/vignettes/v00_teaching_hydrology.rmd
index 5137a5a2df9ef66b7fc11d24e5b660f10b77da75..9efe01fec977f63201ae7befb0b96d94c27d521e 100644
--- a/vignettes/v00_teaching_hydrology.rmd
+++ b/vignettes/v00_teaching_hydrology.rmd
@@ -219,7 +219,7 @@ Note that if `WupPer` is not set by the user, by default `WupPer = NULL`, and in
 ```{r, v00_fig_wup, echo=TRUE, eval=TRUE, fig.width=3*3, fig.height=3*1.7, out.width='98%'}
 # Warm-up and simulation periods
 per_wup1m <- c("2002-12-01", "2002-12-31")
-per_wup4y <- c("1999-01-01", "2002-12-31")
+per_wup1y <- c("2002-01-01", "2002-12-31")
 per_sim   <- c("2003-01-01", "2006-12-31")
 
 # Parameter set
@@ -238,9 +238,9 @@ sim_wup1m <- SimGR(PrepGR  = prep_no_q,
                    SimPer  = per_sim)
 
 # Simulation with a 4-year warm-up period
-sim_wup4y <- SimGR(PrepGR  = prep_no_q, 
+sim_wup1y <- SimGR(PrepGR  = prep_no_q, 
                    Param   = param_gr4j,
-                   WupPer  = per_wup4y, 
+                   WupPer  = per_wup1y, 
                    SimPer  = per_sim)
 
 # Graphical comparison
@@ -248,12 +248,12 @@ col_wup <- c("orchid", "orange2", "green3")
 matplot(x = as.POSIXct(sim_wup0d$OutputsModel$DatesR),
         y = cbind(sim_wup0d$OutputsModel$Qsim,
                   sim_wup1m$OutputsModel$Qsim,
-                  sim_wup4y$OutputsModel$Qsim),
-        xlab = "time [days]", ylab = "flow [mm/day]", 
+                  sim_wup1y$OutputsModel$Qsim),
+        xlab = "time [d]", ylab = "flow [mm/d]", 
         type = "l", lty = 1, lwd = 2, col = col_wup,
         xlim = as.POSIXct(x = c("2003-01-01", "2003-09-01"), tz = "UTC"))
 legend("topright", 
-       legend = c("no warm-up", "1-month warm-up", "4-year warm-up"), 
+       legend = c("no warm-up", "1-month warm-up", "1-year warm-up"), 
        col = col_wup, lwd = 2,
        bg = "white")
 ```