!-----------------------------------------------------------------------------
Subroutine Geographic_Extent_1
Integer :: it, ilon, ilat, ireg, idrou1,idrou2,idrou3, itotal
! Area_Reg1(Nreg,Ntime)
Area_Reg1=Miss; Area_Reg2=Miss; Area_Reg3=Miss
Do ireg=1, Nreg
   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
      Area_Reg1(ireg,it)=real(idrou1)/real(itotal)*100.   ! Phan tram so tram bi anh huong
      Area_Reg2(ireg,it)=real(idrou2)/real(itotal)*100.   ! Phan tram so tram bi anh huong
      Area_Reg3(ireg,it)=real(idrou3)/real(itotal)*100.   ! Phan tram so tram bi anh huong
   Enddo
Enddo
End Subroutine Geographic_Extent_1
!-----------------------------------------------------------------------------