historyautoput.f90

Go to the documentation of this file.

!--
! *** Caution!! ***
!
! This file is generated from "historyautoput.rb2f90" by Ruby 3.3.8.
! Please do not edit this file directly.
!
! [JAPANESE]
!
! ※※※ 注意!!! ※※※
!
! このファイルは "historyautoput.rb2f90" から Ruby 3.3.8
! によって自動生成されたファイルです.
! このファイルを直接編集しませんようお願い致します.
!
!
!++
 
  subroutine historyautoputdouble0( &
    & time, varname, value, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                                            !
    ! データの出力を行います.
    ! このサブルーチンを用いる前に, "HistoryAutoCreate"
    ! による初期設定が必要です.
    !
    ! *varname* は HistoryAutoAddVariable で指定されている必要があります.
    !
    ! *HistoryAutoPut* は複数のサブルーチンの総称名です. *array* には
    ! 0 〜 7 次元のデータを与えることが可能です.
    ! (以下の同名のサブルーチンを参照ください).
    ! また, 整数, 単精度実数, 倍精度実数を与えることが可能です.
    ! ただし, 0 次元のデータを与える際の引数キーワードは
    ! *value* を用いてください.
    !
    ! Output data.
    ! Initialization by "HistoryAutoCreate" is needed
    ! before use of this subroutine.
    !
    ! "varname" must be specified by "HistoryAutoAddVariable".
    !
    ! "HistoryAutoPut" is a generic name of multiple subroutines.
    ! Then 0 -- 7 dimensional data can be given to "array".
    ! (See bellow subroutines with the same name).
    ! And, integer, sinble or double precision can be given.
    ! However, if 0 dimensional data is given, use "value" as a
    ! keyword argument.
    !
    !
    ! * 時間平均について
    !
    !   時間平均については HistoryAutoAddVariable を参照ください。
    !
    ! * About time average
    !
    !   See "HistoryAutoAddVariable" for details of time average
    !
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                                                                      ! データの時刻.
                              ! Time of data
                    
    character(*), intent(in):: varname
                                                                      ! 変数の名前.
                              !
                              ! ただし, ここで指定するものは,
                              ! HistoryAutoAddVariable の
                              ! *varname* で既に指定されてい
                              ! なければなりません.
                              !
                              ! Name of a variable.
                              !
                              ! This must be specified
                              ! *varname* in "HistoryAutoAddVariable".
                    
    real(DP), intent(in), target:: value
                                                                      ! 出力データ.
                              !
                              ! データ型は整数, 単精度実数型,
                              ! 倍精度実数型のどれでもかまいません.
                              ! ただし, ファイルへ出力される際には,
                              ! HistoryAutoAddVariable の *xtypes* で指定した
                              ! データ型へ変換されます.
                              !
                              ! Output data.
                              !
                              ! Integer, single or double precision are
                              ! acceptable as data type.
                              ! Note that when this is output to a file,
                              ! data type is converted into "xtype"
                              ! specified in "HistoryAutoAddVariable"
                              !
                    
    logical, intent(out), optional:: err
                                                                      ! 例外処理用フラグ.
                              ! デフォルトでは, この手続き内でエラーが
                              ! 生じた場合, プログラムは強制終了します.
                              ! 引数 *err* が与えられる場合,
                              ! プログラムは強制終了せず, 代わりに
                              ! *err* に .true. が代入されます.
                              !
                              ! Exception handling flag.
                              ! By default, when error occur in
                              ! this procedure, the program aborts.
                              ! If this *err* argument is given,
                              ! .true. is substituted to *err* and
                              ! the program does not abort.
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                                        
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble0"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        ! array only
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        ! array only
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                                        
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputdouble0
 
 
  subroutine historyautoputdouble1( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:) =>null()
    type(slice_info), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble1"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputdouble1
 
 
  subroutine historyautoputdouble2( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble2"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputdouble2
 
 
  subroutine historyautoputdouble3( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble3"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputdouble3
 
 
  subroutine historyautoputdouble4( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble4"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputdouble4
 
 
  subroutine historyautoputdouble5( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble5"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputdouble5
 
 
  subroutine historyautoputdouble6( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble6"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputdouble6
 
 
  subroutine historyautoputdouble7( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble7"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
!!$        write(*,*) '  sv%st(7), sv%ed(7), sv%sd(7)=', sv%st(7), sv%ed(7), sv%sd(7)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & , sv%st(7):sv%ed(7):sv%sd(7) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt7( sv%st(7):sv%ed(7):sv%sd(7) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputdouble7
 
 
  subroutine historyautoputreal0( &
    & time, varname, value, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                                        
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                                        
    character(*), intent(in):: varname
                                        
    real, intent(in), target:: value
                                        
    logical, intent(out), optional:: err
                                        
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                                        
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutReal0"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        ! array only
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        ! array only
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                                        
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputreal0
 
 
  subroutine historyautoputreal1( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real, intent(in), target:: array(:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real, pointer:: array_slice(:) =>null()
    type(slice_info), pointer:: sv =>null()
    real, pointer:: array_avr(:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutReal1"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputreal1
 
 
  subroutine historyautoputreal2( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real, intent(in), target:: array(:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real, pointer:: array_slice(:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real, pointer:: array_avr(:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutReal2"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputreal2
 
 
  subroutine historyautoputreal3( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real, intent(in), target:: array(:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real, pointer:: array_slice(:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real, pointer:: array_avr(:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutReal3"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputreal3
 
 
  subroutine historyautoputreal4( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real, intent(in), target:: array(:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real, pointer:: array_slice(:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real, pointer:: array_avr(:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutReal4"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputreal4
 
 
  subroutine historyautoputreal5( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real, intent(in), target:: array(:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real, pointer:: array_slice(:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real, pointer:: array_avr(:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutReal5"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputreal5
 
 
  subroutine historyautoputreal6( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real, intent(in), target:: array(:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real, pointer:: array_slice(:,:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real, pointer:: array_avr(:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutReal6"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputreal6
 
 
  subroutine historyautoputreal7( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real, intent(in), target:: array(:,:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real, pointer:: array_slice(:,:,:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    real, pointer:: array_avr(:,:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutReal7"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
!!$        write(*,*) '  sv%st(7), sv%ed(7), sv%sd(7)=', sv%st(7), sv%ed(7), sv%sd(7)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & , sv%st(7):sv%ed(7):sv%sd(7) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt7( sv%st(7):sv%ed(7):sv%sd(7) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputreal7
 
 
  subroutine historyautoputint0( &
    & time, varname, value, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                                        
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                                        
    character(*), intent(in):: varname
                                        
    integer, intent(in), target:: value
                                        
    logical, intent(out), optional:: err
                                        
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                                        
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt0"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        ! array only
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        ! array only
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                                        
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputint0
 
 
  subroutine historyautoputint1( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:) =>null()
    type(slice_info), pointer:: sv =>null()
    integer, pointer:: array_avr(:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt1"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputint1
 
 
  subroutine historyautoputint2( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt2"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputint2
 
 
  subroutine historyautoputint3( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt3"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputint3
 
 
  subroutine historyautoputint4( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt4"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputint4
 
 
  subroutine historyautoputint5( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt5"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputint5
 
 
  subroutine historyautoputint6( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt6"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputint6
 
 
  subroutine historyautoputint7( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(gt_history), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:,:,:,:,:) =>null()
    type(slice_info), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt7"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
!!$        write(*,*) '  sv%st(7), sv%ed(7), sv%sd(7)=', sv%st(7), sv%ed(7), sv%sd(7)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & , sv%st(7):sv%ed(7):sv%sd(7) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt7( sv%st(7):sv%ed(7):sv%sd(7) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)
  end subroutine historyautoputint7
 
 
  subroutine historyautoputold1double0( &
    & time, varname, value, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: value
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, value, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1double0
 
 
  subroutine historyautoputold1double1( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: array(:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1double1
 
 
  subroutine historyautoputold1double2( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: array(:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1double2
 
 
  subroutine historyautoputold1double3( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: array(:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1double3
 
 
  subroutine historyautoputold1double4( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: array(:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1double4
 
 
  subroutine historyautoputold1double5( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: array(:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1double5
 
 
  subroutine historyautoputold1double6( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: array(:,:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1double6
 
 
  subroutine historyautoputold1double7( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: array(:,:,:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1double7
 
 
  subroutine historyautoputold1real0( &
    & time, varname, value, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real, intent(in), target:: value
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, value, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1real0
 
 
  subroutine historyautoputold1real1( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real, intent(in), target:: array(:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1real1
 
 
  subroutine historyautoputold1real2( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real, intent(in), target:: array(:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1real2
 
 
  subroutine historyautoputold1real3( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real, intent(in), target:: array(:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1real3
 
 
  subroutine historyautoputold1real4( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real, intent(in), target:: array(:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1real4
 
 
  subroutine historyautoputold1real5( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real, intent(in), target:: array(:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1real5
 
 
  subroutine historyautoputold1real6( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real, intent(in), target:: array(:,:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1real6
 
 
  subroutine historyautoputold1real7( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    real, intent(in), target:: array(:,:,:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1real7
 
 
  subroutine historyautoputold1int0( &
    & time, varname, value, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    integer, intent(in), target:: value
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, value, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1int0
 
 
  subroutine historyautoputold1int1( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    integer, intent(in), target:: array(:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1int1
 
 
  subroutine historyautoputold1int2( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    integer, intent(in), target:: array(:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1int2
 
 
  subroutine historyautoputold1int3( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    integer, intent(in), target:: array(:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1int3
 
 
  subroutine historyautoputold1int4( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    integer, intent(in), target:: array(:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1int4
 
 
  subroutine historyautoputold1int5( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    integer, intent(in), target:: array(:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1int5
 
 
  subroutine historyautoputold1int6( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    integer, intent(in), target:: array(:,:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1int6
 
 
  subroutine historyautoputold1int7( &
    & time, varname, array, &  ! (in)
    & err &                                   ! (out) optional
    & )
    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(dc_difftime), intent(in):: time
    character(*), intent(in):: varname
    integer, intent(in), target:: array(:,:,:,:,:,:,:)
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, array, &  ! (in)
      & err )                                    ! (out) optional
 
  end subroutine historyautoputold1int7
 
!--
! vi:set readonly sw=4 ts=8:
!
!Local Variables:
!mode: f90
!buffer-read-only: t
!End:
!
!++