diff --git a/R/utils.RunModel.R b/R/utils.RunModel.R
index 9f5bb678da6ecfcbdca1107a461ad83338d6e7a2..71b0512fbd7f586e2f0911881dd1793677e1d388 100644
--- a/R/utils.RunModel.R
+++ b/R/utils.RunModel.R
@@ -148,6 +148,9 @@ merge.OutputsModel <- function(x, y, ...) {
   for (item in items) {
     y[[item]] <- c(x[[item]], y[[item]])
   }
+  # We keep original warm-up data
+  if (!is.null(x$RunOptions$WarmUpQsim)) y$RunOptions$WarmUpQsim <- x$RunOptions$WarmUpQsim
+  if (!is.null(x$RunOptions$WarmUpQsim_m3)) y$RunOptions$WarmUpQsim_m3 <- x$RunOptions$WarmUpQsim_m3
   return(y)
 }