!-----------------------------------------------------------------------------
Subroutine Cal_Freq_All_Timeseries(Dr_IDX, Models)
Integer :: it, ilon, ilat, ireg, ibin, ic, N
Real :: tmp
Character (Len=*) :: Dr_IDX, Models
Do ireg=1, Nreg
   Do ibin=1, Nbin
      ic=0; N=0
      Do ilon=1, Nlon
         Do ilat=1, Nlat
            if (Int(Reg(ilon,ilat))==ireg) then
               Do it=1, Ntime
                  if (PDSI(ilon,ilat,it) /= Miss) then
                     N=N+1
                     if (PDSI(ilon,ilat,it)>=BIN(ibin) .AND. PDSI(ilon,ilat,it)<BIN(ibin+1)) ic=ic+1
                  endif
               Enddo
            endif
         Enddo
      Enddo
      Freq(ireg,ibin) = Real(ic)/Real(N)*100; Xbin(ibin)=(BIN(ibin)+BIN(ibin+1))/2.
   !   print "(I2,2F6.2)",ireg, Xbin(ibin), Freq(ireg,ibin)
   Enddo
Enddo
! write to netCDF file:
!Fname="/home/tanpv/SEACLID_New/Post_Pro_Stn_No_Rain_Center/PDSI_Stat.nc"
Fname=trim(Models)//trim(Dr_IDX)//"_Reg_Fre_Distribution.nc"
Do ireg=1, Nreg; Xreg(ireg)=ireg; Enddo
Call wr_ncdf_2d (trim(Fname), Xreg, Xbin, Nreg, Nbin, Freq, "Freq", "Freq_Drought_Index", "-", "Sub-Region", "BIN", Miss)

End Subroutine Cal_Freq_All_Timeseries
!-----------------------------------------------------------------------------