SUBROUTINE matrix_out (a, nrow, ncol,matfile)
USE Numeric_Kinds_Module
USE CommonInfo_Module
IMPLICIT NONE

! This routine outputs the input matrix to the file matfile

CHARACTER(LEN=*) matfile
INTEGER nrow, ncol, irow, jcol, istat
REAL(KIND=dp) a(nrow,ncol)

irow = 1

OPEN(UNIT=output_unit0,FILE=TRIM(outdir)//TRIM(matfile),IOSTAT=istat,STATUS='unknown',ACTION='write')
IF (istat.ne.0) STOP 'OPEN failed - matrix_out'
WRITE(output_unit0,230)(jcol,jcol=1,min(ncol,13))
WRITE(output_unit0,240)irow, (a(1,jcol),jcol=1,min(ncol,13))

DO irow = 2, min(nrow,13)
   WRITE(output_unit0,240)irow, (a(irow,jcol),jcol=1,min(ncol,13))
END DO

 230  FORMAT (1h0, 3x, 8h  row   , 13(3x, i3, 6x))
 240  FORMAT (1h , 1x, 4hcol , i3, 2x, 13(1pe12.4))
 250  FORMAT (1h0, 1x, 4hcol , i3, 2x, 13(1pe12.4))

CLOSE(UNIT=output_unit0,IOSTAT=istat,STATUS='keep')
IF (istat.ne.0) STOP 'CLOSE failed - matrix_out'

END SUBROUTINE matrix_out