cd /home/tanpv/Progs/562_VNU/GS

LEADS='c("02")'
MODEL=RegCM
MEMS='c("m5")'
VAR=pr
cat > my_bpl_1.R << EOF
library(ncdf4); library(ggplot2)
library(reshape2); library(ggpubr)
clead <- ${LEADS} # c("04") list of leastimes
model <- "${MODEL}" # "RegCM"
c.mem <- ${MEMS}
idir <- "/home/tanpv/Data/562_VNU/0.25_VN/Dat_Reg/"
odir <- "/home/tanpv/Progs/562_VNU/PDF/"
for (lead in c(clead)){
    ifile <- paste("${VAR}_",model,"_lead_",lead,sep="")

EOF
cat my_bpl.R >> my_bpl_1.R

Rscript my_bpl_1.R