historyautoput.f90

Go to the documentation of this file.

!--
! *** Caution!! ***
!
! This file is generated from "historyautoput.rb2f90" by Ruby 2.3.3.
! Please do not edit this file directly.
!
! [JAPANESE]
!
! ※※※ 注意!!! ※※※
!
! このファイルは "historyautoput.rb2f90" から Ruby 2.3.3
! によって自動生成されたファイルです.
! このファイルを直接編集しませんようお願い致します.
!
!
!++
!
!= データ出力
!= Output data
!
! Authors::   Yasuhiro MORIKAWA
! Version::   $Id: historyautoput.rb2f90,v 1.6 2010-07-04 22:01:50 morikawa Exp $
! Tag Name::  $Name:  $
! Copyright:: Copyright (C) GFD Dennou Club, 2008-2009. All rights reserved. 
! License::   See COPYRIGHT[link:../../../COPYRIGHT]
!

  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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                                                                      ! データの時刻. 
                              ! 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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                                        
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                                        
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, numdims, &
      & 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, gt_eargsizemismatch, &
      & 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, token
    implicit none
    real(DP), intent(in):: time
                    
    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 ( mod(time, interval_time_vars(vnum) ) == 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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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, string, token
    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:
!
!++