Merge pull request #160 from thielema/cleanup-comments

more cleaned up comments

67 files changed, 140 insertions, 148 deletions

diff --git a/SRC/cgebd2.f b/SRC/cgebd2.f
index e8741202..7ce92a91 100644
--- a/SRC/cgebd2.f
+++ b/SRC/cgebd2.f
@@ -100,7 +100,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is COMPLEX array dimension (min(M,N))
+*>          TAUQ is COMPLEX array, dimension (min(M,N))
 *>          The scalar factors of the elementary reflectors which
 *>          represent the unitary matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/cgebrd.f b/SRC/cgebrd.f
index d01e228a..a15a2aa0 100644
--- a/SRC/cgebrd.f
+++ b/SRC/cgebrd.f
@@ -102,7 +102,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is COMPLEX array dimension (min(M,N))
+*>          TAUQ is COMPLEX array, dimension (min(M,N))
 *>          The scalar factors of the elementary reflectors which
 *>          represent the unitary matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/cgesvj.f b/SRC/cgesvj.f
index 0168c5c3..8cce26b4 100644
--- a/SRC/cgesvj.f
+++ b/SRC/cgesvj.f
@@ -52,7 +52,7 @@
 *
 *> \param[in] JOBA
 *> \verbatim
-*>          JOBA is CHARACTER* 1
+*>          JOBA is CHARACTER*1
 *>          Specifies the structure of A.
 *>          = 'L': The input matrix A is lower triangular;
 *>          = 'U': The input matrix A is upper triangular;
@@ -206,7 +206,7 @@
 *>
 *> \param[in,out] CWORK
 *> \verbatim
-*>          CWORK is COMPLEX array, dimension max(1,LWORK).
+*>          CWORK is COMPLEX array, dimension (max(1,LWORK))
 *>          Used as workspace.
 *>          If on entry LWORK .EQ. -1, then a workspace query is assumed and
 *>          no computation is done; CWORK(1) is set to the minial (and optimal)
@@ -221,7 +221,7 @@
 *>
 *> \param[in,out] RWORK
 *> \verbatim
-*>          RWORK is REAL array, dimension max(6,LRWORK).
+*>          RWORK is REAL array, dimension (max(6,LRWORK))
 *>          On entry,
 *>          If JOBU .EQ. 'C' :
 *>          RWORK(1) = CTOL, where CTOL defines the threshold for convergence.
diff --git a/SRC/cla_heamv.f b/SRC/cla_heamv.f
index 6a3eef1b..2b630631 100644
--- a/SRC/cla_heamv.f
+++ b/SRC/cla_heamv.f
@@ -89,7 +89,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is COMPLEX array of DIMENSION ( LDA, n ).
+*>          A is COMPLEX array, dimension ( LDA, n ).
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/cla_syamv.f b/SRC/cla_syamv.f
index 362d4559..dc8bd5dc 100644
--- a/SRC/cla_syamv.f
+++ b/SRC/cla_syamv.f
@@ -90,7 +90,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is COMPLEX array of DIMENSION ( LDA, n ).
+*>          A is COMPLEX array, dimension ( LDA, n ).
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/clabrd.f b/SRC/clabrd.f
index 87bcb1bc..0972010e 100644
--- a/SRC/clabrd.f
+++ b/SRC/clabrd.f
@@ -111,7 +111,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is COMPLEX array dimension (NB)
+*>          TAUQ is COMPLEX array, dimension (NB)
 *>          The scalar factors of the elementary reflectors which
 *>          represent the unitary matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/claqr1.f b/SRC/claqr1.f
index b129bbab..b448ddb2 100644
--- a/SRC/claqr1.f
+++ b/SRC/claqr1.f
@@ -56,7 +56,7 @@
 *>
 *> \param[in] H
 *> \verbatim
-*>          H is COMPLEX array of dimension (LDH,N)
+*>          H is COMPLEX array, dimension (LDH,N)
 *>              The 2-by-2 or 3-by-3 matrix H in (*).
 *> \endverbatim
 *>
@@ -81,7 +81,7 @@
 *>
 *> \param[out] V
 *> \verbatim
-*>          V is COMPLEX array of dimension N
+*>          V is COMPLEX array, dimension (N)
 *>              A scalar multiple of the first column of the
 *>              matrix K in (*).
 *> \endverbatim
diff --git a/SRC/claqr5.f b/SRC/claqr5.f
index d5c532e4..5b88d7c4 100644
--- a/SRC/claqr5.f
+++ b/SRC/claqr5.f
@@ -47,14 +47,14 @@
 *
 *> \param[in] WANTT
 *> \verbatim
-*>          WANTT is logical scalar
+*>          WANTT is LOGICAL
 *>             WANTT = .true. if the triangular Schur factor
 *>             is being computed.  WANTT is set to .false. otherwise.
 *> \endverbatim
 *>
 *> \param[in] WANTZ
 *> \verbatim
-*>          WANTZ is logical scalar
+*>          WANTZ is LOGICAL
 *>             WANTZ = .true. if the unitary Schur factor is being
 *>             computed.  WANTZ is set to .false. otherwise.
 *> \endverbatim
@@ -107,14 +107,14 @@
 *>
 *> \param[in,out] S
 *> \verbatim
-*>          S is COMPLEX array of size (NSHFTS)
+*>          S is COMPLEX array, dimension (NSHFTS)
 *>             S contains the shifts of origin that define the multi-
 *>             shift QR sweep.  On output S may be reordered.
 *> \endverbatim
 *>
 *> \param[in,out] H
 *> \verbatim
-*>          H is COMPLEX array of size (LDH,N)
+*>          H is COMPLEX array, dimension (LDH,N)
 *>             On input H contains a Hessenberg matrix.  On output a
 *>             multi-shift QR sweep with shifts SR(J)+i*SI(J) is applied
 *>             to the isolated diagonal block in rows and columns KTOP
@@ -142,7 +142,7 @@
 *>
 *> \param[in,out] Z
 *> \verbatim
-*>          Z is COMPLEX array of size (LDZ,IHIZ)
+*>          Z is COMPLEX array, dimension (LDZ,IHIZ)
 *>             If WANTZ = .TRUE., then the QR Sweep unitary
 *>             similarity transformation is accumulated into
 *>             Z(ILOZ:IHIZ,ILOZ:IHIZ) from the right.
@@ -158,7 +158,7 @@
 *>
 *> \param[out] V
 *> \verbatim
-*>          V is COMPLEX array of size (LDV,NSHFTS/2)
+*>          V is COMPLEX array, dimension (LDV,NSHFTS/2)
 *> \endverbatim
 *>
 *> \param[in] LDV
@@ -170,8 +170,7 @@
 *>
 *> \param[out] U
 *> \verbatim
-*>          U is COMPLEX array of size
-*>             (LDU,3*NSHFTS-3)
+*>          U is COMPLEX array, dimension (LDU,3*NSHFTS-3)
 *> \endverbatim
 *>
 *> \param[in] LDU
@@ -190,7 +189,7 @@
 *>
 *> \param[out] WH
 *> \verbatim
