!**********************************************************************
!     subroutine sb2out
!**********************************************************************
      subroutine sb2out (anemht, wzoflg, awu, kbr, o_unit)
!
!     + + + PURPOSE + + +
!     To print to file tst.out some key variables used in erosion
!     use wind direction of 270 to see output along downwind direction
!
!     + + + ARGUEMENT DECLARATIONS + + +
      real anemht, awu
      integer  wzoflg, kbr, o_unit
!
!     + + + ARGUMENT DEFINITIONS + + +
!     anemht = anemometer height (m)
!     awzzo =  aerodynamic roughness at anemometer (mm)
!     wzz0  =  aerodynamic roughness length (mm) 
!     awu   =  wind speed (m/s)
!     wus   =  friction velocity (m/s)
!     wust  =  threshold friction velocity (m/s)
!     wzoflg = flag to showing anemometer at field (1) or wx sta (0)
!     o_unit= Unit number for output file
!
!     + + + GLOBAL COMMON BLOCKS + + +

      include 'p1werm.inc'
      include 'h1db1.inc'
!
!     + + + LOCAL COMMON BLOCKS + + +
!

      include 'erosion/s2agg.inc'
      include 'erosion/s2surf.inc'
      include 'erosion/s2sgeo.inc'
      include 'erosion/w2wind.inc'
      include 'erosion/m2geo.inc'
      include 'erosion/e2erod.inc'
!
!     + + + LOCAL VARIABLES + + +

      real egavg(mngdpt)
      integer outflag, m, n, k, icsr, j
      save outflag

!     outflag = 0 - print heading output, 1 - no more heading

!     + + + END SPECIFICATIONS + + +
!
!     output headings?
      if (outflag .eq. 0) then
        write (o_unit,*)
        write (o_unit,*) 'OUT PUT from sb2out'
!        write (o_unit, 15)
        write (o_unit,*)
!     turn off heading output
        outflag = 1
      endif

!     define index of current subregions
      icsr = 1

      write (o_unit,20)  anemht,wzoflg,awu,kbr
!
!     set output increment
      m = (imax - 1)/8
      m = max0(m,1)
      n = 1
!
!     initialize avg erosion variable
      do 3 j = 1, imax
        egavg(j) = 0.0
    3 continue
!
!     calc. avg erosion over a given field length
      do 5 j = 1, (imax-1), m
      do 4 k = 1, j
         egavg(j) = egavg(j) + egt(k,n)
    4 continue
         egavg(j) = egavg(j)/j
    5 continue

      write (o_unit,*)  'sb2out output'

      write (o_unit,18)  (k , k=1,(imax-1),m), n
      write (o_unit,21) (egt(k,n),k=1,(imax-1),m)
      write (o_unit,22) (egtss(k,n),k=1,(imax-1),m)
      write (o_unit,23) ((egtss(k,n)/(egt(k,n)+0.00001)),               &
     &                   k=1,(imax-1),m)
      write (o_unit,35) (egavg(k), k=1,(imax-1),m)
      write (o_unit,*) '----------------------------------------------'

!



!

!
!     output formats
!   10 format (1x, 'anemht wzoflg   awu    kbr')
   18 format(1x, 'i..n,j', 20i6)
!   15 format (1x, ' (m)          (m/s)    ')
   20 format (1x, 'anemht wzoflg awu kbr',                              &
     &  f6.0, i6, f8.2, i6)

   21 format (1x, 'egt=   ', 20f6.2)
   22 format (1x, 'egtss= ', 20f6.2)
   23 format (1x, 'egtss/egt=', 20f6.2)
   35 format (1x, 'egavg = ', 20f6.2)

   24 format (1x, 'sf84=', 20f6.2)

   25 format (1x, 'szcr=', 20f6.2)
   26 format (1x, 'sfcr=', 20f6.2)
   27 format (1x, 'smlos=', 20f6.2)
   28 format (1x, 'sflos=', 20f6.2)

   29 format (1x, 'szrgh=', 20f6.2)
   30 format (1x, 'slrr=', 20f6.2)

   31 format (1x, 'ahrwc0(icsr,12)', f6.2)
   32 format (1x, 'wus=', 20f6.2)
   33 format (1x, 'wusp=', 20f6.2)
   34 format (1x, 'wust=', 20f6.2)
!
      return
      end
!++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++