c**********************************************************************
c     subroutine sbwus
c**********************************************************************
      subroutine sbwus
     i (anemht, awzzo, awu, wzzo,
     o  wus)
c
c     +++ PURPOSE +++
c     To calculate subregion, friction velocity, given station
c     anemometer height, surface roughness, wind speed; and subregion
c     aerodynamic roughness.
c
c     +++ ARGUMENT DECLARATIONS +++
      real*4 anemht, awzzo, awu, wzzo, wus
c
c     +++ ARGUMENT DEFINITIONS +++
c     anemht - parameter, anemometer height of input wind speed (m).
c     awzzo - parameter, surface aerodynamic roughness at input wind
c             speed location (mm).
c     awu - input wind speed driving EROSION submodel (m/s).
c     wzzo - subregion aerodynamic roughness (mm).
c     wus - subregion soil surface friction velocity (m/s).
c
c     +++ LOCAL VARIABLES +++
      real*4 wusst
c
c     +++ END SPECIFICATIONS +++
c
c     Calc station (input wind speed location) friction velocity
      wusst = awu*0.4/alog(anemht*1000./awzzo)
c
c     calc subregion friction velocity
      wus = wusst * (wzzo/awzzo)**0.067
c
c     temp code
      return
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++