mcmcsave = function(rep,R,name_pattern="mcmc_", save_pattern="dMcmc_") {
	nmp			= name_pattern
	snp			= save_pattern
	mobjects 	= ls(pattern = nmp, envir=.GlobalEnv)

if(rep==1) {
	for(ii in 1:(length(mobjects))) {
		namei 	= substr(mobjects[ii] ,nchar(nmp)+1,nchar(mobjects[ii]))

		if(is.null(dim(get(mobjects[ii])))) {dims=length(get(mobjects[ii]))} else {
		dims	= dim(get(mobjects[ii]))}

		assign(paste(snp, namei, sep="")  , array(dim=c(R,dims)), envir = .GlobalEnv)
}}

for(ii in 1:(length(mobjects))) {
	namei 	= substr(mobjects[ii] ,nchar(nmp)+1,nchar(mobjects[ii]))
	dims	= dim(get(mobjects[ii]))

	if(length(dims)==1) 	{
		if(dims==1) {
		tmp = get((mobjects[ii]))
		eval(parse(text=paste(snp,namei,"[rep]","<-tmp",sep="") ))
		} else {
		tmp = get((mobjects[ii]))
		eval(parse(text=paste(snp,namei,"[rep,]","<-tmp",sep="") ))
		}
	}

	if(length(dims)==2) 			{
		tmp = get((mobjects[ii]))
		eval(parse(text=paste(snp,namei,"[rep,,]","<-tmp",sep="") ))
	}
	
	if(length(dims)==3) 			{
		tmp = get((mobjects[ii]))
		eval(parse(text=paste(snp,namei,"[rep,,,]","<-tmp",sep="") ))
	}
}
}