SUBROUTINE dmat_prt (name, dmat, nchn, nq, ndim)
USE Numeric_Kinds_Module
USE FileUnits_Module
IMPLICIT NONE
!==============================================================================
LOGICAL, PARAMETER :: debug=.False.
CHARACTER (len=*) :: name
INTEGER :: ix, iq, nchn, nq, ndim
REAL (dp) :: dmat (ndim, nq)
!==============================================================================
! Simple matrix printing routine.

IF (.not.debug) RETURN
WRITE (61, * )
WRITE (61, * ) 'Matrix ', name
WRITE (61,  * ) '     Actual dimension=', nchn, ' by ', nq
DO ix = 1, min (nchn, 10)
WRITE (61, 10) (dmat (ix, iq), iq = 1, min (nq, 10) )
ENDDO
WRITE (61, * ) 'Competed printing matrix ', name

WRITE (61, * )
RETURN
   10 FORMAT(6x,1p10e10.3)
END SUBROUTINE dmat_prt