diff --git a/inst/doc/usecase/ex1_12_rotated_coordinates.R b/inst/doc/usecase/ex1_12_rotated_coordinates.R
index 3da4b794a376c71de0e6848e927b78e7824cf86e..f657a70418aead4f1bc7770e430f4a887f37b03a 100644
--- a/inst/doc/usecase/ex1_12_rotated_coordinates.R
+++ b/inst/doc/usecase/ex1_12_rotated_coordinates.R
@@ -5,11 +5,11 @@
 #-------------------------------
 library(startR)
 
-obs_path <- '/esarchive/exp/monarch-dustclim/3hourly/$var$/$var$_$date$03_av_an.nc'
+path <- '/esarchive/oper/thredds-dust/monarch-dustclim/3hourly/$var$-av_an/$var$_$date$03_av_an.nc'
 
 date <- c('20131229', '20131230')
 # two temporal dimensions: one for days and another four hours
-data_split <- Start(dataset = obs_path,
+data_split <- Start(dataset = path,
               var = 'od550du',
               lev = 'all',
               date = date,
@@ -21,7 +21,7 @@ data_split <- Start(dataset = obs_path,
                                  lat = NULL, lon = NULL),
               retrieve = TRUE, num_procs = 1)
 # one temporal dimension
-data <- Start(dataset = obs_path,
+data <- Start(dataset = path,
               var = 'od550du',
               lev = 'all',
               date = date,
@@ -48,7 +48,7 @@ attributes(data)$Variables$common$time
 
 
 # PLOT:
-map <- data[1,1,1,1,1,,] # Time-step 1
+map <- data[1,1,1,1,,] # Time-step 1
 library(sf)
 library(ggplot2)
 map <- data.frame(rlat = rep(rlat, length(rlon)),
diff --git a/inst/doc/usecase/ex2_1_timedim.R b/inst/doc/usecase/ex2_1_timedim.R
index f06b766dd1fff0db1ef3a02b0c1005ad376e2752..f8253c85162db28d4471e56c47d1da855de42861 100644
--- a/inst/doc/usecase/ex2_1_timedim.R
+++ b/inst/doc/usecase/ex2_1_timedim.R
@@ -13,16 +13,14 @@ library(startR)
                 return_vars = list(latitude = 'dat', longitude = 'dat', time = 'sdate'),
                 retrieve = FALSE)
   
-  library(multiApply)
   fun_spring <- function(x) {
-      source("/esarchive/scratch/nperez/Season_v2.R")
-      y <- Season_v2(x,  monini = 1, moninf = 3, monsup = 5)
+      y <- s2dv::Season(x,  time_dim = 'time', monini = 1, moninf = 3, monsup = 5)
       return(y)
   }
   
   step1 <- Step(fun = fun_spring,
-               target_dims = c('time'),
-               output_dims = c('time'))
+               target_dims = c('var', 'time'),
+               output_dims = c('var', 'time'))
   
   wf1 <- AddStep(data, step1)