dcstringputline.f90

Go to the documentation of this file.

! -*- mode: f90; coding: utf-8 -*-
!-----------------------------------------------------------------------
! Copyright (c) 2000-2026 Gtool Development Group. All rights reserved.
!-----------------------------------------------------------------------
!>
!> @author    Youhei SASAKI, Yasuhiro MORIKAWA, Eizi TOYODA
!> @copyright Copyright (C) GFD Dennou Club, 2000-2026. All rights reserved. <br/>
!>            License is BSD-2-Clause. see [COPYRIGHT](@ref COPYRIGHT) in detail
!>
!> **Important**
!>
!> This file is generated from ../../../../src/dc_utils/dcstringputline.erb by ERB included Ruby 3.3.8.
!> Please do not edit this file directly.
!>
!> @en
!> @brief Print array summary (PutLine subroutines)
!> @details
!> This file contains PutLine subroutines for printing array summaries.
!> Each subroutine prints the array size, maximum value, minimum value,
!> and average value. Optionally, standard deviation can also be printed.
!>
!> Generated subroutines:
!> - PutLineInt1 to PutLineInt7: for integer arrays of 1 to 7 dimensions
!> - PutLineReal1 to PutLineReal7: for single precision real arrays of 1 to 7 dimensions
!> - PutLineDouble1 to PutLineDouble7: for double precision real arrays of 1 to 7 dimensions
!>
!> @enden
!>
!> @ja
!> @brief 配列の要約印字（PutLine サブルーチン）
!> @details
!> 配列の要約を印字する PutLine サブルーチン群を含むファイルです。
!> 各サブルーチンは配列のサイズ、最大値、最小値、平均値を印字します。
!> オプションで標準偏差も印字できます。
!>
!> 生成されるサブルーチン:
!> - PutLineInt1 〜 PutLineInt7: 1〜7次元の整数配列用
!> - PutLineReal1 〜 PutLineReal7: 1〜7次元の単精度実数配列用
!> - PutLineDouble1 〜 PutLineDouble7: 1〜7次元の倍精度実数配列用
!>
!> @endja
!>
!> @en
!> @brief Print summary of 1-dimensional int array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 1-dimensional int array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 1次元整数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する1次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlineint1( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  integer, intent(in):: array(:)
  integer, intent(in), optional:: lbounds(1)
  integer, intent(in), optional:: ubounds(1)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(1), ubound_nums(1)
  character(STRING):: size_str, sd_str
  integer:: max_value, min_value
  integer, allocatable:: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = array
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<INT-ARRAY:: @size=%c, @max=%d, @min=%d, @avg=%f%c>', &
    & i = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlineint1
!> @en
!> @brief Print summary of 2-dimensional int array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 2-dimensional int array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 2次元整数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する2次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlineint2( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  integer, intent(in):: array(:,:)
  integer, intent(in), optional:: lbounds(2)
  integer, intent(in), optional:: ubounds(2)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(2), ubound_nums(2)
  character(STRING):: size_str, sd_str
  integer:: max_value, min_value
  integer, allocatable:: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<INT-ARRAY:: @size=%c, @max=%d, @min=%d, @avg=%f%c>', &
    & i = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlineint2
!> @en
!> @brief Print summary of 3-dimensional int array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 3-dimensional int array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 3次元整数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する3次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlineint3( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  integer, intent(in):: array(:,:,:)
  integer, intent(in), optional:: lbounds(3)
  integer, intent(in), optional:: ubounds(3)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(3), ubound_nums(3)
  character(STRING):: size_str, sd_str
  integer:: max_value, min_value
  integer, allocatable:: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<INT-ARRAY:: @size=%c, @max=%d, @min=%d, @avg=%f%c>', &
    & i = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlineint3
!> @en
!> @brief Print summary of 4-dimensional int array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 4-dimensional int array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 4次元整数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する4次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlineint4( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  integer, intent(in):: array(:,:,:,:)
  integer, intent(in), optional:: lbounds(4)
  integer, intent(in), optional:: ubounds(4)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(4), ubound_nums(4)
  character(STRING):: size_str, sd_str
  integer:: max_value, min_value
  integer, allocatable:: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<INT-ARRAY:: @size=%c, @max=%d, @min=%d, @avg=%f%c>', &
    & i = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlineint4
!> @en
!> @brief Print summary of 5-dimensional int array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 5-dimensional int array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 5次元整数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する5次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlineint5( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  integer, intent(in):: array(:,:,:,:,:)
  integer, intent(in), optional:: lbounds(5)
  integer, intent(in), optional:: ubounds(5)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(5), ubound_nums(5)
  character(STRING):: size_str, sd_str
  integer:: max_value, min_value
  integer, allocatable:: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<INT-ARRAY:: @size=%c, @max=%d, @min=%d, @avg=%f%c>', &
    & i = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlineint5
!> @en
!> @brief Print summary of 6-dimensional int array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 6-dimensional int array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 6次元整数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する6次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlineint6( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  integer, intent(in):: array(:,:,:,:,:,:)
  integer, intent(in), optional:: lbounds(6)
  integer, intent(in), optional:: ubounds(6)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(6), ubound_nums(6)
  character(STRING):: size_str, sd_str
  integer:: max_value, min_value
  integer, allocatable:: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
    lbound_nums(6) = lbound( array, 6 )
    ubound_nums(6) = ubound( array, 6 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(6)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(6)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<INT-ARRAY:: @size=%c, @max=%d, @min=%d, @avg=%f%c>', &
    & i = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlineint6
!> @en
!> @brief Print summary of 7-dimensional int array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 7-dimensional int array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 7次元整数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する7次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlineint7( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  integer, intent(in):: array(:,:,:,:,:,:,:)
  integer, intent(in), optional:: lbounds(7)
  integer, intent(in), optional:: ubounds(7)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(7), ubound_nums(7)
  character(STRING):: size_str, sd_str
  integer:: max_value, min_value
  integer, allocatable:: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
    lbound_nums(6) = lbound( array, 6 )
    ubound_nums(6) = ubound( array, 6 )
    lbound_nums(7) = lbound( array, 7 )
    ubound_nums(7) = ubound( array, 7 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(6)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(6)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(7)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(7)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<INT-ARRAY:: @size=%c, @max=%d, @min=%d, @avg=%f%c>', &
    & i = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlineint7
!> @en
!> @brief Print summary of 1-dimensional real array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 1-dimensional real array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 1次元単精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する1次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinereal1( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout,sp
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(SP), intent(in):: array(:)
  integer, intent(in), optional:: lbounds(1)
  integer, intent(in), optional:: ubounds(1)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(1), ubound_nums(1)
  character(STRING):: size_str, sd_str
  real(SP) :: max_value, min_value
  real(SP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = array
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<SP-ARRAY:: @size=%c, @max=%r, @min=%r, @avg=%f%c>', &
    & r = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinereal1
!> @en
!> @brief Print summary of 2-dimensional real array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 2-dimensional real array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 2次元単精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する2次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinereal2( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout,sp
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(SP), intent(in):: array(:,:)
  integer, intent(in), optional:: lbounds(2)
  integer, intent(in), optional:: ubounds(2)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(2), ubound_nums(2)
  character(STRING):: size_str, sd_str
  real(SP) :: max_value, min_value
  real(SP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<SP-ARRAY:: @size=%c, @max=%r, @min=%r, @avg=%f%c>', &
    & r = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinereal2
!> @en
!> @brief Print summary of 3-dimensional real array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 3-dimensional real array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 3次元単精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する3次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinereal3( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout,sp
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(SP), intent(in):: array(:,:,:)
  integer, intent(in), optional:: lbounds(3)
  integer, intent(in), optional:: ubounds(3)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(3), ubound_nums(3)
  character(STRING):: size_str, sd_str
  real(SP) :: max_value, min_value
  real(SP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<SP-ARRAY:: @size=%c, @max=%r, @min=%r, @avg=%f%c>', &
    & r = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinereal3
!> @en
!> @brief Print summary of 4-dimensional real array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 4-dimensional real array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 4次元単精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する4次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinereal4( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout,sp
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(SP), intent(in):: array(:,:,:,:)
  integer, intent(in), optional:: lbounds(4)
  integer, intent(in), optional:: ubounds(4)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(4), ubound_nums(4)
  character(STRING):: size_str, sd_str
  real(SP) :: max_value, min_value
  real(SP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<SP-ARRAY:: @size=%c, @max=%r, @min=%r, @avg=%f%c>', &
    & r = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinereal4
!> @en
!> @brief Print summary of 5-dimensional real array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 5-dimensional real array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 5次元単精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する5次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinereal5( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout,sp
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(SP), intent(in):: array(:,:,:,:,:)
  integer, intent(in), optional:: lbounds(5)
  integer, intent(in), optional:: ubounds(5)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(5), ubound_nums(5)
  character(STRING):: size_str, sd_str
  real(SP) :: max_value, min_value
  real(SP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<SP-ARRAY:: @size=%c, @max=%r, @min=%r, @avg=%f%c>', &
    & r = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinereal5
!> @en
!> @brief Print summary of 6-dimensional real array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 6-dimensional real array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 6次元単精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する6次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinereal6( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout,sp
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(SP), intent(in):: array(:,:,:,:,:,:)
  integer, intent(in), optional:: lbounds(6)
  integer, intent(in), optional:: ubounds(6)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(6), ubound_nums(6)
  character(STRING):: size_str, sd_str
  real(SP) :: max_value, min_value
  real(SP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
    lbound_nums(6) = lbound( array, 6 )
    ubound_nums(6) = ubound( array, 6 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(6)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(6)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<SP-ARRAY:: @size=%c, @max=%r, @min=%r, @avg=%f%c>', &
    & r = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinereal6
!> @en
!> @brief Print summary of 7-dimensional real array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 7-dimensional real array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 7次元単精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する7次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinereal7( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout,sp
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(SP), intent(in):: array(:,:,:,:,:,:,:)
  integer, intent(in), optional:: lbounds(7)
  integer, intent(in), optional:: ubounds(7)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(7), ubound_nums(7)
  character(STRING):: size_str, sd_str
  real(SP) :: max_value, min_value
  real(SP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
    lbound_nums(6) = lbound( array, 6 )
    ubound_nums(6) = ubound( array, 6 )
    lbound_nums(7) = lbound( array, 7 )
    ubound_nums(7) = ubound( array, 7 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(6)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(6)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(7)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(7)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<SP-ARRAY:: @size=%c, @max=%r, @min=%r, @avg=%f%c>', &
    & r = (/max_value, min_value/), d = (/avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinereal7
!> @en
!> @brief Print summary of 1-dimensional double array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 1-dimensional double array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 1次元倍精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する1次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinedouble1( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(DP), intent(in):: array(:)
  integer, intent(in), optional:: lbounds(1)
  integer, intent(in), optional:: ubounds(1)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(1), ubound_nums(1)
  character(STRING):: size_str, sd_str
  real(DP) :: max_value, min_value
  real(DP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = array
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<DP-ARRAY:: @size=%c, @max=%f, @min=%f, @avg=%f%c>', &
    & d = (/max_value, min_value, avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinedouble1
!> @en
!> @brief Print summary of 2-dimensional double array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 2-dimensional double array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 2次元倍精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する2次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinedouble2( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(DP), intent(in):: array(:,:)
  integer, intent(in), optional:: lbounds(2)
  integer, intent(in), optional:: ubounds(2)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(2), ubound_nums(2)
  character(STRING):: size_str, sd_str
  real(DP) :: max_value, min_value
  real(DP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<DP-ARRAY:: @size=%c, @max=%f, @min=%f, @avg=%f%c>', &
    & d = (/max_value, min_value, avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinedouble2
!> @en
!> @brief Print summary of 3-dimensional double array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 3-dimensional double array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 3次元倍精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する3次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinedouble3( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(DP), intent(in):: array(:,:,:)
  integer, intent(in), optional:: lbounds(3)
  integer, intent(in), optional:: ubounds(3)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(3), ubound_nums(3)
  character(STRING):: size_str, sd_str
  real(DP) :: max_value, min_value
  real(DP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<DP-ARRAY:: @size=%c, @max=%f, @min=%f, @avg=%f%c>', &
    & d = (/max_value, min_value, avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinedouble3
!> @en
!> @brief Print summary of 4-dimensional double array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 4-dimensional double array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 4次元倍精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する4次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinedouble4( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(DP), intent(in):: array(:,:,:,:)
  integer, intent(in), optional:: lbounds(4)
  integer, intent(in), optional:: ubounds(4)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(4), ubound_nums(4)
  character(STRING):: size_str, sd_str
  real(DP) :: max_value, min_value
  real(DP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<DP-ARRAY:: @size=%c, @max=%f, @min=%f, @avg=%f%c>', &
    & d = (/max_value, min_value, avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinedouble4
!> @en
!> @brief Print summary of 5-dimensional double array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 5-dimensional double array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 5次元倍精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する5次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinedouble5( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(DP), intent(in):: array(:,:,:,:,:)
  integer, intent(in), optional:: lbounds(5)
  integer, intent(in), optional:: ubounds(5)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(5), ubound_nums(5)
  character(STRING):: size_str, sd_str
  real(DP) :: max_value, min_value
  real(DP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<DP-ARRAY:: @size=%c, @max=%f, @min=%f, @avg=%f%c>', &
    & d = (/max_value, min_value, avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinedouble5
!> @en
!> @brief Print summary of 6-dimensional double array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 6-dimensional double array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 6次元倍精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する6次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinedouble6( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(DP), intent(in):: array(:,:,:,:,:,:)
  integer, intent(in), optional:: lbounds(6)
  integer, intent(in), optional:: ubounds(6)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(6), ubound_nums(6)
  character(STRING):: size_str, sd_str
  real(DP) :: max_value, min_value
  real(DP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
    lbound_nums(6) = lbound( array, 6 )
    ubound_nums(6) = ubound( array, 6 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(6)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(6)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<DP-ARRAY:: @size=%c, @max=%f, @min=%f, @avg=%f%c>', &
    & d = (/max_value, min_value, avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinedouble6
!> @en
!> @brief Print summary of 7-dimensional double array
!> @details
!> Prints array size, maximum, minimum, and average values.
!> If sd is .true., also prints standard deviation.
!> @param[in] array 7-dimensional double array to summarize
!> @param[in] lbounds Lower bounds of array dimensions (optional)
!> @param[in] ubounds Upper bounds of array dimensions (optional)
!> @param[in] unit Output unit number (optional, default: STDOUT)
!> @param[in] indent Indentation string (optional)
!> @param[in] sd If .true., print standard deviation (optional)
!> @enden
!>
!> @ja
!> @brief 7次元倍精度実数配列の要約を印字
!> @details
!> 配列のサイズ、最大値、最小値、平均値を印字します。
!> sd が .true. の場合、標準偏差も印字します。
!> @param[in] array 要約を印字する7次元配列
!> @param[in] lbounds 配列次元の下限（省略可能）
!> @param[in] ubounds 配列次元の上限（省略可能）
!> @param[in] unit 出力装置番号（省略可能、デフォルト: STDOUT）
!> @param[in] indent インデント文字列（省略可能）
!> @param[in] sd .true.の場合、標準偏差を印字（省略可能）
!> @endja
subroutine putlinedouble7( array, lbounds, ubounds, unit, indent, sd )
  use dc_types, only: dp, string, stdout
  use dc_string, only: tochar, printf, cprintf
  use dc_present, only: present_and_true
  implicit none
  real(DP), intent(in):: array(:,:,:,:,:,:,:)
  integer, intent(in), optional:: lbounds(7)
  integer, intent(in), optional:: ubounds(7)
  integer, intent(in), optional:: unit
  character(*), intent(in), optional:: indent
  logical, intent(in), optional:: sd
  integer:: out_unit
  integer:: indent_len
  character(STRING):: indent_str
  integer:: i
  integer:: alldim_size, lbound_nums(7), ubound_nums(7)
  character(STRING):: size_str, sd_str
  real(DP) :: max_value, min_value
  real(DP), allocatable :: array_packed(:)
  real(DP) :: avg_value, variance_value, sd_value
  continue
  !-----------------------------------------------------------------
  !  オプショナル引数のチェック
  !  Check optional arguments
  !-----------------------------------------------------------------
  if ( present(unit) ) then
    out_unit = unit
  else
    out_unit = stdout
  end if
  indent_len = 0
  indent_str = ''
  if ( present(indent) ) then
    if (len(indent) /= 0) then
      indent_len = len(indent)
      indent_str(1:indent_len) = indent
    end if
  end if
  !-------------------------------------------------------------------
  !  配列サイズ
  !  Array size
  !-------------------------------------------------------------------
  if ( present(lbounds) .and. present(ubounds) ) then
    lbound_nums = lbounds
    ubound_nums = ubounds
  else
    lbound_nums(1) = lbound( array, 1 )
    ubound_nums(1) = ubound( array, 1 )
    lbound_nums(2) = lbound( array, 2 )
    ubound_nums(2) = ubound( array, 2 )
    lbound_nums(3) = lbound( array, 3 )
    ubound_nums(3) = ubound( array, 3 )
    lbound_nums(4) = lbound( array, 4 )
    ubound_nums(4) = ubound( array, 4 )
    lbound_nums(5) = lbound( array, 5 )
    ubound_nums(5) = ubound( array, 5 )
    lbound_nums(6) = lbound( array, 6 )
    ubound_nums(6) = ubound( array, 6 )
    lbound_nums(7) = lbound( array, 7 )
    ubound_nums(7) = ubound( array, 7 )
  end if
  size_str = '('
  size_str = trim(size_str) // trim(tochar(lbound_nums(1))) // ':'
  size_str = trim(size_str) // trim(tochar(ubound_nums(1)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(2)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(2)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(3)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(3)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(4)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(4)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(5)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(5)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(6)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(6)))
  size_str = trim(size_str) // ',' // trim(tochar(lbound_nums(7)))
  size_str = trim(size_str) // ':' // trim(tochar(ubound_nums(7)))
  size_str = trim(size_str) // ')'
  !-------------------------------------------------------------------
  !  最大値
  !  Maximum value
  !-------------------------------------------------------------------
  max_value = maxval(array)
  !-------------------------------------------------------------------
  !  最小値
  !  Minimum value
  !-------------------------------------------------------------------
  min_value = minval(array)
  !-------------------------------------------------------------------
  !  平均値
  !  Average value
  !-------------------------------------------------------------------
  alldim_size = size(array)
  avg_value = sum(array) / dble(alldim_size)
  !-------------------------------------------------------------------
  !  標準偏差
  !  Standard deviation
  !-------------------------------------------------------------------
  sd_value = 0.0d0
  variance_value = 0.0d0
  sd_str = ''
  if ( present_and_true( sd ) ) then
    if ( alldim_size > 1 ) then
      if (allocated(array_packed)) then
        deallocate(array_packed)
      end if
      allocate( array_packed(alldim_size) )
      array_packed = pack(array, .true.)
      do i = 1, alldim_size
        variance_value = variance_value + &
          & (array_packed(i) - avg_value) * (array_packed(i) - avg_value)
      end do
      variance_value = variance_value / dble(alldim_size)
      sd_value = sqrt( variance_value )
      sd_str = cprintf( ' @sd=%f', d = (/ sd_value /) )
    end if
  end if
  !-------------------------------------------------------------------
  !  印字
  !  Print
  !-------------------------------------------------------------------
  call printf(out_unit, &
    & indent_str(1:indent_len) // &
    & '#<DP-ARRAY:: @size=%c, @max=%f, @min=%f, @avg=%f%c>', &
    & d = (/max_value, min_value, avg_value/), &
    & c1 = trim(size_str), c2 = trim(sd_str) )
end subroutine putlinedouble7