CDORemap: it should automatically put 'time' dimension at the end of each chunk sent to cdo.

The following crashes, because the sent chunks have dims (time, lon, lat):

library(startR)

repos_path <- '/esarchive/exp/ecearth/i00k/monthly_mean/heatc/ohc_2d_avg_0-300m_i00k_$sdate$_fc0-4_*.nc'
heatc <- Start(repos = repos_path,
               sdate = 'first', 
               var = 'heatc_sl', 
               time = 'all',
               ensemble = 'all', 
               x = 'all', 
               y = 'all', 
               return_vars = list(var_names = NULL), 
               var_var = 'var_names')

mask_path <- '/esnas/autosubmit/con_files/mask.regions.Ec3.0_O1L46.nc'
lon <- Start(repos = mask_path,
             var = 'nav_lon', 
             t = 'first',
             z = 'first', 
             x = 'all', 
             y = 'all', 
             return_vars = list(var_names = NULL), 
             var_var = 'var_names')
lat <- Start(repos = mask_path,
             var = 'nav_lat', 
             t = 'first',
             z = 'first', 
             x = 'all', 
             y = 'all', 
             return_vars = list(var_names = NULL), 
             var_var = 'var_names')


library(s2dverification)
heatc2 <- heatc
lon2 <- lon
lat2 <- lat
attr(heatc2$Data, 'variables') <- NULL
names(dim(heatc2$Data))[4] <- 't'
attr(lon2$Data, 'variables') <- NULL
attr(lat2$Data, 'variables') <- NULL
heatc3 <- CDORemap(Subset(heatc2$Data, 1:3, list(1, 1, 1), drop = 'selected'),
                   Subset(lon2$Data, 1:4, list(1, 1, 1, 1), drop = 'selected'),
                   Subset(lat2$Data, 1:4, list(1, 1, 1, 1), drop = 'selected'), 
                   't106grid', 'bil')

Error in R_nc4_def_var_float: NetCDF: NC_UNLIMITED in the wrong index
Name of variable that the error occurred on: "var"
[1] "----------------------"
[1] "Var:  var"
[1] "Ndims:  3"
[1] "Dimids: "
[1] 2 1 0
Error in ncvar_add(nc, vars[[ivar]], verbose = verbose, indefine = TRUE) : 
  Error in ncvar_add, defining var var