# Calculate drought indices using observation data (VnGP+APH):
#make clean
#make aph
#./aph.out

# Calculate drought indices for the baseline period using data from 6 models
#make clean
#make bl
#./bl.out

# Calculate drought indices for the baseline period using ensemble data from 6 models
#make clean
#make blens
#./blens.out

# Calculate drought indices for the future periods and 2 scenarios using data from 6 models
make clean
make sce
make scens

PER=" 1 2 "
for peri in ${PER} ; do
    if [ ${peri} -eq 1 ] ; then
       YR1=2046
       YR2=2065
       PERIOD=MF
    else
       YR1=2080
       YR2=2099
       PERIOD=FF
    fi
    #echo $PERIOD " " $YR1 " " $YR2
    cat > Namelist_SCE.txt << EOF

&File_Dir_Names
   PDSI_Dir     = "/home/tanpv/SEACLID_New/Dr_VN_2020/my_result"       ! Thu muc chua KQ tinh PDSI cho tung tram (tung lan tinh)
   PDSI_WDir    = "/home/tanpv/SEACLID_New/Dr_VN_2020/my_dat"          ! Thu muc chua SL xu ly trung gian khi tinh PDSI
   Data_dir     = "/home/tanpv/SEACLID_New/Data/"                      ! Data directory
   Period       = ${PERIOD}                     ! Must be: YBas1=2080, YBas2=2099
   Model_nm     = "CNRM", "CSIR", "ECEA", "GFDL", "HadG", "MPI"
   MaskFname    = "/work/projects/G09/SEACLID/Land_Mask_VN.nc"
   AWC_Fname    = "/work/projects/G09/SEACLID/AWC_VN_FAO.nc"           ! Ten file SL AWC_VN noi suy tu FAO
   Out_Dir      = "/home/tanpv/SEACLID_New/Dr_VN_2020/Results/"        ! Output directory
   TVar         = "tas"             ! Ten bien nhiet do trong file NC
   RVar         = "pr"              ! Ten bien luong mua trong file NC
/

&Data_Info
   YBas1        = ${YR1}              ! Nam dau tien cua chuoi SL T & R
   YBas2        = ${YR2}              ! Nam cuoi cung cua chuoi SL T & R
   NLen         = 3                 ! So quy mo thoi gian han (Drought timescale)
/

&Dr_Timescale
   Len_SPI      = 3, 6, 12          ! Cac drought timescale (thang)
/

EOF
    # Run program:
   ./sce.out
   ./scens.out

done