-*>          WH is COMPLEX array of size (LDWH,NH)
+*>          WH is COMPLEX array, dimension (LDWH,NH)
 *> \endverbatim
 *>
 *> \param[in] LDWH
@@ -209,8 +208,7 @@
 *>
 *> \param[out] WV
 *> \verbatim
-*>          WV is COMPLEX array of size
-*>             (LDWV,3*NSHFTS-3)
+*>          WV is COMPLEX array, dimension (LDWV,3*NSHFTS-3)
 *> \endverbatim
 *>
 *> \param[in] LDWV
diff --git a/SRC/clarrv.f b/SRC/clarrv.f
index d6503013..72ef0459 100644
--- a/SRC/clarrv.f
+++ b/SRC/clarrv.f
@@ -199,7 +199,7 @@
 *>
 *> \param[out] Z
 *> \verbatim
-*>          Z is array, dimension (LDZ, max(1,M) )
+*>          Z is COMPLEX array, dimension (LDZ, max(1,M) )
 *>          If INFO = 0, the first M columns of Z contain the
 *>          orthonormal eigenvectors of the matrix T
 *>          corresponding to the input eigenvalues, with the i-th
diff --git a/SRC/cstegr.f b/SRC/cstegr.f
index 3209f27e..77a2c210 100644
--- a/SRC/cstegr.f
+++ b/SRC/cstegr.f
@@ -184,7 +184,7 @@
 *>
 *> \param[out] ISUPPZ
 *> \verbatim
-*>          ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*>          ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
 *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th computed eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
diff --git a/SRC/cstemr.f b/SRC/cstemr.f
index f7e0abbe..1e9df5c3 100644
--- a/SRC/cstemr.f
+++ b/SRC/cstemr.f
@@ -239,7 +239,7 @@
 *>
 *> \param[out] ISUPPZ
 *> \verbatim
-*>          ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*>          ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
 *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th computed eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
diff --git a/SRC/ctgex2.f b/SRC/ctgex2.f
index c487e8f1..63ca4781 100644
--- a/SRC/ctgex2.f
+++ b/SRC/ctgex2.f
@@ -76,7 +76,7 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is COMPLEX arrays, dimensions (LDA,N)
+*>          A is COMPLEX array, dimension (LDA,N)
 *>          On entry, the matrix A in the pair (A, B).
 *>          On exit, the updated matrix A.
 *> \endverbatim
@@ -89,7 +89,7 @@
 *>
 *> \param[in,out] B
 *> \verbatim
-*>          B is COMPLEX arrays, dimensions (LDB,N)
+*>          B is COMPLEX array, dimension (LDB,N)
 *>          On entry, the matrix B in the pair (A, B).
 *>          On exit, the updated matrix B.
 *> \endverbatim
diff --git a/SRC/cuncsd.f b/SRC/cuncsd.f
index 829c4bcf..a8637e35 100644
--- a/SRC/cuncsd.f
+++ b/SRC/cuncsd.f
@@ -188,7 +188,7 @@
 *>
 *> \param[out] U1
 *> \verbatim
-*>          U1 is COMPLEX array, dimension (P)
+*>          U1 is COMPLEX array, dimension (P,P)
 *>          If JOBU1 = 'Y', U1 contains the P-by-P unitary matrix U1.
 *> \endverbatim
 *>
@@ -201,7 +201,7 @@
 *>
 *> \param[out] U2
 *> \verbatim
-*>          U2 is COMPLEX array, dimension (M-P)
+*>          U2 is COMPLEX array, dimension (M-P,M-P)
 *>          If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) unitary
 *>          matrix U2.
 *> \endverbatim
@@ -215,7 +215,7 @@
 *>
 *> \param[out] V1T
 *> \verbatim
-*>          V1T is COMPLEX array, dimension (Q)
+*>          V1T is COMPLEX array, dimension (Q,Q)
 *>          If JOBV1T = 'Y', V1T contains the Q-by-Q matrix unitary
 *>          matrix V1**H.
 *> \endverbatim
@@ -229,7 +229,7 @@
 *>
 *> \param[out] V2T
 *> \verbatim
-*>          V2T is COMPLEX array, dimension (M-Q)
+*>          V2T is COMPLEX array, dimension (M-Q,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) unitary
 *>          matrix V2**H.
 *> \endverbatim
diff --git a/SRC/cunm22.f b/SRC/cunm22.f
index 85c2269a..f62939bc 100644
--- a/SRC/cunm22.f
+++ b/SRC/cunm22.f
@@ -52,8 +52,8 @@
 *>  N2-by-N2 upper triangular matrix.
 *> \endverbatim
 *
-*  Arguments
-*  =========
+*  Arguments:
+*  ==========
 *
 *> \param[in] SIDE
 *> \verbatim
diff --git a/SRC/dgebal.f b/SRC/dgebal.f
index 93efd289..0527220a 100644
--- a/SRC/dgebal.f
+++ b/SRC/dgebal.f
@@ -67,7 +67,7 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is DOUBLE array, dimension (LDA,N)
+*>          A is DOUBLE PRECISION array, dimension (LDA,N)
 *>          On entry, the input matrix A.
 *>          On exit,  A is overwritten by the balanced matrix.
 *>          If JOB = 'N', A is not referenced.
@@ -94,7 +94,7 @@
 *>
 *> \param[out] SCALE
 *> \verbatim
-*>          SCALE is DOUBLE array, dimension (N)
+*>          SCALE is DOUBLE PRECISION array, dimension (N)
 *>          Details of the permutations and scaling factors applied to
 *>          A.  If P(j) is the index of the row and column interchanged
 *>          with row and column j and D(j) is the scaling factor
diff --git a/SRC/dgebd2.f b/SRC/dgebd2.f
index bb4035db..c7a7d787 100644
--- a/SRC/dgebd2.f
+++ b/SRC/dgebd2.f
@@ -100,7 +100,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is DOUBLE PRECISION array dimension (min(M,N))
+*>          TAUQ is DOUBLE PRECISION array, dimension (min(M,N))
 *>          The scalar factors of the elementary reflectors which
 *>          represent the orthogonal matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/dgebrd.f b/SRC/dgebrd.f
index 885ad9bb..c8780be3 100644
--- a/SRC/dgebrd.f
+++ b/SRC/dgebrd.f
@@ -101,7 +101,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is DOUBLE PRECISION array dimension (min(M,N))
+*>          TAUQ is DOUBLE PRECISION array, dimension (min(M,N))
 *>          The scalar factors of the elementary reflectors which
 *>          represent the orthogonal matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/dgesvj.f b/SRC/dgesvj.f
index 2b259942..09025c23 100644
--- a/SRC/dgesvj.f
+++ b/SRC/dgesvj.f
@@ -54,7 +54,7 @@
 *
 *> \param[in] JOBA
 *> \verbatim
-*>          JOBA is CHARACTER* 1
+*>          JOBA is CHARACTER*1
 *>          Specifies the structure of A.
 *>          = 'L': The input matrix A is lower triangular;
 *>          = 'U': The input matrix A is upper triangular;
