Output data . More...

Functions/Subroutines
subroutine	historyautoputdouble0 (time, varname, value, err)
subroutine	historyautoputdouble1 (time, varname, array, err)
subroutine	historyautoputdouble2 (time, varname, array, err)
subroutine	historyautoputdouble3 (time, varname, array, err)
subroutine	historyautoputdouble4 (time, varname, array, err)
subroutine	historyautoputdouble5 (time, varname, array, err)
subroutine	historyautoputdouble6 (time, varname, array, err)
subroutine	historyautoputdouble7 (time, varname, array, err)
subroutine	historyautoputreal0 (time, varname, value, err)
subroutine	historyautoputreal1 (time, varname, array, err)
subroutine	historyautoputreal2 (time, varname, array, err)
subroutine	historyautoputreal3 (time, varname, array, err)
subroutine	historyautoputreal4 (time, varname, array, err)
subroutine	historyautoputreal5 (time, varname, array, err)
subroutine	historyautoputreal6 (time, varname, array, err)
subroutine	historyautoputreal7 (time, varname, array, err)
subroutine	historyautoputint0 (time, varname, value, err)
subroutine	historyautoputint1 (time, varname, array, err)
subroutine	historyautoputint2 (time, varname, array, err)
subroutine	historyautoputint3 (time, varname, array, err)
subroutine	historyautoputint4 (time, varname, array, err)
subroutine	historyautoputint5 (time, varname, array, err)
subroutine	historyautoputint6 (time, varname, array, err)
subroutine	historyautoputint7 (time, varname, array, err)
subroutine	historyautoputold1double0 (time, varname, value, err)
subroutine	historyautoputold1double1 (time, varname, array, err)
subroutine	historyautoputold1double2 (time, varname, array, err)
subroutine	historyautoputold1double3 (time, varname, array, err)
subroutine	historyautoputold1double4 (time, varname, array, err)
subroutine	historyautoputold1double5 (time, varname, array, err)
subroutine	historyautoputold1double6 (time, varname, array, err)
subroutine	historyautoputold1double7 (time, varname, array, err)
subroutine	historyautoputold1real0 (time, varname, value, err)
subroutine	historyautoputold1real1 (time, varname, array, err)
subroutine	historyautoputold1real2 (time, varname, array, err)
subroutine	historyautoputold1real3 (time, varname, array, err)
subroutine	historyautoputold1real4 (time, varname, array, err)
subroutine	historyautoputold1real5 (time, varname, array, err)
subroutine	historyautoputold1real6 (time, varname, array, err)
subroutine	historyautoputold1real7 (time, varname, array, err)
subroutine	historyautoputold1int0 (time, varname, value, err)
subroutine	historyautoputold1int1 (time, varname, array, err)
subroutine	historyautoputold1int2 (time, varname, array, err)
subroutine	historyautoputold1int3 (time, varname, array, err)
subroutine	historyautoputold1int4 (time, varname, array, err)
subroutine	historyautoputold1int5 (time, varname, array, err)
subroutine	historyautoputold1int6 (time, varname, array, err)
subroutine	historyautoputold1int7 (time, varname, array, err)

Detailed Description

Output data

.

Author: Yasuhiro MORIKAWA

Copyright: Copyright (C) GFD Dennou Club, 2008-2026. All rights reserved.
License is BSD-2-Clause. See COPYRIGHT in detail

Definition in file historyautoput.f90.

Function/Subroutine Documentation

◆ historyautoputdouble0()

