Data_DIR=/home/tanpv/Progs/Drou_SEA/EOF_New/NC/
PDF_DIR=/home/tanpv/Progs/Drou_SEA/EOF_New/PDF/
# Plot Scenarios:
SCENARIOS=" SCE_rcp45_ SCE_rcp85_ "
MODELS=" CNRM-CM5 CSIRO-Mk3-6-0 EC-EARTH GFDL-ESM2M HadGEM2-AO IPSL-CM5A-LR MPI-ESM-MR Ens "
NT=1116 
Y1=2007 
Y2=2099

for sce in $SCENARIOS ; do
   BB=`echo $sce |cut -c 5-9`
   for model in $MODELS ; do
      cat > tt.gs << EOF
   "reinit"

    MODEL="${model}"
    PCA_File="${Data_DIR}pca_${sce}${model}"
    EOF_File="${Data_DIR}eof_${sce}${model}"
    PDF_PCAs="${PDF_DIR}pca_${sce}${model}"
    PDF_EOFs="${PDF_DIR}eof_${sce}${model}"
    Ntime=${NT}; Yr1=${Y1}; Yr2=${Y2}
    SCENA="Scenario ${BB}"
    PERIOD="${Y1}-${Y2}"
    Vrange1=-0.1; Vrange2=0.1; YInt=0.05
EOF
       cat tt.gs PL_PCA_EOF_Temp.gs >> tt.gs
       grads -pcb tt.gs

    done
done
# Plot Baseline:
SCENARIOS=" BL_ "
MODELS=" CNRM-CM5 CSIRO-Mk3-6-0 EC-EARTH GFDL-ESM2M HadGEM2-AO IPSL-CM5A-LR MPI-ESM-MR Ens "
NT=348
Y1=1977
Y2=2005

for sce in $SCENARIOS ; do
   for model in $MODELS ; do
      cat > tt.gs << EOF
   "reinit"

    MODEL="${model}"
    PCA_File="${Data_DIR}pca_${sce}${model}"
    EOF_File="${Data_DIR}eof_${sce}${model}"
    PDF_PCAs="${PDF_DIR}pca_${sce}${model}"
    PDF_EOFs="${PDF_DIR}eof_${sce}${model}"
    Ntime=${NT}; Yr1=${Y1}; Yr2=${Y2}
    SCENA="Baseline"
    PERIOD="${Y1}-${Y2}"
    Vrange1=-0.2; Vrange2=0.2; YInt=0.1
*    Vrange1=-6; Vrange2=6; YInt=2
EOF
       cat tt.gs PL_PCA_EOF_Temp.gs >> tt.gs
       grads -pcb tt.gs

    done
done
# Plot APH:
MODELS=" APH " 
NT=348
Y1=1977
Y2=2005

 for model in $MODELS ; do
      cat > tt.gs << EOF
   "reinit"

    MODEL="${model}"
    PCA_File="${Data_DIR}pca_${model}"
    EOF_File="${Data_DIR}eof_${model}"
    PDF_PCAs="${PDF_DIR}pca_${model}"
    PDF_EOFs="${PDF_DIR}eof_${model}"
    Ntime=${NT}; Yr1=${Y1}; Yr2=${Y2}
    SCENA="APHRODITE data"
    PERIOD="${Y1}-${Y2}"
    Vrange1=-0.2; Vrange2=0.2; YInt=0.1
*    Vrange1=-6; Vrange2=6; YInt=2
EOF
       cat tt.gs PL_PCA_EOF_Temp.gs >> tt.gs
       grads -pcb tt.gs

 done


echo "==================================="