!----------------------------------------------------------------------------------!
Subroutine Read_Clim_IDX
Integer, Parameter :: N_ClimIDX=16
Character (Len=50), Parameter :: Fname_idx(N_ClimIDX) = &
!                        1        2          3          4            5          6        7         8
!                   (/"DMI.txt","BEST.txt","MEI.txt","Nino3.4.txt","ONI.txt","SOI.txt","WP.txt","WPC850.txt", &
!                        9          10                 11         12            13     14       15           16  
!                     "QBO.txt","RINDO_SLPA.txt","CPAC850.txt","PDO.txt","TNI.txt","PNA.txt","REQSOI.txt","WHWP.txt"/)
!                        1        2          3             4               5            6             7            8
                   (/"ONI.txt", "SOI.txt","Nino3.4.txt","MEI.txt","RINDO_SLPA.txt","CPAC850.txt","REQSOI.txt","WPC850.txt", &
                     "BEST.txt","PDO.txt","DMI.txt",    "TNI.txt","PNA.txt",       "WHWP.txt",    "WP.txt",   "QBO.txt"/)
Integer :: iidx, iyr, imon, year, IO, it
Character (Len=50) :: Fname
Real :: X(Nmon)
Logical :: Exists
!Allocate(ClimIDX(Nidx, Nyear, Nmon), Dat_IDX(Ntime, Nidx))
! X1(Nmax, Mmax, Mmon)

!----------------------------------------------------------------------------
CL_Name(1)="ONI"; CL_Name(2)= "SOI"; CL_Name(3) ="Nino3.4"; CL_Name(4) ="MEI"; CL_Name(5)= "RINDO_SLPA"; CL_Name(6)="CPAC850"; CL_Name(7)="REQSOI"; CL_Name(8)="WPC850"
CL_Name(9)="BEST";CL_Name(10)="PDO"; CL_Name(11)="DMI";     CL_Name(12)="TNI"; CL_Name(13)="PNA";       CL_Name(14)="WHWP";   CL_Name(15)="WP";    CL_Name(16)="QBO"

!CL_Name(1)="DMI"; CL_Name(2)="BEST"; CL_Name(3)="MEI"; CL_Name(4)="Nino3.4"; CL_Name(5)="ONI"; CL_Name(6)="SOI"; CL_Name(7)="WP"; CL_Name(8)="WPC850"
!CL_Name(9)="QBO"; CL_Name(10)="RINDO_SLPA"; CL_Name(11)="CPAC850"; CL_Name(12)="PDO"; CL_Name(13)="TNI"; CL_Name(14)="PNA"; CL_Name(15)="REQSOI"; CL_Name(16)="WHWP"
!----------------------------------------------------------------------------


Do iidx=1, N_ClimIDX
   Fname="/home/tanpv/Nafosted/Data/ClimIDX/"//trim(Fname_idx(iidx))
   Inquire(file=trim(Fname), Exist=Exists)
   if (.Not.Exists) then
      print*,"File ",trim(Fname)," does not exist"
      cycle
   endif
   Open(1, file=trim(Fname), Status="Old")
   read(1,*)         ! Bo dong dau
   it=0
   Do
      read(1,*,IOStat=IO) year, (X(imon),imon=1,Nmon)
      if (IO /= 0) exit
      Do iyr=Yr1, Yr2
         if (year==iyr) then 
            YY(iyr-Yr1+1) = iyr 
            Do imon=1, Nmon
               it=it+1
               ClimIDX(iidx,it)=X(imon) 
            Enddo
         endif
      Enddo
   Enddo
   Close (1)
Enddo

End Subroutine Read_Clim_IDX
!-----------------------------------------------------------------------