diff --git a/inst/doc/paper-figure-PlotForecastPDF.R b/inst/doc/paper-figure-PlotForecastPDF.R
index b8867a61c882f4cb6453a64226c0539aae74aa8a..bebc49725a02dec4c256a398974c532705d3833a 100644
--- a/inst/doc/paper-figure-PlotForecastPDF.R
+++ b/inst/doc/paper-figure-PlotForecastPDF.R
@@ -34,9 +34,9 @@ c(hcst, hcst_ref) %<-% CST_Load(var = 'prlr',
                                 leadtimemin = 1, leadtimemax = 214, 
                                 nmember = 25, output = "lonlat") 
 hcst$data <- hcst$data * 3600 * 24 * 1000
-attributes(hcst$attrs$Variable)$units <- 'mm'
+hcst$attrs$Variable$metadata$prlr$units <- 'mm'
 hcst_ref$data <- hcst_ref$data * 3600 * 24 * 1000
-attributes(hcst_ref$attrs$Variable)$units <- 'mm'
+hcst_ref$attrs$Variable$metadata$prlr$units <- 'mm'
 
 
 c(fcst, obs) %<-% CST_Load(var = 'prlr',
@@ -49,9 +49,9 @@ c(fcst, obs) %<-% CST_Load(var = 'prlr',
                            leadtimemin = 1, leadtimemax = 214,
                            nmember = 50, output = "lonlat")
 fcst$data <- fcst$data * 1000 * 3600 * 24
-attributes(fcst$attrs$Variable)$units <- 'mm'
+fcst$attrs$Variable$metadata$prlr$units <- 'mm'
 obs$data <- obs$data  * 1000 * 3600 * 24
-attributes(obs$attrs$Variable)$units <- 'mm'
+obs$attrs$Variable$metadata$prlr$units <- 'mm'
 
 
 fcst_QM <- CST_QuantileMapping(exp = hcst,