@@ -208,7 +208,7 @@
 *>
 *> \param[in,out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION array, dimension MAX(6,M+N).
+*>          WORK is DOUBLE PRECISION array, dimension (max(6,M+N))
 *>          On entry :
 *>          If JOBU .EQ. 'C' :
 *>          WORK(1) = CTOL, where CTOL defines the threshold for convergence.
diff --git a/SRC/dla_gbamv.f b/SRC/dla_gbamv.f
index 577866cf..707177e9 100644
--- a/SRC/dla_gbamv.f
+++ b/SRC/dla_gbamv.f
@@ -106,7 +106,7 @@
 *>
 *> \param[in] AB
 *> \verbatim
-*>          AB is DOUBLE PRECISION array of DIMENSION ( LDAB, n )
+*>          AB is DOUBLE PRECISION array, dimension ( LDAB, n )
 *>           Before entry, the leading m by n part of the array AB must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/dla_geamv.f b/SRC/dla_geamv.f
index 9a91f6ff..5b6690b8 100644
--- a/SRC/dla_geamv.f
+++ b/SRC/dla_geamv.f
@@ -94,7 +94,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is DOUBLE PRECISION array of DIMENSION ( LDA, n )
+*>          A is DOUBLE PRECISION array, dimension ( LDA, n )
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/dla_syamv.f b/SRC/dla_syamv.f
index 1f948a2d..127371ba 100644
--- a/SRC/dla_syamv.f
+++ b/SRC/dla_syamv.f
@@ -88,7 +88,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is DOUBLE PRECISION array of DIMENSION ( LDA, n ).
+*>          A is DOUBLE PRECISION array, dimension ( LDA, n ).
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/dlabrd.f b/SRC/dlabrd.f
index 36c2e85b..0cad5089 100644
--- a/SRC/dlabrd.f
+++ b/SRC/dlabrd.f
@@ -110,7 +110,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is DOUBLE PRECISION array dimension (NB)
+*>          TAUQ is DOUBLE PRECISION array, dimension (NB)
 *>          The scalar factors of the elementary reflectors which
 *>          represent the orthogonal matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/dlaqr1.f b/SRC/dlaqr1.f
index 25021d5a..c3e6e8ab 100644
--- a/SRC/dlaqr1.f
+++ b/SRC/dlaqr1.f
@@ -61,7 +61,7 @@
 *>
 *> \param[in] H
 *> \verbatim
-*>          H is DOUBLE PRECISION array of dimension (LDH,N)
+*>          H is DOUBLE PRECISION array, dimension (LDH,N)
 *>              The 2-by-2 or 3-by-3 matrix H in (*).
 *> \endverbatim
 *>
@@ -95,7 +95,7 @@
 *>
 *> \param[out] V
 *> \verbatim
-*>          V is DOUBLE PRECISION array of dimension N
+*>          V is DOUBLE PRECISION array, dimension (N)
 *>              A scalar multiple of the first column of the
 *>              matrix K in (*).
 *> \endverbatim
diff --git a/SRC/dlaqr5.f b/SRC/dlaqr5.f
index 9b8763f0..a5a49610 100644
--- a/SRC/dlaqr5.f
+++ b/SRC/dlaqr5.f
@@ -48,14 +48,14 @@
 *
 *> \param[in] WANTT
 *> \verbatim
-*>          WANTT is logical scalar
+*>          WANTT is LOGICAL
 *>             WANTT = .true. if the quasi-triangular Schur factor
 *>             is being computed.  WANTT is set to .false. otherwise.
 *> \endverbatim
 *>
 *> \param[in] WANTZ
 *> \verbatim
-*>          WANTZ is logical scalar
+*>          WANTZ is LOGICAL
 *>             WANTZ = .true. if the orthogonal Schur factor is being
 *>             computed.  WANTZ is set to .false. otherwise.
 *> \endverbatim
@@ -108,12 +108,12 @@
 *>
 *> \param[in,out] SR
 *> \verbatim
-*>          SR is DOUBLE PRECISION array of size (NSHFTS)
+*>          SR is DOUBLE PRECISION array, dimension (NSHFTS)
 *> \endverbatim
 *>
 *> \param[in,out] SI
 *> \verbatim
-*>          SI is DOUBLE PRECISION array of size (NSHFTS)
+*>          SI is DOUBLE PRECISION array, dimension (NSHFTS)
 *>             SR contains the real parts and SI contains the imaginary
 *>             parts of the NSHFTS shifts of origin that define the
 *>             multi-shift QR sweep.  On output SR and SI may be
@@ -122,7 +122,7 @@
 *>
 *> \param[in,out] H
 *> \verbatim
-*>          H is DOUBLE PRECISION array of size (LDH,N)
+*>          H is DOUBLE PRECISION array, dimension (LDH,N)
 *>             On input H contains a Hessenberg matrix.  On output a
 *>             multi-shift QR sweep with shifts SR(J)+i*SI(J) is applied
 *>             to the isolated diagonal block in rows and columns KTOP
@@ -150,7 +150,7 @@
 *>
 *> \param[in,out] Z
 *> \verbatim
-*>          Z is DOUBLE PRECISION array of size (LDZ,IHIZ)
+*>          Z is DOUBLE PRECISION array, dimension (LDZ,IHIZ)
 *>             If WANTZ = .TRUE., then the QR Sweep orthogonal
 *>             similarity transformation is accumulated into
 *>             Z(ILOZ:IHIZ,ILOZ:IHIZ) from the right.
@@ -166,7 +166,7 @@
 *>
 *> \param[out] V
 *> \verbatim
-*>          V is DOUBLE PRECISION array of size (LDV,NSHFTS/2)
+*>          V is DOUBLE PRECISION array, dimension (LDV,NSHFTS/2)
 *> \endverbatim
 *>
 *> \param[in] LDV
@@ -178,8 +178,7 @@
 *>
 *> \param[out] U
 *> \verbatim
-*>          U is DOUBLE PRECISION array of size
-*>             (LDU,3*NSHFTS-3)
+*>          U is DOUBLE PRECISION array, dimension (LDU,3*NSHFTS-3)
 *> \endverbatim
 *>
 *> \param[in] LDU
@@ -198,7 +197,7 @@
 *>
 *> \param[out] WH
 *> \verbatim
-*>          WH is DOUBLE PRECISION array of size (LDWH,NH)
+*>          WH is DOUBLE PRECISION array, dimension (LDWH,NH)
 *> \endverbatim
 *>
 *> \param[in] LDWH
@@ -217,8 +216,7 @@
 *>
 *> \param[out] WV
 *> \verbatim
-*>          WV is DOUBLE PRECISION array of size
-*>             (LDWV,3*NSHFTS-3)
+*>          WV is DOUBLE PRECISION array, dimension (LDWV,3*NSHFTS-3)
 *> \endverbatim
 *>
 *> \param[in] LDWV
diff --git a/SRC/dlasd0.f b/SRC/dlasd0.f
index ca0b3b98..832481a4 100644
--- a/SRC/dlasd0.f
+++ b/SRC/dlasd0.f
@@ -112,14 +112,12 @@
 *>
 *> \param[out] IWORK
 *> \verbatim
-*>          IWORK is INTEGER work array.
-*>         Dimension must be at least (8 * N)
+*>          IWORK is INTEGER array, dimension (8*N)
 *> \endverbatim
 *>
 *> \param[out] WORK
 *> \verbatim
-*>          WORK is DOUBLE PRECISION work array.
-*>         Dimension must be at least (3 * M**2 + 2 * M)
+*>          WORK is DOUBLE PRECISION array, dimension (3*M**2+2*M)
 *> \endverbatim
 *>
 *> \param[out] INFO
diff --git a/SRC/dlasd2.f b/SRC/dlasd2.f
index a7ced418..cbc10c42 100644
--- a/SRC/dlasd2.f
+++ b/SRC/dlasd2.f
@@ -190,7 +190,7 @@
 *>
 *> \param[out] IDXP
 *> \verbatim
-*>          IDXP is INTEGER array dimension(N)
+*>          IDXP is INTEGER array, dimension(N)
 *>         This will contain the permutation used to place deflated
 *>         values of D at the end of the array. On output IDXP(2:K)
 *>         points to the nondeflated D-values and IDXP(K+1:N)
@@ -199,14 +199,14 @@
 *>
 *> \param[out] IDX
 *> \verbatim
-*>          IDX is INTEGER array dimension(N)
+*>          IDX is INTEGER array, dimension(N)
 *>         This will contain the permutation used to sort the contents of
 *>         D into ascending order.
 *> \endverbatim
 *>
 *> \param[out] IDXC
 *> \verbatim
-*>          IDXC is INTEGER array dimension(N)
+*>          IDXC is INTEGER array, dimension(N)
 *>         This will contain the permutation used to arrange the columns
 *>         of the deflated U matrix into three groups:  the first group
 *>         contains non-zero entries only at and above NL, the second
@@ -216,7 +216,7 @@
 *>
 *> \param[in,out] IDXQ
 *> \verbatim
-*>          IDXQ is INTEGER array dimension(N)
+*>          IDXQ is INTEGER array, dimension(N)
 *>         This contains the permutation which separately sorts the two
 *>         sub-problems in D into ascending order.  Note that entries in
 *>         the first hlaf of this permutation must first be moved one
@@ -226,7 +226,7 @@
 *>
 *> \param[out] COLTYP
 *> \verbatim
-*>          COLTYP is INTEGER array dimension(N)
+*>          COLTYP is INTEGER array, dimension(N)
 *>         As workspace, this will contain a label which will indicate
 *>         which of the following types a column in the U2 matrix or a
 *>         row in the VT2 matrix is:
diff --git a/SRC/dlasq5.f b/SRC/dlasq5.f
index 99d4f678..0eb1fd30 100644
--- a/SRC/dlasq5.f
+++ b/SRC/dlasq5.f
@@ -121,7 +121,7 @@
 *>          IEEE is LOGICAL
 *>        Flag for IEEE or non IEEE arithmetic.
 *> \endverbatim
-*
+*>
 *> \param[in] EPS
 *> \verbatim
 *>          EPS is DOUBLE PRECISION
diff --git a/SRC/dorcsd.f b/SRC/dorcsd.f
index 340e16a5..34c1dbba 100644
--- a/SRC/dorcsd.f
+++ b/SRC/dorcsd.f
@@ -186,7 +186,7 @@
 *>
 *> \param[out] U1
 *> \verbatim
-*>          U1 is DOUBLE PRECISION array, dimension (P)
+*>          U1 is DOUBLE PRECISION array, dimension (P,P)
 *>          If JOBU1 = 'Y', U1 contains the P-by-P orthogonal matrix U1.
 *> \endverbatim
 *>
@@ -199,7 +199,7 @@
 *>
 *> \param[out] U2
 *> \verbatim
-*>          U2 is DOUBLE PRECISION array, dimension (M-P)
+*>          U2 is DOUBLE PRECISION array, dimension (M-P,M-P)
 *>          If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) orthogonal
 *>          matrix U2.
 *> \endverbatim
