{
  library(ncdf4)
  library(ggplot2)
  library(reshape2)
  library(ggpubr)
}

{
  #idir <- "/home/tanpv/Data/562_VNU/0.25_VN/Dat_Reg/"
  #odir <- "/home/tanpv/Progs/562_VNU/PDF/"
  idir <- c("/home/tanpv/Data/562_VNU/0.25_VN/Dat_Reg")
  odir <- c("/home/tanpv/Progs/562_VNU/PDF")

  clead <- c("05") # c("04") list of leastimes
  model <- "RegCM" # "RegCM"
  var <- "t2"
}

{
  xmins <- c(10,10,10,10,10,10,20)
  xmaxs <- c(35,35,35,35,35,35,30)
  xdels <- c(5,5,5,5,5,5,5)
}


for (lead in c(clead)){
  ifile <- paste0(var,"_",model,"_lead_",lead,"_ens")
  iname <- paste0(idir,"/", ifile,".nc")
  nc_file <- nc_open(iname)
  cvar.name <- c("obs","mod","cli","ann", "bjp", "q_m")    # Order lai, add more or remove
  boxcolor <- c("obs"="#FF0000","mod"="#ffff00","cli"="#0000FF","ann"="#009900", "bjp"="orange", "q_m"="dodgerblue")
  for (ivar in c(cvar.name)){
    tmp.cmd <- paste0(ivar,".dat<-ncvar_get(nc_file,'",ivar,"')")
    eval(parse(text=tmp.cmd))
  }
  
  nreg <- 7
  nmon <- 12
  itmp <- 0
  for (ivar in c(cvar.name)){
    itmp <- itmp + 1
    tmp.cmd <- paste0("tmp.dat<-",ivar,".dat")
    eval(parse(text=tmp.cmd))
    nsample <- dim(tmp.dat)[3]
    my.dat <- tmp.dat
    dimnames(my.dat) <- list(c(seq(7)),
                             c(seq(12)),
                             c(seq(nsample)))
    tmp.mel <- melt(my.dat)
    colnames(tmp.mel) <- c("reg","mon","sample","value")
    tmp.mel$dat <- ivar
    if (itmp == 1) { 
        odat <- tmp.mel
    } else {
        odat <- rbind(odat,tmp.mel)
    }
    remove(tmp.mel,my.dat)
  }
  
  odat.1 <- odat
  for (ireg in seq(nreg)){
    draw.dat <- subset(odat.1,odat.1$reg==ireg)
    a <- data.frame(draw.dat$value,draw.dat$mon,draw.dat$dat)
    colnames(a) <- c("value","mon","dat")
    a$mon <- factor(a$mon, levels = seq(12))
    a$dat <- as.character(a$dat)
    remove(p)
    a$dat <- factor(a$dat, levels=c(cvar.name))
    xmin <- xmins[ireg]
    xmax <- xmaxs[ireg]
    xdel <- xdels[ireg]
    p <- ggplot(a,aes(x=mon,y=value))+
        geom_boxplot(aes(fill=dat),
                     outlier.size=0.1,outlier.shape = 1,outlier.color ="gray50" )+
        scale_y_continuous(breaks = c(seq(xmin,xmax,xdel)),limits = c(xmin,xmax))+
        scale_fill_manual(name="",values = boxcolor)+ 
        ggtitle(paste0("R",ireg,"-",model,"-LT=",lead)) + xlab("") + ylab("(mm)") + labs(x=NULL)
    tmp.cmd <- paste("p",ireg,"<-p",sep="")
    eval(parse(text=tmp.cmd))+theme(legend.position = "none")
  }
  
  save.image <- ggarrange(p1,p2,p3,p4,p5,p6,p7,ncol = 1, nrow = 7,
                          common.legend = TRUE, legend="bottom")
  save.image
  
  png_file <- paste0(odir,"/",ifile,"-","_ens.png")
  png(png_file, width = 8.0, height = 14.0, units = 'in', res = 300)
  
  #  pdf(paste(odir,ifile,".pdf",sep=""), width = 4.5, height = 14.0)
  #  plot(save.image) 
  print(save.image)
  dev.off()
  remove(save.image)
}