subroutine historyautoputdouble0	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), intent(in), target	value,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 109 of file historyautoput.f90.

    !
                                            !
    ! データの出力を行います.
    ! このサブルーチンを用いる前に, "HistoryAutoCreate"
    ! による初期設定が必要です.
    !
    ! *varname* は HistoryAutoAddVariable で指定されている必要があります.
    !
    ! *HistoryAutoPut* は複数のサブルーチンの総称名です. *array* には
    ! 0 〜 7 次元のデータを与えることが可能です.
    ! (以下の同名のサブルーチンを参照ください).
    ! また, 整数, 単精度実数, 倍精度実数を与えることが可能です.
    ! ただし, 0 次元のデータを与える際の引数キーワードは
    ! *value* を用いてください.
    !
    ! Output data.
    ! Initialization by "HistoryAutoCreate" is needed
    ! before use of this subroutine.
    !
    ! "varname" must be specified by "HistoryAutoAddVariable".
    !
    ! "HistoryAutoPut" is a generic name of multiple subroutines.
    ! Then 0 -- 7 dimensional data can be given to "array".
    ! (See bellow subroutines with the same name).
    ! And, integer, sinble or double precision can be given.
    ! However, if 0 dimensional data is given, use "value" as a
    ! keyword argument.
    !
    !
    ! * 時間平均について
    !
    !   時間平均については HistoryAutoAddVariable を参照ください。
    !
    ! * About time average
    !
    !   See "HistoryAutoAddVariable" for details of time average
    !
    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                                                                      ! データの時刻.
                              ! Time of data
                    
    character(*), intent(in):: varname
                                                                      ! 変数の名前.
                              !
                              ! ただし, ここで指定するものは,
                              ! HistoryAutoAddVariable の
                              ! *varname* で既に指定されてい
                              ! なければなりません.
                              !
                              ! Name of a variable.
                              !
                              ! This must be specified
                              ! *varname* in "HistoryAutoAddVariable".
                    
    real(DP), intent(in), target:: value
                                                                      ! 出力データ.
                              !
                              ! データ型は整数, 単精度実数型,
                              ! 倍精度実数型のどれでもかまいません.
                              ! ただし, ファイルへ出力される際には,
                              ! HistoryAutoAddVariable の *xtypes* で指定した
                              ! データ型へ変換されます.
                              !
                              ! Output data.
                              !
                              ! Integer, single or double precision are
                              ! acceptable as data type.
                              ! Note that when this is output to a file,
                              ! data type is converted into "xtype"
                              ! specified in "HistoryAutoAddVariable"
                              !
                    
    logical, intent(out), optional:: err
                                                                      ! 例外処理用フラグ.
                              ! デフォルトでは, この手続き内でエラーが
                              ! 生じた場合, プログラムは強制終了します.
                              ! 引数 *err* が与えられる場合,
                              ! プログラムは強制終了せず, 代わりに
                              ! *err* に .true. が代入されます.
                              !
                              ! Exception handling flag.
                              ! By default, when error occur in
                              ! this procedure, the program aborts.
                              ! If this *err* argument is given,
                              ! .true. is substituted to *err* and
                              ! the program does not abort.
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                                        
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble0"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        ! array only
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        ! array only
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                                        
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputdouble1()

subroutine historyautoputdouble1	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 521 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble1"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputdouble2()

subroutine historyautoputdouble2	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 877 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble2"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputdouble3()

subroutine historyautoputdouble3	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 1238 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble3"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputdouble4()

subroutine historyautoputdouble4	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 1605 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble4"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputdouble5()

subroutine historyautoputdouble5	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 1978 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble5"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputdouble6()

subroutine historyautoputdouble6	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 2357 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:,:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble6"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputdouble7()

