diff --git a/modules/Anomalies/Anomalies.R b/modules/Anomalies/Anomalies.R
index f8b49da8b4743e6c5505d5da9144e655210a0869..552f895a8ae3944c45aad0c96d169e7edffbce17 100644
--- a/modules/Anomalies/Anomalies.R
+++ b/modules/Anomalies/Anomalies.R
@@ -29,7 +29,8 @@ compute_anomalies <- function(recipe, data) {
 			memb_dim = 'ensemble',
 			dim_anom = 'syear',
 			dat_dim = c('dat', 'ensemble'),
-			ftime_dim = 'time')
+			ftime_dim = 'time',
+			ncores = recipe$Analysis$ncores)
     # Reorder dims
     anom$exp$data <- Reorder(anom$exp$data, names(original_dims))
     anom$obs$data <- Reorder(anom$obs$data, names(original_dims))
@@ -59,7 +60,8 @@ compute_anomalies <- function(recipe, data) {
         	         dat_dim = c("dat", "ensemble"),
         	         memb = FALSE,
         	         memb_dim = "ensemble",
-        	         ftime_dim = "time")
+        	         ftime_dim = "time",
+			 ncores = recipe$Analysis$ncores)
       clim_hcst <- InsertDim(clim$clim_exp, posdim = 1, lendim = 1,
                              name = "syear")
       dims <- dim(clim_hcst)
@@ -95,9 +97,3 @@ compute_anomalies <- function(recipe, data) {
 	      hcst.full_val = hcst_fullvalue, obs.full_val = obs_fullvalue))
 
 }
-
-
-
-
-
-
diff --git a/modules/Anomalies/tmp/CST_Anomaly.R b/modules/Anomalies/tmp/CST_Anomaly.R
index a84b6fc8538b03f4113b96aa8b3126189a0bdee9..f38e39b050f7c46be452ac6e6571542c465264b9 100644
--- a/modules/Anomalies/tmp/CST_Anomaly.R
+++ b/modules/Anomalies/tmp/CST_Anomaly.R
@@ -53,12 +53,7 @@
 #'attr(exp, 'class') <- 's2dv_cube'
 #'attr(obs, 'class') <- 's2dv_cube'
 #'
-#'anom1 <- CST_Anomaly(exp = exp, obs = obs, cross = FALSE, memb = TRUE)
-#'anom2 <- CST_Anomaly(exp = exp, obs = obs, cross = TRUE, memb = TRUE)
-#'anom3 <- CST_Anomaly(exp = exp, obs = obs, cross = TRUE, memb = FALSE)
-#'anom4 <- CST_Anomaly(exp = exp, obs = obs, cross = FALSE, memb = FALSE)
-#'anom5 <- CST_Anomaly(lonlat_temp$exp)
-#'anom6 <- CST_Anomaly(obs = lonlat_temp$obs)
+#'anom <- CST_Anomaly(exp = exp, obs = obs, cross = FALSE, memb = TRUE)
 #'
 #'@seealso \code{\link[s2dv]{Ano_CrossValid}}, \code{\link[s2dv]{Clim}} and \code{\link{CST_Load}}
 #'
@@ -178,11 +173,21 @@ CST_Anomaly <- function(exp = NULL, obs = NULL, dim_anom = 'sdate', cross = FALS
   
   # With cross-validation
   if (cross) {
-    ano <- Ano_CrossValid(exp = exp$data, obs = obs$data, time_dim = dim_anom, memb_dim = memb_dim, memb = memb, dat_dim = dat_dim)
+    ano <- Ano_CrossValid(exp = exp$data, obs = obs$data,
+			  time_dim = dim_anom,
+			  memb_dim = memb_dim,
+			  memb = memb,
+			  dat_dim = dat_dim,
+			  ncores = ncores)
    
-    #  Without cross-validation 
+  # Without cross-validation 
   } else {
-    tmp <- Clim(exp = exp$data, obs = obs$data, time_dim = dim_anom, memb_dim = memb_dim, memb = memb, dat_dim = dat_dim)
+    tmp <- Clim(exp = exp$data, obs = obs$data,
+		time_dim = dim_anom,
+		memb_dim = memb_dim,
+		memb = memb,
+		dat_dim = dat_dim,
+		ncores = ncores)
     if (!is.null(filter_span)) {
       tmp$clim_exp <- Apply(tmp$clim_exp, 
                             target_dims = c(ftime_dim),
diff --git a/modules/Loading/testing_recipes/recipe_seasonal-tests.yml b/modules/Loading/testing_recipes/recipe_seasonal-tests.yml
index e1857ac06d77f649c29dff13756a89b02a27c81c..61177b7162e55b07e8471e59a2f4a493e7cfb42e 100644
--- a/modules/Loading/testing_recipes/recipe_seasonal-tests.yml
+++ b/modules/Loading/testing_recipes/recipe_seasonal-tests.yml
@@ -30,6 +30,9 @@ Analysis:
   Workflow:
     Calibration:
       method: mse_min
+    Anomalies:
+      compute: yes
+      cross_validation: yes
     Skill:
       metric: RPS RPSS CRPS CRPSS FRPSS BSS10 BSS90 EnsCorr Corr mean_bias mean_bias_SS
     Probabilities: