dc_clock.f90

Go to the documentation of this file.

! -*- mode: f90; coding: utf-8 -*-
!-----------------------------------------------------------------------
! Copyright (c) 2006-2026 Gtool Development Group. All rights reserved.
!-----------------------------------------------------------------------
module dc_clock
  use dc_types, only: string, dp
  use dc_trace, only: beginsub, endsub, dbgmessage
  use dc_date_types, only: dc_difftime, dc_datetime
  ! MPI ライブラリ
  ! MPI library
  !
  use mpi
  implicit none
  private
  public:: clock
  public:: dcclockcreate, dcclockclose, dcclockstart, dcclockstop
  public:: dcclockputline, dcclockresult, dcclocksetname
  public:: operator(+), operator(-)
  public:: dcclockget, dcclockevalsec, dcclocktochar, dcclockpredict
  !-----------------------------------------------
  ! 後方互換用
  ! For backward compatibility
  public:: create, Close, start, stop, putline, result, set_name
  public:: get, evalsec, tochar, predict
  type clock
    private
    character(STRING):: name
    real(dp):: start_time           ! 計測を開始した時間
                                    ! (計測の一時停止中には負の値が設定される)
    real(dp):: elapsed_time         ! 経過時間の累計値
    type(dc_datetime):: start_date  ! 計測を開始した日時
    logical:: initialized = .false. ! CLOCK 構造体の初期化チェック用フラグ
  end type clock
  interface dcclockcreate
    module procedure dcclockcreate0
  end interface
  interface dcclockclose
    module procedure dcclockclose0
  end interface
  interface dcclockstart
    module procedure dcclockstart0
  end interface
  interface dcclockstop
    module procedure dcclockstop0
  end interface
  interface dcclockputline
    module procedure dcclockputline0
  end interface
  interface dcclockget
    module procedure dcclockgetr
    module procedure dcclockgetd
  end interface
  interface dcclockevalsec
    module procedure dcclockevalsecd
  end interface
  interface dcclocktochar
    module procedure dcclocktochar0
  end interface
  interface dcclockresult
    module procedure dcclockresult0
  end interface
  interface operator(+)
    module procedure dcclockadd
  end interface
  interface operator(-)
    module procedure dcclocksubtract
  end interface
  interface dcclocksetname
    module procedure dcclocksetname0
  end interface
  interface dcclockpredict
    module procedure dcclockpredict0
  end interface
  !-----------------------------------------------
  ! 後方互換用
  ! For backward compatibility
  interface create
    module procedure dcclockcreate0
  end interface
  interface close
    module procedure dcclockclose0
  end interface
  interface start
    module procedure dcclockstart0
  end interface
  interface stop
    module procedure dcclockstop0
  end interface
  interface putline
    module procedure dcclockputline0
  end interface
  interface get
    module procedure dcclockgetr
    module procedure dcclockgetd
  end interface
  interface evalsec
    module procedure dcclockevalsecd
  end interface
  interface tochar
    module procedure dcclocktochar0
  end interface
  interface result
    module procedure dcclockresult0
  end interface
  interface set_name
    module procedure dcclocksetname0
  end interface
  interface predict
    module procedure dcclockpredict0
  end interface
  character(*), parameter:: version = &
    & '$Name:  $' // &
    & '$Id: dc_clock.F90,v 1.1 2009-03-20 09:09:53 morikawa Exp $'
contains
  subroutine dcclockcreate0(clk, name)
    use dc_message, only: messagenotify
    use dc_date, only: dcdatetimecreate
    implicit none
    type(clock), intent(out):: clk
    character(*), intent(in):: name
    character(*), parameter:: subname = 'DCClockCreate'
  continue
    call beginsub(subname, 'name=%c', c1=trim(name), version=version)
    if (clk % initialized) then
      call messagenotify('W', subname, 'This argument (type CLOCK) is already initialized.')
      call dbgmessage('already initialized')
      goto 999
    end if
    clk % name = name
    clk % elapsed_time = 0.0
    clk % start_time = - 1.0
    clk % initialized = .true.
    call dcdatetimecreate(clk % start_date)
    call dbgmessage('normal initialized')
999 continue
    call endsub(subname)
  end subroutine dcclockcreate0
  subroutine dcclockclose0(clk)
    implicit none
    type(clock), intent(inout):: clk
    character(*), parameter:: subname = 'DCClockClose'
  continue
    call beginsub(subname)
    if (clk % initialized) then
      clk % initialized = .false.
      clk % name = ''
    end if
    call endsub(subname)
  end subroutine dcclockclose0
  subroutine dcclockstart0(clk, err)
    use dc_message, only: messagenotify
    use dc_string, only: tochar
    use dc_error, only: storeerror, dc_enotinit, dc_noerr
    use dc_date, only: evalsec
    !$ use omp_lib
    implicit none
    type(clock), intent(inout):: clk
    logical, intent(out), optional:: err
    character(STRING):: cause_c
    integer:: stat
    character(*), parameter:: subname = 'DCClockStart'
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = 'CLOCK'
    if (.not. clk % initialized) then
      call messagenotify('W', subname, 'Call Create before Start in dc_clock.')
      call dbgmessage('Ignored because input argument was not initialized.')
      stat = dc_enotinit
      goto 999
    end if
    call cpu_time(clk % start_time) ! (out)
    !$ clk % start_time = omp_get_wtime()
    call dbgmessage('name=%c, cpu_time=%f', &
      & c1=trim(clk % name), d=(/clk % start_time/) )
999 continue
    call storeerror(stat, subname, err, cause_c)
    call endsub(subname)
  end subroutine dcclockstart0
  subroutine dcclockstop0(clk, err)
    use dc_message, only: messagenotify
    use dc_string, only: tochar
    use dc_error, only: storeerror, dc_enotinit, dc_noerr
    use dc_date, only: evalsec, operator(+), operator(-)
    use dc_date_types, only: dc_difftime
    use dc_types, only: dp
    !$ use omp_lib
    implicit none
    type(clock), intent(inout):: clk
    logical, intent(out), optional:: err
    character(STRING):: cause_c
    real(DP):: stop_time
    integer:: stat
    character(*), parameter:: subname = 'DCClockStop'
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = 'CLOCK'
    if (.not. clk % initialized) then
      call messagenotify('W', subname, 'Call Create before Stop in dc_clock.')
      call dbgmessage('Ignored because input argument was not initialized.')
      stat = dc_enotinit
      goto 999
    elseif (clk % start_time < 0.0_dp) then
      call messagenotify('W', subname, 'Call Start before Stop in dc_clock.')
      call dbgmessage('Ignored because input argument was not started.')
      goto 999
    end if
    call cpu_time(stop_time)
    !$ stop_time = omp_get_wtime()
    clk % elapsed_time = clk % elapsed_time + stop_time - clk % start_time
    clk % start_time = - 1.0
    call dbgmessage('name=%c, cpu_time=%f, elapsed_time=%f', &
      & c1=trim(clk % name), d=(/stop_time, clk % elapsed_time/))
999 continue
    call storeerror(stat, subname, err, cause_c)
    call endsub(subname)
  end subroutine dcclockstop0
  subroutine dcclockgetr(clk, sec, err) !:doc-priority 40:
    use dc_message, only: messagenotify
    use dc_date, only: evalsec
    use dc_string, only: cprintft
    use dc_error, only: storeerror, dc_enotinit, dc_noerr
    implicit none
    type(clock), intent(in):: clk
    real, intent(out):: sec
    logical, intent(out), optional:: err
    character(STRING):: cause_c
    integer:: stat
    character(*), parameter:: subname = 'DCClockGetR'
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = 'CLOCK'
    if (.not. clk % initialized) then
      call messagenotify('W', subname, 'Call Create before Get in dc_clock.')
      call dbgmessage('Ignored because input argument was not initialized.')
      stat = dc_enotinit
      goto 999
    end if
    sec = real(clk % elapsed_time, kind=kind(sec))
    call dbgmessage('name=%c, return sec=<%r>', &
      & c1=trim(clk % name), r=(/sec/))
999 continue
    call storeerror(stat, subname, err, cause_c)
    call endsub(subname)
  end subroutine dcclockgetr
  subroutine dcclockgetd(clk, sec, err) !:doc-priority 60:
    use dc_types, only: dp
    use dc_string, only: cprintf
    use dc_message, only: messagenotify
    use dc_date, only: evalsec
    use dc_error, only: storeerror, dc_enotinit, dc_noerr
    implicit none
    type(clock), intent(in):: clk
    real(DP), intent(out):: sec
    logical, intent(out), optional:: err
    character(STRING):: cause_c
    integer:: stat
    character(*), parameter:: subname = 'DCClockGetD'
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = 'CLOCK'
    if (.not. clk % initialized) then
      call messagenotify('W', subname, 'Call Create before Get in dc_clock.')
      call dbgmessage('Ignored because input argument was not initialized.')
      stat = dc_enotinit
      goto 999
    end if
    sec = clk % elapsed_time
    call dbgmessage('name=%c, return sec=<%f>', &
      & c1=trim(clk % name), d=(/sec/))
999 continue
    call storeerror(stat, subname, err, cause_c)
    call endsub(subname)
  end subroutine dcclockgetd
  function dcclockevalsecd(clk) result(result)
    use dc_types, only: dp
    implicit none
    type(clock), intent(in):: clk
    real(dp):: result
    logical:: err
  continue
    call dcclockgetd(clk, result, err)
    if (err) result = -1.0_dp
  end function dcclockevalsecd
  function dcclocktochar0(clk) result(result)
    use dc_string, only: cprintf
    use dc_date, only: evalsec
    implicit none
    type(clock), intent(in):: clk
    character(STRING):: result
    character(20):: clk_name
    integer:: name_len
  continue
    clk_name = ''
    name_len = min(len(clk_name), len_trim(clk % name))
    if (name_len > 0) then
      clk_name(1:name_len) = clk % name(1:name_len)
    end if
    if (clk % initialized) then
      result = cprintf(' %c%c  %c', c1 = clk_name, &
        & c2=trim(result_value_form(clk % elapsed_time)), &
        & c3=trim(fit_unit_value(clk % elapsed_time)))
    else
      result = ''
    end if
  end function dcclocktochar0
  subroutine dcclockputline0(clk, unit, indent, err)
    use dc_types, only: stdout
    use dc_message, only: messagenotify
    use dc_string, only: printf, tochar, cprintf
    use dc_date, only: evalsec, evalday, tochar
    use dc_error, only: storeerror, dc_enotinit, dc_noerr
    use dc_types, only: dp
    implicit none
    type(clock), intent(in):: clk
    integer, intent(in), optional:: unit
    character(*), intent(in), optional:: indent
    logical, intent(out), optional:: err
    integer:: out_unit
    character(STRING):: cause_c
    integer:: stat
    integer:: indent_len
    character(STRING):: indent_str
    character(*), parameter:: subname = 'DCClockPutLine'
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = 'CLOCK'
    if (.not. clk % initialized) then
      call messagenotify('W', subname, 'Call Create before PutLine in dc_clock.')
      call dbgmessage('Ignored because input argument was not initialized.')
      stat = dc_enotinit
      goto 999
    end if
    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
    call printf(out_unit, &
      & indent_str(1:indent_len) // &
      & '#<CLOCK:: @name=%c @clocking=%y @elapsed_time=%f sec. %c @start_date=%c>', &
      & c1=trim(clk % name), l=(/clk % start_time > 0.0_dp/), &
      & d=(/clk % elapsed_time/), &
      & c2=trim(fit_unit_value(clk % elapsed_time)), &
      & c3=trim(tochar(clk % start_date)))
    call dbgmessage('name=%c, output to device number %d', &
      & c1=trim(clk % name), i=(/out_unit/))
999 continue
    call storeerror(stat, subname, err, cause_c)
    call endsub(subname)
  end subroutine dcclockputline0
  subroutine dcclockresult0(clks, unit, total_auto, clk_total, total_name, err)
    use dc_types, only: stdout, string, dp
    use dc_message, only: messagenotify
    use dc_string, only: printf, tochar, cprintf
    use dc_date, only: evalsec
    use dc_error, only: storeerror, dc_enotinit, dc_noerr
    implicit none
    type(clock), intent(in):: clks(:)
    integer, intent(in), optional:: unit
    logical, intent(in), optional:: total_auto
    type(clock), intent(in), optional:: clk_total
    logical, intent(out), optional:: err
    character(*), intent(in), optional:: total_name
    integer:: out_unit, i, clks_size, ra
    character(20):: clk_name
    integer:: name_len
    character(STRING):: cause_c
    character(STRING):: total_name_work
    type(clock):: clk_auto_total
    logical:: total_print_complete
    real(DP):: elapsed_time_val_cor
    integer:: stat
    character(*), parameter:: total_time_mes = '       TOTAL TIME = '
    logical:: initflag_mpi
    integer:: err_mpi
    integer:: myrank_mpi, nprocs_mpi
    character(*), parameter:: subname = 'DCClockResult'
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = 'CLOCK'
    clks_size = size(clks)
    do i = 1, clks_size
      if (.not. clks(i) % initialized) then
        call messagenotify('W', subname, 'Call Create before Result in dc_clock.')
        call dbgmessage('Ignored because input argument was not initialized.')
        stat = dc_enotinit
        goto 999
      end if
    end do
    if (present(unit)) then
      out_unit = unit
    else
      out_unit = stdout
    end if
    if (present(total_name)) then
      total_name_work = ' (' // trim(total_name) // ')'
    else
      total_name_work = ''
    end if
    myrank_mpi = -1
    nprocs_mpi = 1
    call mpi_initialized(initflag_mpi, err_mpi)
    if ( initflag_mpi ) then
      call mpi_comm_rank(mpi_comm_world, myrank_mpi, err_mpi)
      call mpi_comm_size(mpi_comm_world, nprocs_mpi, err_mpi)
    end if
    do ra = 0, nprocs_mpi - 1
      if ( initflag_mpi ) call mpi_barrier(mpi_comm_world, err_mpi)
      if ( myrank_mpi > -1 .and. ra /= myrank_mpi ) cycle
      call printf(out_unit, '')
      if ( myrank_mpi < 0 ) then
        call printf(out_unit, &
          & ' ############## CPU TIME SUMMARY%c################', &
          & c1=trim(total_name_work) // ' ')
      else
        call printf(out_unit, &
          & ' ####### CPU TIME SUMMARY%c#### [rank=%06d] ####', &
          & c1=trim(total_name_work) // ' ', &
          & i = (/myrank_mpi/) )
      end if
      do i = 1, clks_size
        clk_name = ''
        name_len = min(len(clk_name), len_trim(clks(i) % name))
        if (name_len > 0) then
          clk_name(1:name_len) = clks(i) % name(1:name_len)
        end if
        elapsed_time_val_cor = clks(i) % elapsed_time
        if (elapsed_time_val_cor < 0.0_dp) elapsed_time_val_cor = 0.0_dp
        call printf(out_unit, &
          & ' %c%c  %c', c1=clk_name, &
          & c2=trim(result_value_form(elapsed_time_val_cor)), &
          & c3=trim(fit_unit_value(clks(i) % elapsed_time)))
      end do
      total_print_complete = .false.
      if (present(clk_total)) then
        if (clk_total % initialized) then
          call printf(out_unit, &
            & ' ------------------------------------------------')
          elapsed_time_val_cor = clk_total % elapsed_time
          if (elapsed_time_val_cor < 0.0_dp) elapsed_time_val_cor = 0.0_dp
          call printf(out_unit, &
            & ' %c%c  %c', c1=total_time_mes, &
            & c2=trim(result_value_form(elapsed_time_val_cor)), &
            & c3=trim(fit_unit_value(clk_total % elapsed_time)))
          total_print_complete = .true.
        end if
      end if
      if (present(total_auto) .and. .not. total_print_complete) then
        if (total_auto) then
          clk_auto_total = clks(1)
          if (clks_size > 1) then
            do i = 2, clks_size
              clk_auto_total = clk_auto_total + clks(i)
            end do
          end if
          call printf(out_unit, &
            & ' ------------------------------------------------')
          elapsed_time_val_cor = clk_auto_total % elapsed_time
          if (elapsed_time_val_cor < 0.0_dp) elapsed_time_val_cor = 0.0_dp
          call printf(out_unit, &
            & ' %c%c  %c', c1=total_time_mes, &
            & c2=trim(result_value_form(elapsed_time_val_cor)), &
            & c3=trim(fit_unit_value(clk_auto_total % elapsed_time)))
        end if
      end if
      call dbgmessage('total results, output to device number %d', &
        & i=(/out_unit/))
    end do
999 continue
    call storeerror(stat, subname, err, cause_c)
    call endsub(subname)
  end subroutine dcclockresult0
  function result_value_form(value) result(result)
    use dc_types, only: dp, token
    implicit none
    character(TOKEN):: result
    real(DP), intent(in):: value
  continue
    write(result, "(e15.6)") value
  end function result_value_form
  function fit_unit_value(sec, diff) result(result)
    use dc_types, only: dp, token
    use dc_date_types, only: dc_difftime
    use dc_date, only: dcdifftimecreate, evalday, evalhour, evalmin, evalsec
    implicit none
    character(TOKEN):: result
    real(DP), intent(in):: sec
    type(dc_difftime), intent(in), optional:: diff
    type(dc_difftime):: diffw
    character(TOKEN):: unit
    real(DP):: value
    character(TOKEN):: cval
  continue
    if ( present(diff) ) then
      diffw = diff
    else
      call dcdifftimecreate( diffw, sec = sec )
    end if
    if (evalday(diffw) > 1.0_dp) then
      unit = ' days'
      value = evalday(diffw)
    elseif (evalhour(diffw) > 1.0_dp) then
      unit = ' hrs.'
      value = evalhour(diffw)
    elseif (evalmin(diffw) > 1.0_dp) then
      unit = ' minutes'
      value = evalmin(diffw)
    else
      result = ''
      return
    end if
    cval = printf_g5_2(value)
    result = '(' // trim(adjustl(cval)) // trim(unit) // ')'
  end function fit_unit_value
  function printf_g5_2(value) result(result)
    use dc_types, only: dp, token
    use dc_string, only: cprintft
    implicit none
    character(TOKEN):: result
    real(DP), intent(in):: value
    character(TOKEN):: int_part, dem_part
    integer:: dem_int
  continue
    write(int_part, "(i20)") int(value)
    dem_int = nint((value - int(value)) * 100)
    if (dem_int < 0) dem_int = - dem_int
    if (dem_int == 100) then
      dem_int = 0
      write(int_part, "(i20)") int(value) + 1
    end if
    dem_part = cprintft('%02d', i=(/dem_int/))
    result = trim(adjustl(int_part)) // '.' // trim(dem_part)
  end function printf_g5_2
  function dcclockadd(clk1, clk2) result(clk_total)
    use dc_string, only: cprintf
    use dc_date, only: operator(+), operator(<)
    implicit none
    type(clock), intent(in):: clk1
    type(clock), intent(in):: clk2
    type(clock):: clk_total
  continue
    if (.not. clk1 % initialized .or. .not. clk2 % initialized) then
      clk_total % initialized = .false.
      return
    end if
    clk_total % name = cprintf('%c+%c', &
      & c1=trim(clk1 % name), c2=trim(clk2 % name))
    clk_total % start_time = - 1.0
    clk_total % initialized = .true.
    clk_total % elapsed_time = 0.0
    if (clk1 % start_date < clk2 % start_date) then
      clk_total % start_date = clk1 % start_date
    else
      clk_total % start_date = clk2 % start_date
    end if
    clk_total % elapsed_time = &
      & clk1 % elapsed_time + clk2 % elapsed_time
  end function dcclockadd
  function dcclocksubtract(clk1, clk2) result(clk_total)
    use dc_string, only: cprintf
    use dc_date, only: operator(-), operator(<)
    implicit none
    type(clock), intent(in):: clk1
    type(clock), intent(in):: clk2
    type(clock):: clk_total
  continue
    if (.not. clk1 % initialized .or. .not. clk2 % initialized) then
      clk_total % initialized = .false.
      return
    end if
    clk_total % name = cprintf('%c-%c', &
      & c1=trim(clk1 % name), c2=trim(clk2 % name))
    clk_total % start_time = - 1.0
    clk_total % initialized = .true.
    clk_total % elapsed_time = 0.0
    if (clk1 % start_date < clk2 % start_date) then
      clk_total % start_date = clk1 % start_date
    else
      clk_total % start_date = clk2 % start_date
    end if
    clk_total % elapsed_time = &
      & clk1 % elapsed_time - clk2 % elapsed_time
  end function dcclocksubtract
  subroutine dcclocksetname0(clk, name, err)
    use dc_message, only: messagenotify
    use dc_string, only: tochar, cprintf
    use dc_error, only: storeerror, dc_enotinit, dc_noerr
    implicit none
    type(clock), intent(inout):: clk
    character(*), intent(in):: name
    logical, intent(out), optional:: err
    character(STRING):: cause_c
    integer:: stat
    character(*), parameter:: subname = 'DCClockSetName'
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = 'CLOCK'
    if (.not. clk % initialized) then
      call messagenotify('W', subname, 'Call Create before Set_Name in dc_clock.')
      call dbgmessage('Ignored because input argument was not initialized.')
      stat = dc_enotinit
      goto 999
    end if
    clk % name = name
    call dbgmessage('set new name "%c"', c1=trim(clk % name))
999 continue
    call storeerror(stat, subname, err, cause_c)
    call endsub(subname)
  end subroutine dcclocksetname0
  subroutine dcclockpredict0(clk, progress, unit, err)
    use dc_types, only: stdout, dp
    use dc_message, only: messagenotify
    use dc_string, only: tochar, cprintf, printf
    use dc_error, only: storeerror, dc_enotinit, dc_noerr
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_date, only: operator(+), dcdatetimecreate, tochar, evalsec, &
      & dcdifftimecreate
    implicit none
    type(clock), intent(in):: clk
    real, intent(in):: progress
    integer, intent(in), optional:: unit
    logical, intent(out), optional:: err
    character(STRING):: cause_c
    integer:: stat, out_unit
    type(dc_difftime):: remain_diff
    type(dc_datetime):: comp_date, cur_date
    character(7):: prog_percent
    character(25):: prog_bar
    integer:: prog_bar_ptr
    real:: prog_valid
    logical:: initflag_mpi
    integer:: err_mpi, myrank_mpi
    character(*), parameter:: subname = 'DCClockPredict'
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = 'CLOCK'
    if (.not. clk % initialized) then
      call messagenotify('W', subname, 'Call Create before Predict in dc_clock.')
      call dbgmessage('Ignored because input argument was not initialized.')
      stat = dc_enotinit
      goto 999
    end if
    if (progress <= 0.0) then
      call messagenotify('W', subname, 'Specify 0.0 -- 1.0 value to "progress"')
      return
    elseif (progress > 1.0) then
      call messagenotify('W', subname, 'Over 1.0 value to "progress" was modified to 1.0')
      prog_valid = 1.0
    else
      prog_valid = progress
    end if
    if (present(unit)) then
      out_unit = unit
    else
      out_unit = stdout
    end if
    call dcdifftimecreate( remain_diff, &
      & sec = real(nint(evalsec(clk) / prog_valid * (1.0 - prog_valid)), dp) )
    call dcdatetimecreate(cur_date)
    comp_date = cur_date + remain_diff
    prog_percent = ''
    prog_percent = adjustr(trim(printf_g5_2(real(prog_valid * 100, dp))) // '%')
    prog_bar = ''
    prog_bar_ptr = int(prog_valid * 25)
    if (prog_bar_ptr > 0) prog_bar(1:prog_bar_ptr) = '*************************'
    call mpi_initialized(initflag_mpi, err_mpi)
    if ( initflag_mpi ) then
      call mpi_comm_rank(mpi_comm_world, myrank_mpi, err_mpi)
      if ( myrank_mpi /= 0 ) goto 999
    end if
    call printf(out_unit, '')
    call printf(out_unit, &
      & ' ########## PREDICTION OF CALCULATION ###########')
    call printf(out_unit, &
      & ' Start Date             %c', c1=trim(tochar(clk % start_date)))
    call printf(out_unit, &
      & ' Current Date           %c', c1=trim(tochar(cur_date)))
    call printf(out_unit, &
      & ' Progress     %c [%c]', c1=prog_percent, c2=prog_bar)
    call printf(out_unit, &
      & ' Remaining CPU TIME   %c %c', &
      & c1=trim(result_value_form(evalsec(remain_diff))), &
      & c2=trim(fit_unit_value(0.0_dp, remain_diff)))
    call printf(out_unit, &
      & ' Completion Date        %c', c1=trim(tochar(comp_date)))
999 continue
    call storeerror(stat, subname, err, cause_c)
    call endsub(subname)
  end subroutine dcclockpredict0
end module dc_clock