diff --git a/vignettes/V03_First_Calibration.Rmd b/vignettes/V03_First_Calibration.Rmd
index 2c470e6f0e4291f4fe5a8e007702ac33231abe40..c0c4f909365fc798a9d4c7601936591a3bfbffbd 100644
--- a/vignettes/V03_First_Calibration.Rmd
+++ b/vignettes/V03_First_Calibration.Rmd
@@ -44,4 +44,34 @@ str(CalibOptions)
 OutputsCalib <- Calibration(InputsModel, RunOptions, InputsCrit, CalibOptions)
 ```
 
+```{r}
+save(OutputsCalib, file = "_cache/V03.RData")
+```
+
+
+## Run model with Michel calibration
+
+```{r}
+giropMichel <- girop
+for(id in giropMichel$id) {
+  giropMichel$params[giropMichel$id == id] <- list(OutputsCalib[[id]]$Param)
+}
+
+OutputsModels <- RunModel(
+  InputsModel = InputsModel,
+  RunOptions = RunOptions,
+  girop = giropMichel
+)
+```
+
+## Plot the result for each basin
+
+```{r, fig.height = 5, fig.width = 8}
+htmltools::tagList(lapply(
+  names(OutputsModels), 
+  function(x) {
+    plot(OutputsModels[[x]], Qobs = gits[[x]]$Qobs[RunOptions[[x]]$IndPeriod_Run] , main = x)
+  }
+))
+```
 
diff --git a/vignettes/v02_First_run.Rmd b/vignettes/v02_First_run.Rmd
index ae3050f458663e97f847d489e5e096e0b490db0a..67ca6f03cfbfd7795051b7cce3ca43853744da5a 100644
--- a/vignettes/v02_First_run.Rmd
+++ b/vignettes/v02_First_run.Rmd
@@ -104,7 +104,7 @@ RunOptions <- CreateRunOptions(
 ## Run the SD model for the whole basin
 
 ```{r}
-OutputsModels <- RunModel(
+OutputsModelsClimAware <- RunModel(
   InputsModel = InputsModel,
   RunOptions = RunOptions,
   girop
@@ -121,9 +121,9 @@ save(RunOptions, file = "_cache/V02.RData")
 
 ```{r, fig.height = 5, fig.width = 8}
 htmltools::tagList(lapply(
-  names(OutputsModels), 
+  names(OutputsModelsClimAware), 
   function(x) {
-    plot(OutputsModels[[x]], Qobs = gits[[x]]$Qobs[IndPeriod_Run] , main = x)
+    plot(OutputsModelsClimAware[[x]], Qobs = gits[[x]]$Qobs[IndPeriod_Run] , main = x)
   }
 ))
 ```