Subroutine Write_3D_Frequancy (Fname)
Integer :: ifile,ibin,ireg
Character (Len=*) :: Fname 
Real :: Xtime(Nfile)

Do ifile=1, Nfile; Xtime(ifile)=ifile; Enddo      ! Mang thoi gian gia tao, Don vi la "Dat_Type"
Do ibin=1, Nbin; Xbin(ibin)=(BIN(ibin)+BIN(ibin+1))/2.; Enddo
Do ireg=1, Nreg+1; Xreg(ireg)=ireg; Enddo

!Call wr_ncdf_3d(trim(Fname),Lon,Lat,Xtime, Nlon,Nlat,Nfile, T2m_Ann,"T2m_Ann","months since 1980-01-15 ","Annual mean temperature","C-deg","Lon","Lat","Dat_Type",Miss)
!Call wr_ncdf_3d_apd(trim(Fname), Pre_Ann,"Pre_Ann", Nlon,Nlat,Nfile, "Annual mean precipitation", "mm/day", "Lon","Lat","Dat_Type",Miss)

Call wr_ncdf_3d(trim(Fname),Xreg,Xbin,Xtime, Nreg+1,Nbin,Nfile, Freq,"Freq","months since 1980-01-15 ","Frequency Distribution of Drought IDX","%","Subreg","BIN","Dat_Type",Miss)

End Subroutine Write_3D_Frequancy
!-------------------------------------