@@ -213,7 +213,7 @@
 *>
 *> \param[out] V1T
 *> \verbatim
-*>          V1T is DOUBLE PRECISION array, dimension (Q)
+*>          V1T is DOUBLE PRECISION array, dimension (Q,Q)
 *>          If JOBV1T = 'Y', V1T contains the Q-by-Q matrix orthogonal
 *>          matrix V1**T.
 *> \endverbatim
@@ -227,7 +227,7 @@
 *>
 *> \param[out] V2T
 *> \verbatim
-*>          V2T is DOUBLE PRECISION array, dimension (M-Q)
+*>          V2T is DOUBLE PRECISION array, dimension (M-Q,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) orthogonal
 *>          matrix V2**T.
 *> \endverbatim
diff --git a/SRC/dorm22.f b/SRC/dorm22.f
index ac79e1e7..3e17d493 100644
--- a/SRC/dorm22.f
+++ b/SRC/dorm22.f
@@ -53,8 +53,8 @@
 *>  N2-by-N2 upper triangular matrix.
 *> \endverbatim
 *
-*  Arguments
-*  =========
+*  Arguments:
+*  ==========
 *
 *> \param[in] SIDE
 *> \verbatim
diff --git a/SRC/dstegr.f b/SRC/dstegr.f
index f3286032..84c08be1 100644
--- a/SRC/dstegr.f
+++ b/SRC/dstegr.f
@@ -184,7 +184,7 @@
 *>
 *> \param[out] ISUPPZ
 *> \verbatim
-*>          ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*>          ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
 *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th computed eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
diff --git a/SRC/dstemr.f b/SRC/dstemr.f
index 924d738d..c50183b9 100644
--- a/SRC/dstemr.f
+++ b/SRC/dstemr.f
@@ -222,7 +222,7 @@
 *>
 *> \param[out] ISUPPZ
 *> \verbatim
-*>          ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*>          ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
 *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th computed eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
diff --git a/SRC/dtgsen.f b/SRC/dtgsen.f
index 9f492393..c61b613c 100644
--- a/SRC/dtgsen.f
+++ b/SRC/dtgsen.f
@@ -222,7 +222,7 @@
 *> \verbatim
 *>          PL is DOUBLE PRECISION
 *> \endverbatim
-
+*>
 *> \param[out] PR
 *> \verbatim
 *>          PR is DOUBLE PRECISION
diff --git a/SRC/ilaclr.f b/SRC/ilaclr.f
index c2e0584b..dcb4fe2d 100644
--- a/SRC/ilaclr.f
+++ b/SRC/ilaclr.f
@@ -53,7 +53,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is array, dimension (LDA,N)
+*>          A is COMPLEX array, dimension (LDA,N)
 *>          The m by n matrix A.
 *> \endverbatim
 *>
diff --git a/SRC/sgebd2.f b/SRC/sgebd2.f
index eb125fba..d00d189d 100644
--- a/SRC/sgebd2.f
+++ b/SRC/sgebd2.f
@@ -100,7 +100,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is REAL array dimension (min(M,N))
+*>          TAUQ is REAL array, dimension (min(M,N))
 *>          The scalar factors of the elementary reflectors which
 *>          represent the orthogonal matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/sgebrd.f b/SRC/sgebrd.f
