Subroutine Read_Data_test
integer :: it, imon, ist
Real :: xlat0, xrain(1:Ntime), xtemp(1:Ntime)

!Real, Dimension (Nst_Max) :: Lat, Lon
!Real, Allocatable, Dimension (:) :: XXLat, XXLon
!Real, Dimension (Ntime, Nst_Max) :: DatBL, DatTMP


Nst=20

Open (1,file="/home/tanpv/SEACLID/.4ME/SPEI/fortran/tampa.txt", status="old")
Read(1,*); 
Read (1,*) xlat0;
Read(1,*); Read(1,*);
Do it=1, Ntime
   Read (1,*) xrain(it), xtemp(it)
   Do ist=1, Nst
      DatBL(it,ist)=xrain(it);
      DatTMP(it,ist)=xtemp(it)
      Lat(ist)=xlat0
   Enddo
Enddo

End Subroutine Read_Data_test
!-----------------------------
Subroutine Read_Data_test_SPI
integer :: it, imon, ist
Real :: xlat0, xrain(1:Ntime)

!Real, Dimension (Nst_Max) :: Lat, Lon
!Real, Allocatable, Dimension (:) :: XXLat, XXLon
!Real, Dimension (Ntime, Nst_Max) :: DatBL, DatTMP


Nst=20

Open (1,file="/home/tanpv/PDSI/Luu/SPI/src_new/stdin", status="old")
Read(1,*);
Do it=1, Ntime
   Read (1,*) xlat0, xlat0, xrain(it)
   Do ist=1, Nst
      DatBL(it,ist)=xrain(it);
   Enddo
Enddo

End Subroutine Read_Data_test_SPI
!-----------------------------
Subroutine Write_Results_test
Integer :: ilen, ist, it
Open (3,file="Test_origin_data.txt", Status="unknown")
Do ilen=1, Nlen
   write (3,*)," Timescale: ", Len_SPI(ilen), " months"
   write (3,"(A4,50I8)") "Stn", (ist,ist=1,Nst)
   Do it=1, Ntime !ist=1, Nst
   !   Do it=1, Ntime
      write (3,"(I4,50F8.4)") it, (SPI(ist,it,ilen), ist=1,Nst)
   !   Enddo
   Enddo
Enddo
End Subroutine Write_Results_test