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