subroutine historyautoputdouble7	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 2742 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    real(DP), intent(in), target:: array(:,:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        real(DP), pointer:: array_slice(:,:,:,:,:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    real(DP), pointer:: array_avr(:,:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutDouble7"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
!!$        write(*,*) '  sv%st(7), sv%ed(7), sv%sd(7)=', sv%st(7), sv%ed(7), sv%sd(7)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & , sv%st(7):sv%ed(7):sv%sd(7) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt7( sv%st(7):sv%ed(7):sv%sd(7) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputint0()

subroutine historyautoputint0	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		integer, intent(in), target	value,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 6078 of file historyautoput.f90.

    !
                                        
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                                        
    character(*), intent(in):: varname
                                        
    integer, intent(in), target:: value
                                        
    logical, intent(out), optional:: err
                                        
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                                        
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt0"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        ! array only
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        ! array only
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & (/value/), &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                                        
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputint1()

subroutine historyautoputint1	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 6411 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    integer, pointer:: array_avr(:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt1"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputint2()

subroutine historyautoputint2	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 6767 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt2"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputint3()

subroutine historyautoputint3	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 7128 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt3"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputint4()

subroutine historyautoputint4	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 7495 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt4"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputint5()

subroutine historyautoputint5	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 7868 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt5"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputint6()

subroutine historyautoputint6	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 8247 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:,:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt6"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputint7()

subroutine historyautoputint7	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Output data.

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 below subroutines with the same name). And, integer, single or double precision can be given. However, if 0 dimensional data is given, use value as a keyword argument.

About time average

See HistoryAutoAddVariable for details of time average.

Parameters

[in]	time	Time of data.
[in]	varname	Name of a variable. This must be specified varname in HistoryAutoAddVariable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array). 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.
[out]	err	(optional) 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.

Definition at line 8632 of file historyautoput.f90.

    !
                    
    !
 
    use gtool_historyauto_internal, only: initialized, &
      & numvars, gthstnml, zero_time, &
      & create_timing_vars, close_timing_vars, &
      & histaddvar_vars, interval_time_vars, interval_unitsym_vars, &
      & newfile_createtime_vars, &
      & output_timing_avr_vars, output_timing_vars, prev_outtime_vars, &
      & renew_timing_vars, varname_vars, gthst_history_vars, slice_vars, &
      & space_avr_vars, weight_vars, cal_save, flag_output_prev_vars, &
      & origin_time_vars
    use gtool_historyauto_internal, only: hstvarsoutputcheck, hstfilecreate, &
      & averagereduce
    use gtool_historyauto_internal, only: slice_info
    use gtool_history_nmlinfo_generic, only: hstnmlinfoenddefine, &
      & hstnmlinfodefinemode
    use gtool_history, only: gt_history, historyput, historyinitialized, &
      & historyaddvariable, historyinquire, historyvarinfoinquire, &
      & historyclose, historysettime
    use dc_string, only: tochar
    use dc_message, only: messagenotify
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, &
      & hst_ebadvarname, dc_enotinit, dc_enegative
    use dc_calendar, only: unit_symbol_sec, dccalconvertbyunit
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline
    use dc_date_types, only: dc_difftime, dc_datetime
    use dc_types, only: dp, string
    implicit none
    real(DP), intent(in):: time
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
                    
    character(*), intent(in):: varname
                    
    integer, intent(in), target:: array(:,:,:,:,:,:,:)
                    
    logical, intent(out), optional:: err
                    
 
    type(GT_HISTORY), pointer:: gthist =>null()
                              ! gtool_history モジュール用構造体. 
                              ! Derived type for "gtool_history" module
 
                        integer, pointer:: array_slice(:,:,:,:,:,:,:) =>null()
    type(SLICE_INFO), pointer:: sv =>null()
    integer, pointer:: array_avr(:,:,:,:,:,:,:) =>null()
                    
    real(DP):: settime
    integer:: stat, i
    integer:: vnum
    character(STRING):: cause_c
    integer, save:: svnum = 1, svtstep
    character(*), parameter:: subname = "HistoryAutoPutInt7"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      stat = dc_enotinit
      cause_c = 'gtool_historyauto'
      goto 999
    end if
 
    ! 時刻に関するエラー処理
    ! Error handling for time
    !
    if ( time < zero_time ) then
      cause_c = tochar( time )
      call messagenotify( 'W', subname, &
        & '"time=<%c>" must be positive value (varname=<%c>).', &
        & c1 = trim( cause_c ), c2 = trim( varname ) )
      stat = dc_enegative
      cause_c = 'time'
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 定義モードからデータモードへ
    ! Transit from define mode to data mode
    !
    if ( hstnmlinfodefinemode( gthstnml ) ) then
      call hstnmlinfoenddefine( gthstnml ) ! (inout)
    end if
 
    ! 出力タイミングのチェックとファイルの作成
    ! Check output timing and create files
    !
    call hstvarsoutputcheck( &
      &        time = time, &    ! (in)
      & stime_index = svtstep )  ! (out)
 
    ! ファイルのオープン・クローズ・再オープン
    ! Open, close, reopen files
    !
    if ( create_timing_vars(vnum, svtstep) &
      &  .and. .not. histaddvar_vars(vnum)   ) then
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      histaddvar_vars(vnum) = .true.
      if ( flag_output_prev_vars(vnum) ) then
        prev_outtime_vars(vnum) = time
      else
        flag_output_prev_vars(vnum) = .true.
        if ( origin_time_vars(vnum) > zero_time ) then
          prev_outtime_vars(vnum) = origin_time_vars(vnum)
        else
          prev_outtime_vars(vnum) = time
        end if
      end if
    end if
 
    if ( close_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
    end if
 
    if ( renew_timing_vars(vnum, svtstep) ) then
      if ( historyinitialized( gthst_history_vars(vnum) % gthist ) ) then
        call historyclose( gthst_history_vars(vnum) % gthist )   ! (inout)
      end if
      call hstfilecreate( gthst_history_vars(vnum) % gthist, & ! (inout)
        & varname_vars(vnum), time )                           ! (in)
      newfile_createtime_vars(vnum) = time
      prev_outtime_vars(vnum) = time
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if (       .not. output_timing_vars(vnum, svtstep) &
      &  .and. .not. output_timing_avr_vars(vnum, svtstep) ) then
 
      goto 999
    end if
 
    ! GT_HISTORY 変数の取得
    ! Get "GT_HISTORY" variable
    !
    gthist => gthst_history_vars(vnum) % gthist
 
 
    ! 空間切り出し
    ! Slice of spaces
    !
                        sv => slice_vars(vnum)
 
                      !!$        write(*,*) '  sv%st(1), sv%ed(1), sv%sd(1)=', sv%st(1), sv%ed(1), sv%sd(1)
                      
!!$        write(*,*) '  sv%st(2), sv%ed(2), sv%sd(2)=', sv%st(2), sv%ed(2), sv%sd(2)
                      
!!$        write(*,*) '  sv%st(3), sv%ed(3), sv%sd(3)=', sv%st(3), sv%ed(3), sv%sd(3)
                      
!!$        write(*,*) '  sv%st(4), sv%ed(4), sv%sd(4)=', sv%st(4), sv%ed(4), sv%sd(4)
                      
!!$        write(*,*) '  sv%st(5), sv%ed(5), sv%sd(5)=', sv%st(5), sv%ed(5), sv%sd(5)
                      
!!$        write(*,*) '  sv%st(6), sv%ed(6), sv%sd(6)=', sv%st(6), sv%ed(6), sv%sd(6)
                      
!!$        write(*,*) '  sv%st(7), sv%ed(7), sv%sd(7)=', sv%st(7), sv%ed(7), sv%sd(7)
                      
 
    array_slice => array(sv%st(1):sv%ed(1):sv%sd(1) &
                            & , sv%st(2):sv%ed(2):sv%sd(2) &
                      
      & , sv%st(3):sv%ed(3):sv%sd(3) &
                      
      & , sv%st(4):sv%ed(4):sv%sd(4) &
                      
      & , sv%st(5):sv%ed(5):sv%sd(5) &
                      
      & , sv%st(6):sv%ed(6):sv%sd(6) &
                      
      & , sv%st(7):sv%ed(7):sv%sd(7) &
                      
      & )
                    
 
 
    ! 空間平均
    ! Spatial average
    !
                        if ( count(space_avr_vars(vnum) % avr) == 0 ) then
      array_avr => array_slice
    else
      call averagereduce( &
        & array_slice, space_avr_vars(vnum) % avr, &   ! (in)
                              & weight_vars(vnum) % wgt1( sv%st(1):sv%ed(1):sv%sd(1) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt2( sv%st(2):sv%ed(2):sv%sd(2) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt3( sv%st(3):sv%ed(3):sv%sd(3) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt4( sv%st(4):sv%ed(4):sv%sd(4) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt5( sv%st(5):sv%ed(5):sv%sd(5) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt6( sv%st(6):sv%ed(6):sv%sd(6) ) , &   ! (in)
                      
        & weight_vars(vnum) % wgt7( sv%st(7):sv%ed(7):sv%sd(7) ) , &   ! (in)
                      
        & array_avr )                       ! (out)
    end if
 
    ! 座標重みを取得 ; Get weights of axes
 
                    
 
    ! 単位に応じて時刻を変換
    ! Convert time according to units
    !
    if (      output_timing_vars(vnum, svtstep) &
      &  .or. output_timing_avr_vars(vnum, svtstep) ) then
 
      settime = &
        & dccalconvertbyunit( time, unit_symbol_sec, interval_unitsym_vars(vnum), cal_save )
    end if
 
    ! 時刻設定
    ! Set time
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      call historysettime( &
        & history = gthist, &   ! (in) optional
        & timed = settime )     ! (in) optional
    end if
 
    ! 出力
    ! OutPut
    !
    if ( output_timing_avr_vars(vnum, svtstep) ) then
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & timed = settime, &      ! (in) optional
        & time_average_store = &
        &   .not. output_timing_vars(vnum, svtstep), & ! (in) optional
        & history = gthist )                  ! (inout) optional
    else
      call historyput( &
        & varname, &              ! (in)
                              & array_avr, &            ! (in)
                      
        & history = gthist )      ! (inout) optional
    end if
 
    ! 最後に出力した時刻を保存
    ! Save last time of output
    !
    if ( output_timing_vars(vnum, svtstep) ) then
      if ( .not. create_timing_vars(vnum, svtstep) .and. &
        &  .not. renew_timing_vars(vnum, svtstep) ) then
 
        if ( abs(mod(time, interval_time_vars(vnum)) - zero_time) <= &
          & time_eps * max(1.0_dp, abs(time), &
          & abs(interval_time_vars(vnum)), abs(zero_time)) ) then
          prev_outtime_vars(vnum) = time
        else
          prev_outtime_vars(vnum) = &
            & prev_outtime_vars(vnum) + interval_time_vars(vnum)
        end if
      end if
    end if
 
    ! 結合解除
    ! Release associations
    !
    nullify( gthist )
                        nullify( array_avr, array_slice )
                    
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c, err = err)
    call endsub(subname)

Here is the call graph for this function:

◆ historyautoputold1double0()

subroutine historyautoputold1double0	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), intent(in), target	value,
		logical, intent(out), optional	err )

Output data (DC_DIFFTIME version).

This is a wrapper that accepts DC_DIFFTIME type for time argument.

Parameters

[in]	time	Time of data (DC_DIFFTIME type).
[in]	varname	Name of a variable.
[in]	value	Output data (for scalar).
[in]	array	Output data (for array).
[out]	err	(optional) Exception handling flag.

Definition at line 8968 of file historyautoput.f90.

    use gtool_historyauto_generic, only: historyautoput
    use dc_date_generic, only: operator(-), operator(+), &
      & operator(>), operator(<), mod, operator(==), tochar, evalsec, &
      & dcdifftimeputline, assignment(=)
    use dc_types, only: dp
    use dc_date_types, only: dc_difftime, dc_datetime
    implicit none
    type(DC_DIFFTIME), intent(in):: time
    character(*), intent(in):: varname
    real(DP), intent(in), target:: value
    logical, intent(out), optional:: err
 
    real(DP):: timed
 
  continue
    timed = evalsec( time )
 
    call historyautoput( &
      & timed, varname, value, &  ! (in)
      & err )                                    ! (out) optional
 

subroutine historyautoputold1double1	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1double2	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1double3	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1double4	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1double5	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1double6	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1double7	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real(dp), dimension(:,:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1int0	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		integer, intent(in), target	value,
		logical, intent(out), optional	err )

subroutine historyautoputold1int1	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1int2	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1int3	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1int4	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1int5	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1int6	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1int7	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		integer, dimension(:,:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1real0	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real, intent(in), target	value,
		logical, intent(out), optional	err )

subroutine historyautoputold1real1	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1real2	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1real3	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1real4	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ historyautoputdouble0()

About time average

◆ historyautoputdouble1()

About time average

◆ historyautoputdouble2()

About time average

◆ historyautoputdouble3()

About time average

◆ historyautoputdouble4()

About time average

◆ historyautoputdouble5()

About time average

◆ historyautoputdouble6()

About time average

◆ historyautoputdouble7()

About time average

◆ historyautoputint0()

About time average

◆ historyautoputint1()

About time average

◆ historyautoputint2()

About time average

◆ historyautoputint3()

About time average

◆ historyautoputint4()

About time average

◆ historyautoputint5()

About time average

◆ historyautoputint6()

About time average

◆ historyautoputint7()

About time average

◆ historyautoputold1double0()

◆ historyautoputold1double1()

◆ historyautoputold1double2()

◆ historyautoputold1double3()

◆ historyautoputold1double4()

◆ historyautoputold1double5()

◆ historyautoputold1double6()

◆ historyautoputold1double7()

◆ historyautoputold1int0()

◆ historyautoputold1int1()

◆ historyautoputold1int2()

◆ historyautoputold1int3()

◆ historyautoputold1int4()

◆ historyautoputold1int5()

◆ historyautoputold1int6()

◆ historyautoputold1int7()

◆ historyautoputold1real0()

◆ historyautoputold1real1()

◆ historyautoputold1real2()

◆ historyautoputold1real3()

◆ historyautoputold1real4()

◆ historyautoputold1real5()

◆ historyautoputold1real6()

◆ historyautoputold1real7()

◆ historyautoputreal0()

About time average

◆ historyautoputreal1()

About time average

◆ historyautoputreal2()

About time average

◆ historyautoputreal3()

About time average

◆ historyautoputreal4()

About time average

◆ historyautoputreal5()

About time average

◆ historyautoputreal6()

About time average

◆ historyautoputreal7()

About time average

subroutine historyautoputold1real5	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1real6	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputold1real7	(	type(dc_difftime), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputreal0	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real, intent(in), target	value,
		logical, intent(out), optional	err )

subroutine historyautoputreal1	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputreal2	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputreal3	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputreal4	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputreal5	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )

subroutine historyautoputreal6	(	real(dp), intent(in)	time,
		character(*), intent(in)	varname,
		real, dimension(:,:,:,:,:,:), intent(in), target	array,
		logical, intent(out), optional	err )