!**********************************************************************
      subroutine erodinit
!
!     +++ PURPOSE +++
!
!     Controls calls to subroutines that:
!       create the Erosion submodel grid (sbgrid)    
!       initialize Erosion submodel output array to zero (sbigrd).
!       calculate normalized effect of hills on friction velocity 
!        on grid for each wind direction (not activated)
!       calculate normalized effect of barriers on friction velocity
!        on grid for each wind direction (sbbr) 
!     + + + GLOBAL COMMON BLOCKS + + +
!
      include  'p1werm.inc'
      include  'm1flag.inc'
      include  'm1geo.inc'
      
!
!     +++ SUBROUTINES CALLED +++
!
!     sbgrid
!     sbigrd
!     sbhill (not activated)
!     sbbr
!     +++ FUNCTION DECLARATIONS +++
!
!     +++ LOCAL VARIABLES +++
!
!    + + + LOCAL VARIABLE DEFINITIONS + + +
!     nbr  = number of barriers (from m1geo.inc)
!
!     +++ END SPECIFICATIONS +++
!
!     Grid is created at least once.
      if (am0eif .eqv. .true.) then
         call sbgrid
!
!     set grid cell output arrays to zero
         call sbigrd
!
!     check for hills - sbhill not implemented
!        if (nhill .gt. 0) then
!        call sbhill
!        endif
!
!     check for barriers
         if (nbr .gt. 0) then
         call sbbr
         endif
!
!     Turn off grid creation flag
      am0eif = .false.
      endif
      return
      end
!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++