MODELS=" RegCM " #  clWRF "
VARS=" pr t2 "
# ------------------------------------------------------------
cd /home/tanpv/Progs/562_VNU/GS
# ----------------------------------------------------------------------------
MEM_1=1
MEM_2=5
DMEM=2
MEMS='c("m1","m3","m5")'
LEADS='c("01","03","05")'

for model in ${MODELS} ; do
   for var in ${VARS} ; do
   # -----------------------------  Boxplots using R - Template file: my_bpl.R
if [ ${var} == pr ] ; then
  XMIN="c(0,0,0,0,0,0,0)"
  XMAX="c(1000,1000,1000,1000,1000,1000,1000)"
  XDEL="c(200,200,200,200,200,200,200)"
fi
if [ ${var} == t2 ] ; then
  XMIN="c(10,10,10,10,10,10,20)"
  XMAX="c(35,35,35,35,35,35,32)"
  XDEL="c(5,5,5,5,5,5,4)"
fi

cat > my_bpl_2.R << EOF
library(ncdf4); library(ggplot2)
library(reshape2); library(ggpubr)
clead <- ${LEADS} # c("04") list of leastimes
model <- "${model}" # "RegCM"
c.mem <- ${MEMS}
var <- "${var}"
xmins <- ${XMIN} 
xmaxs <- ${XMAX}
xdels <- ${XDEL}
idir <- "/home/tanpv/Data/562_VNU/0.25_VN/Dat_Reg/"
odir <- "/home/tanpv/Progs/562_VNU/PDF/"

EOF

cat my_bpl.R >> my_bpl_2.R

Rscript my_bpl_2.R

   # -----------------------------  Scatterplots - Template file: 4b.disp_Scatter_Template.gs 
   cat > 4b1.gs << EOF
   "reinit"
    Modelnm="${model}"; Ele="${var}"; Mem1=${MEM_1}; Mem2=${MEM_2}; Dmem=${DMEM} 
EOF
   cat 4b.disp_Scatter_Template.gs >> 4b1.gs
   # -----------------------------  Verification scores for Training - Template file: 3.veri_BC_ANN_Training_Template.gs
   cat > 3b1.gs << EOF
   "reinit"
   Model="${model}"; Varnm="${var}"
EOF
   cat 3.veri_BC_ANN_Training_Template.gs >> 3b1.gs 
   # -----------------------------  Verification scores for testing - Template file: 3.veri_BC_ANN_testing_Template.gs
   cat > 3b2.gs << EOF
   "reinit"
   Model="${model}"; Varnm="${var}"; Mem1=${MEM_1}; Mem2=${MEM_2}; Dmem=${DMEM}
EOF
   cat 3.veri_BC_ANN_testing_Template.gs >> 3b2.gs
   #######################################################

   grads -pbc 4b1.gs 
   grads -lbc 3b1.gs
   grads -lbc 3b2.gs

   done
done

#######################################################