C #MXOTA V1A 17-APR-72. C LAST UPDATE: C SUBROUTINE MXOTA(U,N,M,ND) DIMENSION U(M,M) C C *PURPOSE. C TO OUTPUT A SQUARE MATRIX IN EASILY UNDERSTOOD FORMAT. C C *PARAMETERS. C C U - REAL ARRAY OF DIMENSION "M" X "M" CONTAINING THE MATRIX C OF SIZE "N" X "N" TO BE OUTPUT. "U" IS UNCHANGED. C N - THE LOGICAL SIZE OF THE MATRIX TO BE OUTPUT. C M - THE PHYSICAL SIZE OF THE ARRAY CONTAINING THE MATRIX. C ND - THE FORTRAN CHANNEL NUMBER TO BE USED FOR THE C OUTPUT. C C INPUT PARAMETERS: C U, N, M, ND C OUTPUT PARAMETERS: C C C *METHOD. C THE MATRIX IS OUTPUT BY ROWS WITH 5 ELEMENTS ACROSS A LINE C AND BLANK LINES BETWEEN DIFFERENT ROWS. ROWS AND COLUMNS ARE C NUMBERED AND THE OUTPUT FORMAT FOR THE ELEMENTS IS "G13.6" WHICH C GUARANTEES 6 SIGNIFICANT DIGITS. THE PRINTOUT SPREADS ACROSS 70 C COLUMNS. C C *ACCURACY. C C *RESTRICTIONS. C C *ERROR CONDITIONS. C C *NON STANDARD ROUTINES CALLED. C C *TYPICAL TIMES. C C *ORIGIN. M.R.MANNING. C C *COMMENTS. C C #END. C C #START. LP=ND WRITE(LP,101) (I,I=1,N) 101 FORMAT(///(I9,4I13)) K=10 IF(K.GT.N) K=N L=K+1 DO 20 J=1,N WRITE(LP,102) J,(U(J,I),I=1,K) IF(N.GT.K) WRITE(LP,103) (U(J,I),I=L,N) 20 CONTINUE RETURN 102 FORMAT(/1X,I5,5G13.6) 103 FORMAT(6X,5G13.6) END