WriteMatrixFloat.F

c**********************************************************************
#include "author.inc"
c*    $Id: WriteMatrixFloat.F,v 1.6 1995/09/25 03:20:10 turner Exp $
c*
c*    Writes out matrix of floating point numbers.
c*
c*    <PARAMETER LIST>
c*
c*     Input:
c*      lu - logical unit to which to write
c*      idim - leading dimension of array a
c*      nrows - number of active rows in array a
c*      ncols - number of active columns in array a
c*      a - matrix
c*      title - character description
c*
c*     Output:
c*      status - return status
c*
#include "copyright.inc"
c**********************************************************************
      subroutine JT_WriteMatrixFloat (lu, idim, nrows, ncols, a,
     &           title, status)
      implicit none
c
c ... Input:
      integer idim, lu, nrows, ncols
      real a(idim,ncols)
      character*(*) title
c
c ... Output:
      integer status
c
c ... Local:
      integer i, j
c
      write(lu,*)
      write(lu,*) title
      write(lu,*)
      if (ncols .le. 7) then
       do i=1,nrows
        write(lu,500) i, (a(i,j),j=1,ncols)
       enddo
      else
       do j=1,ncols
        write(lu,*)
        write(lu,*) 'Column: ',j
        write(lu,*)
        do i=1,nrows
         write(lu,500) i, a(i,j)
        enddo
       enddo
      endif
c
      status = 0
      return
c
  500 format (t2, i7, 1p, 7(2x,e13.6))
      end