!$Author: jhudd@weru.ksu.edu $
!$Date: 2008-08-27 09:26:46 -0500 (Wed, 27 Aug 2008) $
!$Revision: 9689 $
!$HeadURL: https://svn.weru.ksu.edu/weru/weps1/trunk/weps.src/erosion/test/plotout.for $
!**********************************************************************
!     subroutine plotout.for
!**********************************************************************

      subroutine plotout (yplot, ycharin, yin)

!     + + +  PURPOSE + + +
!     1. to create headings for tsterode.eplt file
!     2. to store dep var (yin) and indep var (xin)
!         and write to tsterode.eplt file for each eros run

!     plotout is called from erodout.for with yin data
!     the xin data come from common/plot/ with erodin.for

!     + + + ARGUMENT DECLARATAIONS + + +
      integer yplot
      character*12 ycharin(30)
      real yin(30)

!     + + + ARGUMENT DEFINITIONS + + +
!     yplot     = number of dep variables to put in tsterode.eplt file
!     ycharin(i)= name(s) of dep. variables
!     yin(i)    = value(s) of dep. variables
!
!     + + + PARAMETERS + + +
!
!     + + + GLOBAL COMMON BLOCKS + + +
!
!     + + + LOCAL COMMON BLOCKS + + +

      integer xplot
      character*12 xcharin(30)
      real xin(30)
      common /plot/ xplot, xcharin, xin
!
!     + + + LOCAL VARIABLES + + +
!
      integer i,j, nline
      logical used
      character*500 line, plotdat(500)
!
!     + + + LOCAL VARIABLE DEFINITIONS + + +
!     i,j    = temp indexes
!     used   = logical for presence of tsterode.eplt file
!
!     + + + SUBROUTINES CALLED + + +
!
!     + + + FORMATS + + +
  200 format(40f12.4)
  201 format(' file:')
  202 format('tsterode.eplt')
!
!     + + + END SPECIFICATIONS + + +
! ^^^ tmp out
!      write(*,*)'^^^'
!       write(*,*)'out from plotout.for'
!       write(*,*) 'yplot=', yplot
!       write(*,*) 'ycharin=', ycharin(1:yplot)
!       write(*,*) 'xplot=', xplot
!       write(*,*) 'xcharin=', xcharin(1:xplot)
!       write(*,*) 'xin=',  xin(1:xplot)
!       write(*,*) 'yin=', yin(1:yplot)
! ^^^ end tmp out
!
!     create tsterode.eplt file
   !

!     update the tsterode.eplt file
   50 nline=j
      rewind (UNIT=80)
      do 55 j=1,nline
         write(80,'(a)') plotdat(j)
   55 continue
!     change sign of erosion components (yin) and write variables
      write (80,200)  (-1)*yin(1:yplot), xin(1:xplot)
      close (UNIT=80)

      return
      end
!**********************************************************************