CST_QuantileMapping() not preserving the "datetime" dimension

Hi @nperez

 I am using CST_QuantileMapping() on several start dates and I am saving the outputs to files with CST_SaveExp(), one file for each start date.  

The output files seem to have an issue in the "date-time" dimension: If you explore the file i.e. with cdo sinfov, you will see that the date-time information is wrong, with any file having the same start date. The date in the filename is correct, however.

I reported below an example to reproduce the data with the date-time issue.

Thanks a lot! Best, Silvia

exp <- list(name = 'experiment', path = infile)
obs <- list(name = 'observation',path = reffile)

exp %<-% CST_Load('prlr', exp=list(exp), obs=NULL, sdates =startDates, nmember = nens, storefreq = "daily", sampleperiod = 1, latmin =latn, latmax = latx, lonmin = lonn, lonmax = lonx, output = 'lonlat',nprocs=6)
obs %<-% CST_Load('prlr', exp=list(obs), obs=NULL, sdates =startDates, nmember = 1, storefreq = "daily", sampleperiod = 1, latmin =latn, latmax = latx, lonmin = lonn, lonmax = lonx, output = 'lonlat',nprocs=6)
dim(exp$data)
dim(obs$data)

# SplitDim 
print("Splitting ftime by month ...")
expm   <- CST_SplitDim(exp, split_dim = c('ftime'))
obsm   <- CST_SplitDim(obs, split_dim = c('ftime'))
dim(expm$data)
dim(obsm$data)

exp.qm  <- CST_QuantileMapping(expm,obsm, method = "QUANT", wet.day=F,sample_dims = c('member', 'sdate', 'ftime'))
dim(exp.qm$data)

# MergeDims
print("Restoring original dimensions ...")
exp.qmm <- CST_MergeDims(exp.qm, merge_dims=c("ftime", "monthly"), na.rm=TRUE)

print (paste("Save Bias-Adjusted data to", outdir))
CST_SaveExp(exp.qmm, outpath)