diff --git a/R/Intlr.R b/R/Intlr.R
index beb8d7ee9ca6d8d11f96c025bf8408c42c3b8975..8b5dabea155f8db79f1ad59c30ad3fcf5ffe2b45 100644
--- a/R/Intlr.R
+++ b/R/Intlr.R
@@ -654,7 +654,7 @@ Intlr <- function(exp, obs, exp_cor = NULL, exp_lats, exp_lons, obs_lats, obs_lo
 
   # restore ensemble dimension in observations if it existed originally
   if (restore_ens) {
-    predictand <- MergeDims (predictand, merge_dims = c(time_dim, sdate_dim), rename_dim = sdate_dim)
+    predictand <- s2dv::InsertDim(predictand, posdim = 1, lendim = 1, name = member_dim)
   }
   
   # Return a list of three elements