Consistent formatting of array parameter comments:

    NAME is TYPE array, dimension (dim0,dim1,...)

51 files changed, 93 insertions, 101 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..695cb45d 100644
--- a/SRC/cgesvj.f
+++ b/SRC/cgesvj.f
@@ -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..194904b0 100644
--- a/SRC/claqr5.f
+++ b/SRC/claqr5.f
@@ -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/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/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..e8f6598f 100644
--- a/SRC/dgesvj.f
+++ b/SRC/dgesvj.f
@@ -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..56435f7c 100644
--- a/SRC/dlaqr5.f
+++ b/SRC/dlaqr5.f
@@ -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/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/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/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/sgesvj.f b/SRC/sgesvj.f
index 5e53cea1..371ec8dc 100644
--- a/SRC/sgesvj.f
+++ b/SRC/sgesvj.f
@@ -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..6a007134 100644
--- a/SRC/slaqr5.f
+++ b/SRC/slaqr5.f
@@ -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/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/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..74c99a22 100644
--- a/SRC/zgesvj.f
+++ b/SRC/zgesvj.f
@@ -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..08b1de5f 100644
--- a/SRC/zlaqr5.f
+++ b/SRC/zlaqr5.f
@@ -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