Subroutine Read_data_origin
! Real, Dimension (NVar, Nst_Max, YBas1:YBas2, NMonth) :: Dat_Mon !
! ivar=1: T2m ivar=2: Rain
! AWC(ist), xlat(ist)
Integer :: ist, iyr, imon, year
Real :: Lat_0, AWC_0, Xtmp(YBas1:YBas2, 1:NMonth), Xrain(YBas1:YBas2, 1:NMonth)

Nst=10   ! Test with 10 Stations
Open (1,file="/home/tanpv/PDSI/Luu/PDSI/example/monthly_T", Status="old");
Open (2,file="/home/tanpv/PDSI/Luu/PDSI/example/monthly_P", Status="old");
Open (3,file="/home/tanpv/PDSI/Luu/PDSI/example/parameter", Status="old");

Read (3,*) AWC_0, Lat_0; Close (3)

Do iyr=YBas1, YBas2
   Read (1,*) year, (Xtmp(iyr,imon), imon=1,NMonth)
   Read (2,*) year, (Xrain(iyr,imon), imon=1,NMonth)
Enddo
Close (1); Close (2)
Do ist=1, Nst
   AWC(ist)=AWC_0; xlat(ist)=Lat_0
   Dat_Mon(1,ist,YBas1:YBas2,1:NMonth)=Xtmp(YBas1:YBas2,1:NMonth)
   Dat_Mon(2,ist,YBas1:YBas2,1:NMonth)=Xrain(YBas1:YBas2,1:NMonth)
Enddo
End Subroutine Read_data_origin
!--------------------------------
Subroutine Write_origin_results
Integer :: ist, iyr, imon, year, it
Real :: Xtmp(1:NMonth)

Open (3,file="Test_with_origin_data.txt", Status="unknown")
Do ist=1, Nst
   it=0
   Do iyr=YBas1, YBas2
      Do imon=1, NMonth
         it=it+1; Xtmp(imon)=mon_ori_pdsi(ist,it)
      Enddo
      write (3,"(I5,12F7.2)") iyr, (Xtmp(imon), imon=1,NMonth)
   Enddo
Enddo
Close (3)
End Subroutine Write_origin_results