# Clean all:
rm -f *.exe
# Make executable file:
make pca
# Set unlimited memory:
ulimit -s unlimited
# 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 "
Y1=2006
Y2=2099
NMN=12
SW=1
STD=1
PCSTD=0
NMOD=5
for sce in $SCENARIOS ; do
    for model in $MODELS ; do
        cat > Namelist_PCA_New.txt << EOF
 &Model_Info
  Fname      = "/home/tanpv/Progs/Dr_SEA/Results/${sce}${model}.nc"
  MaskFname  = "/work/users/hmo/tanpv/CORDEX_SEACLID/Land_mask_sea.nc"
  pca_out    = "pca_${sce}${model}.nc"
  eof_out    = "eof_${sce}${model}.nc"
  Yr1        =  ${Y1}
  Yr2        =  ${Y2}
  Nmon       =  ${NMN}
  switch     =  ${SW}            ! =1: Use X(nx,nt) instead of X(nt,nx) ; =0: Origin X(nt,nx)
  standard   =  ${STD}           ! =1: Su dung ma tran Correlation;       =0: Su dung ma tran Covariance
  PC_Stard   =  ${PCSTD}         ! =1: PCs chuan hoa (PCs=EVtors)         =0: EOFs chuan hoa
  Nmode      =  ${NMOD}          ! So mode duoc giu lai (Nmode <= Min{nt,nx})
/

EOF

./pca.exe

    done
done
echo "==================================="
# Baseline (Models)
MODELS=" CNRM-CM5 CSIRO-Mk3-6-0 EC-EARTH GFDL-ESM2M HadGEM2-AO IPSL-CM5A-LR MPI-ESM-MR Ens "
Y1=1976
Y2=2005

 for model in $MODELS ; do
     cat > Namelist_PCA_New.txt << EOF
 &Model_Info
  Fname      = "/home/tanpv/Progs/Dr_SEA/Results/BL_${model}.nc"
  MaskFname  = "/work/users/hmo/tanpv/CORDEX_SEACLID/Land_mask_sea.nc"
  pca_out    = "pca_BL_${model}.nc"
  eof_out    = "eof_BL_${model}.nc"
  Yr1        =  ${Y1}
  Yr2        =  ${Y2}
  Nmon       =  ${NMN}
  switch     =  ${SW}            ! =1: Use X(nx,nt) instead of X(nt,nx) ; =0: Origin X(nt,nx)
  standard   =  ${STD}           ! =1: Su dung ma tran Correlation;       =0: Su dung ma tran Covariance
  PC_Stard   =  ${PCSTD}         ! =1: PCs chuan hoa (PCs=EVtors)         =0: EOFs chuan hoa
  Nmode      =  ${NMOD}          ! So mode duoc giu lai (Nmode <= Min{nt,nx})
/

EOF

./pca.exe

 done
echo "==================================="
# Baseline (APHRODITE)
###########################
Y1=1976
Y2=2005

cat > Namelist_PCA_New.txt << EOF
 &Model_Info
  Fname      = "/home/tanpv/Progs/Dr_SEA/Results/BL_APH.nc"
  MaskFname  = "/work/users/hmo/tanpv/CORDEX_SEACLID/Land_mask_sea.nc"
  pca_out    = "pca_APH.nc"
  eof_out    = "eof_APH.nc"
  Yr1        =  ${Y1}
  Yr2        =  ${Y2}
  Nmon       =  ${NMN}
  switch     =  ${SW}            ! =1: Use X(nx,nt) instead of X(nt,nx) ; =0: Origin X(nt,nx)
  standard   =  ${STD}           ! =1: Su dung ma tran Correlation;       =0: Su dung ma tran Covariance
  PC_Stard   =  ${PCSTD}         ! =1: PCs chuan hoa (PCs=EVtors)         =0: EOFs chuan hoa
  Nmode      =  ${NMOD}          ! So mode duoc giu lai (Nmode <= Min{nt,nx})
/

EOF

./pca.exe

echo "All have been done!"