index 0f38156d..2666c233 100644
--- a/SRC/sgebrd.f
+++ b/SRC/sgebrd.f
@@ -101,7 +101,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is REAL array dimension (min(M,N))
+*>          TAUQ is REAL array, dimension (min(M,N))
 *>          The scalar factors of the elementary reflectors which
 *>          represent the orthogonal matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/sgees.f b/SRC/sgees.f
index 06319bf7..ad65a0c2 100644
--- a/SRC/sgees.f
+++ b/SRC/sgees.f
@@ -80,7 +80,7 @@
 *>
 *> \param[in] SELECT
 *> \verbatim
-*>          SELECT is LOGICAL FUNCTION of two REAL arguments
+*>          SELECT is a LOGICAL FUNCTION of two REAL arguments
 *>          SELECT must be declared EXTERNAL in the calling subroutine.
 *>          If SORT = 'S', SELECT is used to select eigenvalues to sort
 *>          to the top left of the Schur form.
diff --git a/SRC/sgesvj.f b/SRC/sgesvj.f
index 5e53cea1..adce3451 100644
--- a/SRC/sgesvj.f
+++ b/SRC/sgesvj.f
@@ -54,7 +54,7 @@
 *
 *> \param[in] JOBA
 *> \verbatim
-*>          JOBA is CHARACTER* 1
+*>          JOBA is CHARACTER*1
 *>          Specifies the structure of A.
 *>          = 'L': The input matrix A is lower triangular;
 *>          = 'U': The input matrix A is upper triangular;
@@ -208,7 +208,7 @@
 *>
 *> \param[in,out] WORK
 *> \verbatim
-*>          WORK is REAL array, dimension MAX(6,M+N).
+*>          WORK is REAL array, dimension (max(6,M+N))
 *>          On entry,
 *>          If JOBU .EQ. 'C' :
 *>          WORK(1) = CTOL, where CTOL defines the threshold for convergence.
diff --git a/SRC/sla_gbamv.f b/SRC/sla_gbamv.f
index 0798bacd..073f5610 100644
--- a/SRC/sla_gbamv.f
+++ b/SRC/sla_gbamv.f
@@ -106,7 +106,7 @@
 *>
 *> \param[in] AB
 *> \verbatim
-*>          AB is REAL array of DIMENSION ( LDAB, n )
+*>          AB is REAL array, dimension ( LDAB, n )
 *>           Before entry, the leading m by n part of the array AB must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/sla_geamv.f b/SRC/sla_geamv.f
index 35ce8b80..083a823b 100644
--- a/SRC/sla_geamv.f
+++ b/SRC/sla_geamv.f
@@ -94,7 +94,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is REAL array of DIMENSION ( LDA, n )
+*>          A is REAL array, dimension ( LDA, n )
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/sla_syamv.f b/SRC/sla_syamv.f
index 962e17ac..aad2fbf2 100644
--- a/SRC/sla_syamv.f
+++ b/SRC/sla_syamv.f
@@ -88,7 +88,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is REAL array of DIMENSION ( LDA, n ).
+*>          A is REAL array, dimension ( LDA, n ).
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/slabrd.f b/SRC/slabrd.f
index 8073d003..6321c191 100644
--- a/SRC/slabrd.f
+++ b/SRC/slabrd.f
@@ -110,7 +110,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is REAL array dimension (NB)
+*>          TAUQ is REAL array, dimension (NB)
 *>          The scalar factors of the elementary reflectors which
 *>          represent the orthogonal matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/slaqr1.f b/SRC/slaqr1.f
index ca8e28f4..318b976e 100644
--- a/SRC/slaqr1.f
+++ b/SRC/slaqr1.f
@@ -61,7 +61,7 @@
 *>
 *> \param[in] H
 *> \verbatim
-*>          H is REAL array of dimension (LDH,N)
+*>          H is REAL array, dimension (LDH,N)
 *>              The 2-by-2 or 3-by-3 matrix H in (*).
 *> \endverbatim
 *>
@@ -95,7 +95,7 @@
 *>
 *> \param[out] V
 *> \verbatim
-*>          V is REAL array of dimension N
+*>          V is REAL array, dimension (N)
 *>              A scalar multiple of the first column of the
 *>              matrix K in (*).
 *> \endverbatim
diff --git a/SRC/slaqr5.f b/SRC/slaqr5.f
index 5198c1fb..59692bd0 100644
--- a/SRC/slaqr5.f
+++ b/SRC/slaqr5.f
@@ -48,14 +48,14 @@
 *
 *> \param[in] WANTT
 *> \verbatim
-*>          WANTT is logical scalar
+*>          WANTT is LOGICAL
 *>             WANTT = .true. if the quasi-triangular Schur factor
 *>             is being computed.  WANTT is set to .false. otherwise.
 *> \endverbatim
 *>
 *> \param[in] WANTZ
 *> \verbatim
-*>          WANTZ is logical scalar
+*>          WANTZ is LOGICAL
 *>             WANTZ = .true. if the orthogonal Schur factor is being
 *>             computed.  WANTZ is set to .false. otherwise.
 *> \endverbatim
@@ -108,12 +108,12 @@
 *>
 *> \param[in,out] SR
 *> \verbatim
-*>          SR is REAL array of size (NSHFTS)
+*>          SR is REAL array, dimension (NSHFTS)
 *> \endverbatim
 *>
 *> \param[in,out] SI
 *> \verbatim
-*>          SI is REAL array of size (NSHFTS)
+*>          SI is REAL array, dimension (NSHFTS)
 *>             SR contains the real parts and SI contains the imaginary
 *>             parts of the NSHFTS shifts of origin that define the
 *>             multi-shift QR sweep.  On output SR and SI may be
@@ -122,7 +122,7 @@
 *>
 *> \param[in,out] H
 *> \verbatim
-*>          H is REAL array of size (LDH,N)
+*>          H is REAL array, dimension (LDH,N)
 *>             On input H contains a Hessenberg matrix.  On output a
 *>             multi-shift QR sweep with shifts SR(J)+i*SI(J) is applied
 *>             to the isolated diagonal block in rows and columns KTOP
@@ -150,7 +150,7 @@
 *>
 *> \param[in,out] Z
 *> \verbatim
-*>          Z is REAL array of size (LDZ,IHIZ)
+*>          Z is REAL array, dimension (LDZ,IHIZ)
 *>             If WANTZ = .TRUE., then the QR Sweep orthogonal
 *>             similarity transformation is accumulated into
 *>             Z(ILOZ:IHIZ,ILOZ:IHIZ) from the right.
@@ -166,7 +166,7 @@
 *>
 *> \param[out] V
 *> \verbatim
-*>          V is REAL array of size (LDV,NSHFTS/2)
+*>          V is REAL array, dimension (LDV,NSHFTS/2)
 *> \endverbatim
 *>
 *> \param[in] LDV
@@ -178,8 +178,7 @@
 *>
 *> \param[out] U
 *> \verbatim
