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