c**********************************************************************
c     subroutine sb2out
c**********************************************************************
      subroutine sb2out
     i (anemht, wzoflg, awu, kbr,
     i  o_unit)
c
c     + + + PURPOSE + + +
c     To print to file tst.out some key variables used in erosion
c
c     + + + ARGUEMENT DECLARATIONS + + +
      real anemht, awu, wzoflg
      integer  kbr, o_unit
c
c     + + + ARGUMENT DEFINITIONS + + +
c     anemht =
c     awzzo =
c     wzz0  =
c     awu   =
c     wus   =
c     wust  =
c     o_unit= Unit number for output file
c
c     + + + GLOBAL COMMON BLOCKS + + +
*$noreference
      include 'p1werm.inc'
      include 'h1db1.inc'
c
c     + + + LOCAL COMMON BLOCKS + + +
c

      include 'erosion/s2agg.inc'
      include 'erosion/s2surf.inc'
      include 'erosion/s2sgeo.inc'
      include 'erosion/w2wind.inc'
      include 'erosion/m2geo.inc'
      include 'erosion/e2erod.inc'
*$reference
c
c     + + + LOCAL VARIABLES + + +

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

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

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

c     define index of current subregions
      icsr = 1

      write (o_unit,20)  anemht,wzoflg,awu,kbr
c
c     set output increment
      m = (imax - 1)/8
      m = max0(m,1)
      n = 1
c
c     initialize avg erosion variable
      do 3 j = 1, imax
        egavg(j) = 0.0
    3 continue
c
c     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,*) '----------------------------------------------'

c



c

c
c     output formats
c   10 format (1x, 'anemht wzoflg   awu    kbr')
   18 format(1x, 'i..n,j', 20i6)
c   15 format (1x, ' (m)          (m/s)    ')
   20 format (1x, 'anemht wzoflg awu kbr',
     &  2f6.0,  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)
c
      return
      end
c++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++