-*>          U is REAL array of size
-*>             (LDU,3*NSHFTS-3)
+*>          U is REAL array, dimension (LDU,3*NSHFTS-3)
 *> \endverbatim
 *>
 *> \param[in] LDU
@@ -198,7 +197,7 @@
 *>
 *> \param[out] WH
 *> \verbatim
-*>          WH is REAL array of size (LDWH,NH)
+*>          WH is REAL array, dimension (LDWH,NH)
 *> \endverbatim
 *>
 *> \param[in] LDWH
@@ -217,8 +216,7 @@
 *>
 *> \param[out] WV
 *> \verbatim
-*>          WV is REAL array of size
-*>             (LDWV,3*NSHFTS-3)
+*>          WV is REAL array, dimension (LDWV,3*NSHFTS-3)
 *> \endverbatim
 *>
 *> \param[in] LDWV
diff --git a/SRC/sorcsd.f b/SRC/sorcsd.f
index 5621d58e..9694ef54 100644
--- a/SRC/sorcsd.f
+++ b/SRC/sorcsd.f
@@ -186,7 +186,7 @@
 *>
 *> \param[out] U1
 *> \verbatim
-*>          U1 is REAL array, dimension (P)
+*>          U1 is REAL array, dimension (P,P)
 *>          If JOBU1 = 'Y', U1 contains the P-by-P orthogonal matrix U1.
 *> \endverbatim
 *>
@@ -199,7 +199,7 @@
 *>
 *> \param[out] U2
 *> \verbatim
-*>          U2 is REAL array, dimension (M-P)
+*>          U2 is REAL array, dimension (M-P,M-P)
 *>          If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) orthogonal
 *>          matrix U2.
 *> \endverbatim
@@ -213,7 +213,7 @@
 *>
 *> \param[out] V1T
 *> \verbatim
-*>          V1T is REAL array, dimension (Q)
+*>          V1T is REAL array, dimension (Q,Q)
 *>          If JOBV1T = 'Y', V1T contains the Q-by-Q matrix orthogonal
 *>          matrix V1**T.
 *> \endverbatim
@@ -227,7 +227,7 @@
 *>
 *> \param[out] V2T
 *> \verbatim
-*>          V2T is REAL array, dimension (M-Q)
+*>          V2T is REAL array, dimension (M-Q,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) orthogonal
 *>          matrix V2**T.
 *> \endverbatim
diff --git a/SRC/sorm22.f b/SRC/sorm22.f
index fdb5cd8b..06f18d68 100644
--- a/SRC/sorm22.f
+++ b/SRC/sorm22.f
@@ -53,8 +53,8 @@
 *>  N2-by-N2 upper triangular matrix.
 *> \endverbatim
 *
-*  Arguments
-*  =========
+*  Arguments:
+*  ==========
 *
 *> \param[in] SIDE
 *> \verbatim
diff --git a/SRC/ssfrk.f b/SRC/ssfrk.f
index b2cc4fca..a0c0e9ce 100644
--- a/SRC/ssfrk.f
+++ b/SRC/ssfrk.f
@@ -117,7 +117,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is REAL array of DIMENSION (LDA,ka)
+*>          A is REAL array, dimension (LDA,ka)
 *>           where KA
 *>           is K  when TRANS = 'N' or 'n', and is N otherwise. Before
 *>           entry with TRANS = 'N' or 'n', the leading N--by--K part of
diff --git a/SRC/sstegr.f b/SRC/sstegr.f
index 0a28c584..99a4be62 100644
--- a/SRC/sstegr.f
+++ b/SRC/sstegr.f
@@ -184,7 +184,7 @@
 *>
 *> \param[out] ISUPPZ
 *> \verbatim
-*>          ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*>          ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
 *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th computed eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
diff --git a/SRC/sstemr.f b/SRC/sstemr.f
index cff89ef9..917e896c 100644
--- a/SRC/sstemr.f
+++ b/SRC/sstemr.f
@@ -222,7 +222,7 @@
 *>
 *> \param[out] ISUPPZ
 *> \verbatim
-*>          ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*>          ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
 *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th computed eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
diff --git a/SRC/stfsm.f b/SRC/stfsm.f
index b8b81a87..4ac3c3f5 100644
--- a/SRC/stfsm.f
+++ b/SRC/stfsm.f
@@ -159,7 +159,7 @@
 *>
 *> \param[in,out] B
 *> \verbatim
-*>          B is REAL array, DIMENSION (LDB,N)
+*>          B is REAL array, dimension (LDB,N)
 *>           Before entry,  the leading  m by n part of the array  B must
 *>           contain  the  right-hand  side  matrix  B,  and  on exit  is
 *>           overwritten by the solution matrix  X.
diff --git a/SRC/stgex2.f b/SRC/stgex2.f
index fc5f9330..d7795b14 100644
--- a/SRC/stgex2.f
+++ b/SRC/stgex2.f
@@ -77,7 +77,7 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is REAL arrays, dimensions (LDA,N)
+*>          A is REAL array, dimension (LDA,N)
 *>          On entry, the matrix A in the pair (A, B).
 *>          On exit, the updated matrix A.
 *> \endverbatim
@@ -90,7 +90,7 @@
 *>
 *> \param[in,out] B
 *> \verbatim
-*>          B is REAL arrays, dimensions (LDB,N)
+*>          B is REAL array, dimension (LDB,N)
 *>          On entry, the matrix B in the pair (A, B).
 *>          On exit, the updated matrix B.
 *> \endverbatim
diff --git a/SRC/zgebal.f b/SRC/zgebal.f
index 601d5431..94b8aadf 100644
--- a/SRC/zgebal.f
+++ b/SRC/zgebal.f
@@ -83,10 +83,12 @@
 *>
 *> \param[out] ILO
 *> \verbatim
+*>          ILO is INTEGER
 *> \endverbatim
 *>
 *> \param[out] IHI
 *> \verbatim
+*>          IHI is INTEGER
 *>          ILO and IHI are set to INTEGER such that on exit
 *>          A(i,j) = 0 if i > j and j = 1,...,ILO-1 or I = IHI+1,...,N.
 *>          If JOB = 'N' or 'S', ILO = 1 and IHI = N.
diff --git a/SRC/zgebd2.f b/SRC/zgebd2.f
index d5752e04..77779d6c 100644
--- a/SRC/zgebd2.f
+++ b/SRC/zgebd2.f
@@ -100,7 +100,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is COMPLEX*16 array dimension (min(M,N))
+*>          TAUQ is COMPLEX*16 array, dimension (min(M,N))
 *>          The scalar factors of the elementary reflectors which
 *>          represent the unitary matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/zgebrd.f b/SRC/zgebrd.f
index 26879a75..47db32f4 100644
--- a/SRC/zgebrd.f
+++ b/SRC/zgebrd.f
@@ -101,7 +101,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is COMPLEX*16 array dimension (min(M,N))
+*>          TAUQ is COMPLEX*16 array, dimension (min(M,N))
 *>          The scalar factors of the elementary reflectors which
 *>          represent the unitary matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/zgesvj.f b/SRC/zgesvj.f
