!-----------------------------------------------------------------------------
Subroutine Drought_Index_Areal_Mean
Integer :: ireg, it
Do it=1, Ntime
   Do ireg=1, Nreg
      PDSI_Reg(it,ireg)=Area_Mean (PDSI, ireg, it)
   Enddo
Enddo
End Subroutine Drought_Index_Areal_Mean
!-----------------------------------------------------------------------------