1      subroutine vib_wmatrix(A,MR,MC,NR,NC)
2* $Id$
3C*****************************************************************
4C                   SUBROUTINE WMATRIX
5C*****************************************************************
6      IMPLICIT none
7c
8#include "stdio.fh"
9c
10      integer mr, mc, nr, nc
11      double precision A(MR,MC)
12c
13      integer i, j, l
14      integer ncl, ncu, nloops
15c
16      NCL=1
17      NCU=10
18      NLOOPS=NC/10+1
19      DO 20 L=1,NLOOPS
20      IF (NCU.GT.NC)NCU=NC
21      WRITE(LUOUT,1100)(J,J=NCL,NCU)
22      DO 10 I=1,NR
23   10 WRITE(LUOUT,1200)I,(A(I,J),J=NCL,NCU)
24      IF(NCU.GE.NC)RETURN
25      NCL=NCU+1
26   20 NCU=NCU+10
27      RETURN
28 1100 FORMAT(/11X,10(I3,9X))
29 1200 FORMAT(1X,I3,10F12.5)
30      END
31