!-----------------------------------------------------------------------------
Subroutine Define_Drought_Events
! D1: Mild or more intense drought, D2: Moderate more intense drought, D3: Severe of more intense drought, D4: Extreme drought
!Real, Dimension(Nlon,Nlat,N_thre,Nmodels+2) :: DThre          ! Nguong chi so han cua OBS va SL Baseline 6 mo hinh + 1 Ens
Integer :: ilon, ilat, ireg, it, DR1, DR0
D1=Miss; D2=Miss; D3=Miss; OD1=Miss; OD2=Miss; OD3=Miss
Do ilon=1, Nlon
   Do ilat=1, Nlat
      if (Reg(ilon,ilat)>=1.AND.Reg(ilon,ilat)<=Nreg) then
         Do it=1, Ntime
            if (PDSI(ilon,ilat,it)/=Miss.AND.PDSI(ilon,ilat,it)<=Thre_1) then; D1(ilon,ilat,it)=1; OD1(ilon,ilat,it)=0; else; D1(ilon,ilat,it)=0; OD1(ilon,ilat,it)=1; endif
            if (PDSI(ilon,ilat,it)/=Miss.AND.PDSI(ilon,ilat,it)<=Thre_2) then; D2(ilon,ilat,it)=1; OD2(ilon,ilat,it)=0; else; D2(ilon,ilat,it)=0; OD2(ilon,ilat,it)=1; endif
            if (PDSI(ilon,ilat,it)/=Miss.AND.PDSI(ilon,ilat,it)<=Thre_3) then; D3(ilon,ilat,it)=1; OD3(ilon,ilat,it)=0; else; D3(ilon,ilat,it)=0; OD3(ilon,ilat,it)=1; endif
         Enddo
      endif
   Enddo
Enddo
End Subroutine Define_Drought_Events
!-----------------------------------------------------------------------------