c**********************************************************************
c     subroutine sb1out
c**********************************************************************
      subroutine sb1out
     i (anemht, wzoflg, kbr,jj,
     i  awu, 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, jj
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     include 'b1geom.inc'
      include 'b1glob.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 + + +

      integer outflag,m, n, k, icsr
      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 sb1out'
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,kbr, jj, awu
c
c     set output increment
      m = (imax - 1)/8
      m = max0(m,1)
      n = 1
c
      write (o_unit,*)  'sb1out output'
c      write (o_unit,*) 'for prior wind speed'
c      write (o_unit,21) (egt(k,n),k=1,(imax-1),m)
c      write (o_unit,22) (egtss(k,n),k=1,(imax-1),m)
c      write (o_unit,23??) ((egtss(k,n)/(egt(k,n)+0.0001)),k=1,(imax-1),m)
c      write (o_unit,*)

      write (o_unit,18)  (k , k=1,(imax-1),m), n

      write (o_unit,23)  (sf10(k,n),k=1,(imax-1),m)
      write (o_unit,24)  (sf84(k,n),k=1,(imax-1),m)
      write (o_unit,35)  (sf200(k,n),k=1,(imax-1),m)

      write (o_unit,25)  (szcr(k,n),k=1,(imax-1),m)
      write (o_unit,26)  (sfcr(k,n),k=1,(imax-1),m)
      write (o_unit,27)  (smlos(k,n),k=1,(imax-1),m)
      write (o_unit,28)  (sflos(k,n),k=1,(imax-1),m)

      write (o_unit,29)  (szrgh(k,n),k=1,(imax-1),m)
      write (o_unit,30)  (slrr(k,n),k=1,(imax-1),m)
      write (o_unit,38)  abzht(icsr),abrlai(icsr),
     &                    abrsai(icsr), abffcv(icsr)

      write (o_unit,31)  ahrwcw(1,1), ahrwc0(12,1)


      write (o_unit,32)  (wus(k,n),k=1,(imax-1),m)
      write (o_unit,33)  (wusp(k,n),k=1,(imax-1),m)
      write (o_unit,34)  (wust(k,n),k=1,(imax-1),m)
      write (o_unit,*)




c

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

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

   23 format (1x, 'sf10= ', 20f6.2)
   24 format (1x, 'sf84= ', 20f6.2)
   35 format (1x, 'sf200=', 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)
   38 format (1x, 'abzht=', f6.2, '  abrlai=', f4.2,
     &             '  abrsai=', f4.2, '  abffcv=',f4.2)
   31 format (1x, 'ahrwcw=',f4.2,'  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++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++