#############################################################################
# Custom function to bind parallelization output.
#############################################################################
comb_rbind <- function(x, ...) {  
  mapply(rbind,x,...,SIMPLIFY=FALSE)
}