1 subroutine printmat(ldim,m,n,iout,key,a)
3 integer ldim,m,n,nlim,iout,i,j,k
5 double precision a(ldim,n)
8 write (iout,1000) (key(k),k=i,nlim)
10 1000 format (/5x,8(6x,a3))
11 1020 format (/80(1h-)/)
13 write (iout,1010) key(j),(a(j,k),k=i,nlim)
16 1010 format (a3,2x,8(f9.4))