index 3670ec66..1e1be50c 100644
--- a/SRC/zgesvj.f
+++ b/SRC/zgesvj.f
@@ -52,7 +52,7 @@
 *
 *> \param[in] JOBA
 *> \verbatim
-*>          JOBA is CHARACTER* 1
+*>          JOBA is CHARACTER*1
 *>          Specifies the structure of A.
 *>          = 'L': The input matrix A is lower triangular;
 *>          = 'U': The input matrix A is upper triangular;
@@ -206,7 +206,7 @@
 *>
 *> \param[in,out] CWORK
 *> \verbatim
-*>          CWORK is COMPLEX*16 array, dimension max(1,LWORK).
+*>          CWORK is COMPLEX*16 array, dimension (max(1,LWORK))
 *>          Used as workspace.
 *>          If on entry LWORK .EQ. -1, then a workspace query is assumed and
 *>          no computation is done; CWORK(1) is set to the minial (and optimal)
@@ -221,7 +221,7 @@
 *>
 *> \param[in,out] RWORK
 *> \verbatim
-*>          RWORK is DOUBLE PRECISION array, dimension max(6,LRWORK).
+*>          RWORK is DOUBLE PRECISION array, dimension (max(6,LRWORK))
 *>          On entry,
 *>          If JOBU .EQ. 'C' :
 *>          RWORK(1) = CTOL, where CTOL defines the threshold for convergence.
diff --git a/SRC/zhfrk.f b/SRC/zhfrk.f
index cfc3e111..a4d55078 100644
--- a/SRC/zhfrk.f
+++ b/SRC/zhfrk.f
@@ -117,7 +117,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is COMPLEX*16 array of DIMENSION (LDA,ka)
+*>          A is COMPLEX*16 array, dimension (LDA,ka)
 *>           where KA
 *>           is K  when TRANS = 'N' or 'n', and is N otherwise. Before
 *>           entry with TRANS = 'N' or 'n', the leading N--by--K part of
diff --git a/SRC/zla_gbamv.f b/SRC/zla_gbamv.f
index 9d5291b8..673535fa 100644
--- a/SRC/zla_gbamv.f
+++ b/SRC/zla_gbamv.f
@@ -107,7 +107,7 @@
 *>
 *> \param[in] AB
 *> \verbatim
-*>          AB is COMPLEX*16 array of DIMENSION ( LDAB, n )
+*>          AB is COMPLEX*16 array, dimension ( LDAB, n )
 *>           Before entry, the leading m by n part of the array AB must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
diff --git a/SRC/zla_geamv.f b/SRC/zla_geamv.f
index 8d221691..ffcf6890 100644
--- a/SRC/zla_geamv.f
+++ b/SRC/zla_geamv.f
@@ -96,7 +96,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is COMPLEX*16 array of DIMENSION ( LDA, n )
+*>          A is COMPLEX*16 array, dimension ( LDA, n )
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
@@ -113,7 +113,7 @@
 *>
 *> \param[in] X
 *> \verbatim
-*>          X is COMPLEX*16 array of DIMENSION at least
+*>          X is COMPLEX*16 array, dimension at least
 *>           ( 1 + ( n - 1 )*abs( INCX ) ) when TRANS = 'N' or 'n'
 *>           and at least
 *>           ( 1 + ( m - 1 )*abs( INCX ) ) otherwise.
diff --git a/SRC/zla_heamv.f b/SRC/zla_heamv.f
index 1fa0e0a8..cc6407d6 100644
--- a/SRC/zla_heamv.f
+++ b/SRC/zla_heamv.f
@@ -89,7 +89,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is COMPLEX*16 array, DIMENSION ( LDA, n ).
+*>          A is COMPLEX*16 array, dimension ( LDA, n ).
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
@@ -106,7 +106,7 @@
 *>
 *> \param[in] X
 *> \verbatim
-*>          X is COMPLEX*16 array, DIMENSION at least
+*>          X is COMPLEX*16 array, dimension at least
 *>           ( 1 + ( n - 1 )*abs( INCX ) )
 *>           Before entry, the incremented array X must contain the
 *>           vector x.
diff --git a/SRC/zla_syamv.f b/SRC/zla_syamv.f
index 01854a9b..0f1f4e05 100644
--- a/SRC/zla_syamv.f
+++ b/SRC/zla_syamv.f
@@ -90,7 +90,7 @@
 *>
 *> \param[in] A
 *> \verbatim
-*>          A is COMPLEX*16 array, DIMENSION ( LDA, n ).
+*>          A is COMPLEX*16 array, dimension ( LDA, n ).
 *>           Before entry, the leading m by n part of the array A must
 *>           contain the matrix of coefficients.
 *>           Unchanged on exit.
@@ -107,7 +107,7 @@
 *>
 *> \param[in] X
 *> \verbatim
-*>          X is COMPLEX*16 array, DIMENSION at least
+*>          X is COMPLEX*16 array, dimension at least
 *>           ( 1 + ( n - 1 )*abs( INCX ) )
 *>           Before entry, the incremented array X must contain the
 *>           vector x.
diff --git a/SRC/zlabrd.f b/SRC/zlabrd.f
index 07b5e9fc..812a973b 100644
--- a/SRC/zlabrd.f
+++ b/SRC/zlabrd.f
@@ -111,7 +111,7 @@
 *>
 *> \param[out] TAUQ
 *> \verbatim
-*>          TAUQ is COMPLEX*16 array dimension (NB)
+*>          TAUQ is COMPLEX*16 array, dimension (NB)
 *>          The scalar factors of the elementary reflectors which
 *>          represent the unitary matrix Q. See Further Details.
 *> \endverbatim
diff --git a/SRC/zlaqr1.f b/SRC/zlaqr1.f
index 50b3a682..cc7200eb 100644
--- a/SRC/zlaqr1.f
+++ b/SRC/zlaqr1.f
@@ -56,7 +56,7 @@
 *>
 *> \param[in] H
 *> \verbatim
-*>          H is COMPLEX*16 array of dimension (LDH,N)
+*>          H is COMPLEX*16 array, dimension (LDH,N)
 *>              The 2-by-2 or 3-by-3 matrix H in (*).
 *> \endverbatim
 *>
@@ -81,7 +81,7 @@
 *>
 *> \param[out] V
 *> \verbatim
-*>          V is COMPLEX*16 array of dimension N
+*>          V is COMPLEX*16 array, dimension (N)
 *>              A scalar multiple of the first column of the
 *>              matrix K in (*).
 *> \endverbatim
diff --git a/SRC/zlaqr5.f b/SRC/zlaqr5.f
index 3e463b4a..422d7b21 100644
--- a/SRC/zlaqr5.f
+++ b/SRC/zlaqr5.f
@@ -47,14 +47,14 @@
 *
 *> \param[in] WANTT
 *> \verbatim
-*>          WANTT is logical scalar
+*>          WANTT is LOGICAL
 *>             WANTT = .true. if the triangular Schur factor
 *>             is being computed.  WANTT is set to .false. otherwise.
 *> \endverbatim
 *>
 *> \param[in] WANTZ
 *> \verbatim
