!------------------------------------------------------------------------------
Subroutine Read_PDSI_Stn (FName)
Character (Len=*) :: FName
Real, Allocatable, Dimension (:) :: tmplon, tmplat, tmpTime, XStn
Integer :: Nlatmp, Nlotmp, Ntmp, Nst, Ntime
Character (Len=50) :: VarName

!VarName="ori_pdsi"
!  wr_ncdf_2d (trim(Fname), XStn, t, Nst, Ntime, X1, "ori_pdsi", "ori_pdsi", "-", "Station", "Time", amssng)
Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, PDSI, "ori_pdsi")
!Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, PDSI, "sel_pdsi")
Deallocate (XStn, tmpTime)
!Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, WPLM, "ori_wplm")
!Deallocate (XStn, tmpTime)
!Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, ZIND, "ori_zind")
!Deallocate (XStn, tmpTime)
!Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, scPDSI, "sel_pdsi")
!Deallocate (XStn, tmpTime)
!Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, scWPLM, "sel_wplm")
!Deallocate (XStn, tmpTime)
!Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, scZIND, "sel_zind")
!Deallocate (XStn, tmpTime)

Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, T2m, "T2m")
Deallocate (XStn, tmpTime)
Call rd_ncdf_2d (trim(FName), XStn, tmpTime, Nst, Ntime, Prec, "Pre")
Deallocate (XStn, tmpTime)

Call rd_ncdf_1d (trim(FName), XStn, Nst, Reg, "reg")
Deallocate (XStn)


!print*,"Num of Stn in "//trim(FName)//": ",Nst, PDSI

! Write out IDX to Text file
!Call Write_IDX_To_TEXT (PDSI, Nst, Ntime, Lon, Lat, Reg, "PDSI", "The PDSI drought index", "PDSI.txt")

End Subroutine Read_PDSI_Stn 
!--------------------------------------------------------------------------------