@@ -327,13 +327,19 @@ Subroutine
 
 `array`: Shall be a 2D array of `integer`, `real`, or `complex` type. It is an `intent(in)` argument.
 
-`unit`: Shall be an integer containing the output unit. It is an `intent(in)` argument. The default is `6` (standard output).
+`unit`: Shall be an integer containing the output unit. It is an `intent(in)` argument. The default is the intrinsic `output_unit` provided by `iso_fortran_env`.
 
 `fmt`: Shall be a character string containing the format for printing the array. It is an `intent(in)` argument. The default is based on [the Formatting constants](#formatting-constants).
 
-`delimiter`: Shall be a character string of length 1 containing the delimiter between array elements. It is an `intent(in)` argument. The default is a `" "` (space).
+`delimiter`: Shall be a character string containing the delimiter between array elements. It is an `intent(in)` argument. The default is a `" "` (space).
 
-`brief`: Shall be a logical flag. If `.true.`, the array is printed in a brief format. The default is `.true.`.
+`brief`: Shall be a logical flag. The default is `.true.`. If `.true.`, the array is printed in a shortened/abridged version
+that shows only the representative portions of large arrays, which is useful for gaining a glimpse of large arrays. Specifically:
++ For arrays with more than 5 rows or columns, it will display:
+- First 3 rows and columns;
+- Last row and column;
+- Eilipsis (`...`) to indicate omitted elements.
++ For arrays with 5 rows or columns or less, it will display the entire array.
 
 ### Example
 

@@ -11,6 +11,6 @@ program example_io_print_array
 call print_array(array, unit=6, fmt='(i3)', delimiter='|', brief=.true.)
 
 print "(a)", "=== print_array 2 ==="
-call print_array(array(:1, :))
+call print_array(array(:1, :), delimiter=", ")
 
 end program example_io_print_array

@@ -112,7 +112,7 @@ module stdlib_io
 ${t1}$, intent(in) :: array(:, :)
 integer, intent(in), optional :: unit
 character(len=*), intent(in), optional :: fmt
-character(len=1), intent(in), optional :: delimiter
+character(len=*), intent(in), optional :: delimiter
 logical, intent(in), optional :: brief
 end subroutine print_array_${t1[0]}$${k1}$
 #:endfor

@@ -12,13 +12,13 @@ contains
 ${t1}$, intent(in) :: array(:, :)
 integer, intent(in), optional :: unit
 character(len=*), intent(in), optional :: fmt
-character(len=1), intent(in), optional :: delimiter
+character(len=*), intent(in), optional :: delimiter
 logical, intent(in), optional :: brief
 
 integer :: i, j, unit_, shape_(2)
 character(len=:), allocatable :: fmt_
-character(len=1) :: delimiter_
-character(len=3) :: delim_str
+character(len=:), allocatable :: delimiter_
+character(len=:), allocatable :: delim_str
 logical :: brief_
 
 shape_ = shape(array)
@@ -48,15 +48,15 @@ contains
 write (unit_, fmt='(a)', advance='no') delimiter_//"..."//delimiter_
 write (unit_, fmt=fmt_) array(i, shape_(2))
 end do
-write (unit_, fmt='(a)') ":"
+write (unit_, fmt='(a)') "..."
 write (unit_, fmt=fmt_, advance='no') array(shape_(1), :3)
 write (unit_, fmt='(a)', advance='no') delimiter_//"..."//delimiter_
 write (unit_, fmt=fmt_) array(shape_(1), shape_(2))
 else
 do i = 1, 3
 write (unit_, fmt=fmt_) array(i, :)
 end do
-write (unit_, fmt='(a)') ":"
+write (unit_, fmt='(a)') "..."
 write (unit_, fmt=fmt_) array(shape_(1), :)
 
 end if

@@ -41,7 +41,7 @@ subroutine test_print_rdp(error)
 line(1) = " 1.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 ... 0.0000000000000000E+000"
 line(2) = " 0.0000000000000000E+000 1.0000000000000000E+000 0.0000000000000000E+000 ... 0.0000000000000000E+000"
 line(3) = " 0.0000000000000000E+000 0.0000000000000000E+000 1.0000000000000000E+000 ... 0.0000000000000000E+000"
-line(4) = ":"
+line(4) = "..."
 line(5) = " 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 ... 1.0000000000000000E+000"
 call print_array(a, fh)
 
@@ -87,7 +87,7 @@ subroutine test_print_rsp(error)
 line(1) = " 1.00000000E+00 0.00000000E+00 0.00000000E+00 ... 0.00000000E+00"
 line(2) = " 0.00000000E+00 1.00000000E+00 0.00000000E+00 ... 0.00000000E+00"
 line(3) = " 0.00000000E+00 0.00000000E+00 1.00000000E+00 ... 0.00000000E+00"
-line(4) = ":"
+line(4) = "..."
 line(5) = " 0.00000000E+00 0.00000000E+00 0.00000000E+00 ... 1.00000000E+00"
 call print_array(a, fh)
 
@@ -129,12 +129,12 @@ subroutine test_print_i1(error)
 a = eye(10)
 open (newunit=fh, status='scratch')
 
-line(1) = "1 0 0 ... 0"
-line(2) = "0 1 0 ... 0"
-line(3) = "0 0 1 ... 0"
-line(4) = ":"
-line(5) = "0 0 0 ... 1"
-call print_array(a, fh)
+line(1) = "1, 0, 0, ..., 0"
+line(2) = "0, 1, 0, ..., 0"
+line(3) = "0, 0, 1, ..., 0"
+line(4) = "..."
+line(5) = "0, 0, 0, ..., 1"
+call print_array(a, fh, delimiter=", ")
 
 rewind (fh)
 do i = 1, 5
@@ -144,11 +144,11 @@ subroutine test_print_i1(error)
 end do
 
 rewind (fh)
-line(1) = "01;00;00;00;00"
-line(2) = "00;01;00;00;00"
-line(3) = "00;00;01;00;00"
-line(4) = "00;00;00;01;00"
-line(5) = "00;00;00;00;01"
+line(1) = "01;00;00;00;00"
+line(2) = "00;01;00;00;00"
+line(3) = "00;00;01;00;00"
+line(4) = "00;00;00;01;00"
+line(5) = "00;00;00;00;01"
 line(6:) = "00;00;00;00;00"
 call print_array(a(:, :5), fh, fmt="(i0.2)", brief=.false., delimiter=";")
 rewind (fh)
@@ -176,7 +176,7 @@ subroutine test_print_i2(error)
 line(1) = "1 0 0 ... 0"
 line(2) = "0 1 0 ... 0"
 line(3) = "0 0 1 ... 0"
-line(4) = ":"
+line(4) = "..."
 line(5) = "0 0 0 ... 1"
 
 call print_array(a, fh)
@@ -188,13 +188,13 @@ subroutine test_print_i2(error)
 end do
 
 rewind (fh)
-line(1) = "01;00;00;00;00"
-line(2) = "00;01;00;00;00"
-line(3) = "00;00;01;00;00"
-line(4) = "00;00;00;01;00"
-line(5) = "00;00;00;00;01"
-line(6:) = "00;00;00;00;00"
-call print_array(a(:, :5), fh, fmt="(i0.2)", brief=.false., delimiter=";")
+line(1) = "01; 00; 00; 00; 00"
+line(2) = "00; 01; 00; 00; 00"
+line(3) = "00; 00; 01; 00; 00"
+line(4) = "00; 00; 00; 01; 00"
+line(5) = "00; 00; 00; 00; 01"
+line(6:) = "00; 00; 00; 00; 00"
+call print_array(a(:, :5), fh, fmt="(i0.2)", brief=.false., delimiter="; ")
 rewind (fh)
 do i = 1, 10
 call get_line(fh, buffer)
@@ -223,7 +223,7 @@ subroutine test_print_cdp(error)
 &0.0000000000000000E+000 0.0000000000000000E+000 ... 0.0000000000000000E+000 0.0000000000000000E+000"
 line(3) = " 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 &
 &1.0000000000000000E+000 0.0000000000000000E+000 ... 0.0000000000000000E+000 0.0000000000000000E+000"
-line(4) = ":"
+line(4) = "..."
 line(5) = " 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 0.0000000000000000E+000 &
 &0.0000000000000000E+000 0.0000000000000000E+000 ... 1.0000000000000000E+000 0.0000000000000000E+000"
 
@@ -236,11 +236,11 @@ subroutine test_print_cdp(error)
 end do
 
 rewind (fh)
-line(1) = "1.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00"
-line(2) = "0.00,0.00|1.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00"
-line(3) = "0.00,0.00|0.00,0.00|1.00,0.00|0.00,0.00|0.00,0.00"
-line(4) = "0.00,0.00|0.00,0.00|0.00,0.00|1.00,0.00|0.00,0.00"
-line(5) = "0.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00|1.00,0.00"
+line(1) = "1.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00"
+line(2) = "0.00,0.00|1.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00"
+line(3) = "0.00,0.00|0.00,0.00|1.00,0.00|0.00,0.00|0.00,0.00"
+line(4) = "0.00,0.00|0.00,0.00|0.00,0.00|1.00,0.00|0.00,0.00"
+line(5) = "0.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00|1.00,0.00"
 line(6:) = "0.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00|0.00,0.00"
 call print_array(a(:, :5), fh, fmt="(f4.2,"","",f4.2)", brief=.false., delimiter="|")
 rewind (fh)
@@ -271,7 +271,7 @@ subroutine test_print_csp(error)
 &0.00000000E+00 0.00000000E+00"
 line(3) = " 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 1.00000000E+00 0.00000000E+00 ... &
 &0.00000000E+00 0.00000000E+00"
-line(4) = ":"
+line(4) = "..."
 line(5) = " 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 0.00000000E+00 ... &
 &1.00000000E+00 0.00000000E+00"
 
@@ -284,11 +284,11 @@ subroutine test_print_csp(error)
 end do
 
 rewind (fh)
-line(1) = " 1.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00"
-line(2) = " 0.00, 0.00; 1.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00"
-line(3) = " 0.00, 0.00; 0.00, 0.00; 1.00, 0.00; 0.00, 0.00; 0.00, 0.00"
-line(4) = " 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 1.00, 0.00; 0.00, 0.00"
-line(5) = " 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 1.00, 0.00"
+line(1) = " 1.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00"
+line(2) = " 0.00, 0.00; 1.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00"
+line(3) = " 0.00, 0.00; 0.00, 0.00; 1.00, 0.00; 0.00, 0.00; 0.00, 0.00"
+line(4) = " 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 1.00, 0.00; 0.00, 0.00"
+line(5) = " 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 1.00, 0.00"
 line(6:) = " 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00; 0.00, 0.00"
 call print_array(a(:, :5), fh, fmt="(1x,f4.2,"","",1x,f4.2)", brief=.false., delimiter=";")
 rewind (fh)