-*>          WANTZ is logical scalar
+*>          WANTZ is LOGICAL
 *>             WANTZ = .true. if the unitary Schur factor is being
 *>             computed.  WANTZ is set to .false. otherwise.
 *> \endverbatim
@@ -107,14 +107,14 @@
 *>
 *> \param[in,out] S
 *> \verbatim
-*>          S is COMPLEX*16 array of size (NSHFTS)
+*>          S is COMPLEX*16 array, dimension (NSHFTS)
 *>             S contains the shifts of origin that define the multi-
 *>             shift QR sweep.  On output S may be reordered.
 *> \endverbatim
 *>
 *> \param[in,out] H
 *> \verbatim
-*>          H is COMPLEX*16 array of size (LDH,N)
+*>          H is COMPLEX*16 array, dimension (LDH,N)
 *>             On input H contains a Hessenberg matrix.  On output a
 *>             multi-shift QR sweep with shifts SR(J)+i*SI(J) is applied
 *>             to the isolated diagonal block in rows and columns KTOP
@@ -142,7 +142,7 @@
 *>
 *> \param[in,out] Z
 *> \verbatim
-*>          Z is COMPLEX*16 array of size (LDZ,IHIZ)
+*>          Z is COMPLEX*16 array, dimension (LDZ,IHIZ)
 *>             If WANTZ = .TRUE., then the QR Sweep unitary
 *>             similarity transformation is accumulated into
 *>             Z(ILOZ:IHIZ,ILOZ:IHIZ) from the right.
@@ -158,7 +158,7 @@
 *>
 *> \param[out] V
 *> \verbatim
-*>          V is COMPLEX*16 array of size (LDV,NSHFTS/2)
+*>          V is COMPLEX*16 array, dimension (LDV,NSHFTS/2)
 *> \endverbatim
 *>
 *> \param[in] LDV
@@ -170,8 +170,7 @@
 *>
 *> \param[out] U
 *> \verbatim
-*>          U is COMPLEX*16 array of size
-*>             (LDU,3*NSHFTS-3)
+*>          U is COMPLEX*16 array, dimension (LDU,3*NSHFTS-3)
 *> \endverbatim
 *>
 *> \param[in] LDU
@@ -190,7 +189,7 @@
 *>
 *> \param[out] WH
 *> \verbatim
-*>          WH is COMPLEX*16 array of size (LDWH,NH)
+*>          WH is COMPLEX*16 array, dimension (LDWH,NH)
 *> \endverbatim
 *>
 *> \param[in] LDWH
@@ -209,8 +208,7 @@
 *>
 *> \param[out] WV
 *> \verbatim
-*>          WV is COMPLEX*16 array of size
-*>             (LDWV,3*NSHFTS-3)
+*>          WV is COMPLEX*16 array, dimension (LDWV,3*NSHFTS-3)
 *> \endverbatim
 *>
 *> \param[in] LDWV
diff --git a/SRC/zstegr.f b/SRC/zstegr.f
index 71f304c1..b8d20729 100644
--- a/SRC/zstegr.f
+++ b/SRC/zstegr.f
@@ -184,7 +184,7 @@
 *>
 *> \param[out] ISUPPZ
 *> \verbatim
-*>          ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*>          ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
 *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th computed eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
diff --git a/SRC/zstemr.f b/SRC/zstemr.f
index 681c87ad..7fa8f3cb 100644
--- a/SRC/zstemr.f
+++ b/SRC/zstemr.f
@@ -239,7 +239,7 @@
 *>
 *> \param[out] ISUPPZ
 *> \verbatim
-*>          ISUPPZ is INTEGER ARRAY, dimension ( 2*max(1,M) )
+*>          ISUPPZ is INTEGER array, dimension ( 2*max(1,M) )
 *>          The support of the eigenvectors in Z, i.e., the indices
 *>          indicating the nonzero elements in Z. The i-th computed eigenvector
 *>          is nonzero only in elements ISUPPZ( 2*i-1 ) through
diff --git a/SRC/ztgex2.f b/SRC/ztgex2.f
index 26b83dec..3f484285 100644
--- a/SRC/ztgex2.f
+++ b/SRC/ztgex2.f
@@ -76,7 +76,7 @@
 *>
 *> \param[in,out] A
 *> \verbatim
-*>          A is COMPLEX*16 arrays, dimensions (LDA,N)
+*>          A is COMPLEX*16 array, dimensions (LDA,N)
 *>          On entry, the matrix A in the pair (A, B).
 *>          On exit, the updated matrix A.
 *> \endverbatim
@@ -89,7 +89,7 @@
 *>
 *> \param[in,out] B
 *> \verbatim
-*>          B is COMPLEX*16 arrays, dimensions (LDB,N)
+*>          B is COMPLEX*16 array, dimensions (LDB,N)
 *>          On entry, the matrix B in the pair (A, B).
 *>          On exit, the updated matrix B.
 *> \endverbatim
diff --git a/SRC/zuncsd.f b/SRC/zuncsd.f
index 77a83c09..18dce987 100644
--- a/SRC/zuncsd.f
+++ b/SRC/zuncsd.f
@@ -188,7 +188,7 @@
 *>
 *> \param[out] U1
 *> \verbatim
-*>          U1 is COMPLEX*16 array, dimension (P)
+*>          U1 is COMPLEX*16 array, dimension (P,P)
 *>          If JOBU1 = 'Y', U1 contains the P-by-P unitary matrix U1.
 *> \endverbatim
 *>
@@ -201,7 +201,7 @@
 *>
 *> \param[out] U2
 *> \verbatim
-*>          U2 is COMPLEX*16 array, dimension (M-P)
+*>          U2 is COMPLEX*16 array, dimension (M-P,M-P)
 *>          If JOBU2 = 'Y', U2 contains the (M-P)-by-(M-P) unitary
 *>          matrix U2.
 *> \endverbatim
@@ -215,7 +215,7 @@
 *>
 *> \param[out] V1T
 *> \verbatim
-*>          V1T is COMPLEX*16 array, dimension (Q)
+*>          V1T is COMPLEX*16 array, dimension (Q,Q)
 *>          If JOBV1T = 'Y', V1T contains the Q-by-Q matrix unitary
 *>          matrix V1**H.
 *> \endverbatim
@@ -229,7 +229,7 @@
 *>
 *> \param[out] V2T
 *> \verbatim
-*>          V2T is COMPLEX*16 array, dimension (M-Q)
+*>          V2T is COMPLEX*16 array, dimension (M-Q,M-Q)
 *>          If JOBV2T = 'Y', V2T contains the (M-Q)-by-(M-Q) unitary
 *>          matrix V2**H.
 *> \endverbatim
diff --git a/SRC/zunm22.f b/SRC/zunm22.f
index 468d7d8c..01f8b73c 100644
--- a/SRC/zunm22.f
+++ b/SRC/zunm22.f
@@ -52,8 +52,8 @@
 *>  N2-by-N2 upper triangular matrix.
 *> \endverbatim
 *
-*  Arguments
-*  =========
+*  Arguments:
+*  ==========
 *
 *> \param[in] SIDE
 *> \verbatim