!-----------------------------------------------------------------------------
Subroutine Calc_Geographic_Extent (ifile)
Integer :: ifile, it, ilon, ilat, ireg, idrou1,idrou2,idrou3, itotal
! Area_Reg1(Nreg,Ntime)
!Real, dimension (Nreg+1,Ntime,Nfile) :: Geo_D1, Geo_D2, Geo_D3                     ! Spatial Extent 7 subregion + Entire Country
!Geo_D1=Miss; Geo_D2=Miss; Geo_D3=Miss

Do ireg=1, Nreg+1
   if (ireg<=Nreg) then
      Do it=1, Ntime
         idrou1=0; idrou2=0; idrou3=0; itotal=0     ! Tong so tram bi han va Tong so tram duoc xet trong tung spell
         Do ilon=1, Nlon
            Do ilat=1, Nlat
               if (Int(Reg(ilon,ilat))==ireg) then
                  itotal=itotal+1        ! Tong so tram duoc xet trong 1 thang  
                  if (D1(ilon,ilat,it)==1) idrou1=idrou1+1   ! Tong so tram bi anh huong
                  if (D2(ilon,ilat,it)==1) idrou2=idrou2+1   ! Tong so tram bi anh huong
                  if (D3(ilon,ilat,it)==1) idrou3=idrou3+1   ! Tong so tram bi anh huong
               endif
            Enddo
         Enddo
         Geo_D1(ireg,it,ifile)=real(idrou1)/real(itotal)*100.   ! Phan tram so tram bi anh huong
         Geo_D2(ireg,it,ifile)=real(idrou2)/real(itotal)*100.   ! Phan tram so tram bi anh huong
         Geo_D3(ireg,it,ifile)=real(idrou3)/real(itotal)*100.   ! Phan tram so tram bi anh huong
      Enddo
   else if (ireg==Nreg+1) then
      Do it=1, Ntime
         idrou1=0; idrou2=0; idrou3=0; itotal=0     ! Tong so tram bi han va Tong so tram duoc xet trong tung spell
         Do ilon=1, Nlon
            Do ilat=1, Nlat
               if (Int(Reg(ilon,ilat))>=1.AND.Int(Reg(ilon,ilat))<=Nreg) then
                  itotal=itotal+1        ! Tong so tram duoc xet trong 1 thang  
                  if (D1(ilon,ilat,it)==1) idrou1=idrou1+1   ! Tong so tram bi anh huong
                  if (D2(ilon,ilat,it)==1) idrou2=idrou2+1   ! Tong so tram bi anh huong
                  if (D3(ilon,ilat,it)==1) idrou3=idrou3+1   ! Tong so tram bi anh huong
               endif
            Enddo
         Enddo
         Geo_D1(ireg,it,ifile)=real(idrou1)/real(itotal)*100.   ! Phan tram so tram bi anh huong
         Geo_D2(ireg,it,ifile)=real(idrou2)/real(itotal)*100.   ! Phan tram so tram bi anh huong
         Geo_D3(ireg,it,ifile)=real(idrou3)/real(itotal)*100.   ! Phan tram so tram bi anh huong
      Enddo
   endif
Enddo
End Subroutine Calc_Geographic_Extent 
!-----------------------------------------------------------------------------