diff --git a/DESCRIPTION b/DESCRIPTION
index 4a54654b964cc7ca487ce948214b9759aed2dadd..a7bdd2bbee6cfae86e2723387005174816e77ed7 100644
--- a/DESCRIPTION
+++ b/DESCRIPTION
@@ -11,7 +11,7 @@ Imports:
     abind,
     bigmemory,
     future,
-    multiApply,
+    multiApply (>= 2.0.0),
     parallel
 Suggests:
     easyNCDF,
diff --git a/R/ByChunks.R b/R/ByChunks.R
index 6b4fd07e598829e189ddf754465c510aa76cf57c..c43f4d14c28fba3b4cc0eda1172e191a92d2cf8b 100644
--- a/R/ByChunks.R
+++ b/R/ByChunks.R
@@ -624,7 +624,7 @@ ByChunks <- function(step_fun, cube_headers, ..., chunks = 'auto',
       t_begin_compute <- Sys.time()
       result <- multiApply::Apply(data, 
                                   target_dims = attr(step_fun, 'TargetDims'), 
-                                  AtomicFun = step_fun, ..., 
+                                  fun = step_fun, ..., 
                                   output_dims = attr(step_fun, 'OutputDims'),
                                   ncores = threads_compute)
       if (!found_first_result) {
diff --git a/inst/chunking/load_process_save_chunk.R b/inst/chunking/load_process_save_chunk.R
index 5606b0019b25f710b1d4c4cdcd3650f9034b576f..c552aa2e2b29bbd1893e01a496cf7f3fabcd0283 100644
--- a/inst/chunking/load_process_save_chunk.R
+++ b/inst/chunking/load_process_save_chunk.R
@@ -81,7 +81,7 @@ res <- do.call("Apply",
                c(
                  list(data,
                       target_dims = attr(fun, 'TargetDims'),
-                      AtomicFun = fun,
+                      fun = fun,
                       output_dims = attr(fun, 'OutputDims'), 
                       ncores = threads_compute),
                  params