ulimit -s unlimited
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 "

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        =  2006
  Yr2        =  2099
  Nmon       =  12
  switch     =  1           ! =1: Use X(nx,nt) instead of X(nt,nx) ; =0: Origin X(nt,nx)
  standard   =  1           ! =1: Su dung ma tran Correlation;       =0: Su dung ma tran Covariance
  PC_Stard   =  1           ! =1: PCs chuan hoa (PCs=EVtors)         =0: EOFs chuan hoa
  Nmode      =  5           ! So mode duoc giu lai (Nmode <= Min{nt,nx})
/

EOF

./pca.exe

    done
done
echo "==================================="

MODELS=" CNRM-CM5 CSIRO-Mk3-6-0 EC-EARTH GFDL-ESM2M HadGEM2-AO IPSL-CM5A-LR MPI-ESM-MR Ens "

 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        =  1976
  Yr2        =  2005
  Nmon       =  12
  switch     =  1           ! =1: Use X(nx,nt) instead of X(nt,nx) ; =0: Origin X(nt,nx)
  standard   =  1           ! =1: Su dung ma tran Correlation;       =0: Su dung ma tran Covariance
  PC_Stard   =  1           ! =1: PCs chuan hoa (PCs=EVtors)         =0: EOFs chuan hoa
  Nmode      =  5           ! So mode duoc giu lai (Nmode <= Min{nt,nx})
/

EOF

./pca.exe

 done
echo "==================================="

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

EOF

./pca.exe

echo "All have been done!"