gtool_historyauto_internal.f90

Go to the documentation of this file.

!--
! *** Caution!! ***
!
! This file is generated from "gtool_historyauto_internal.rb2f90" by Ruby 3.3.8.
! Please do not edit this file directly.
!
! [JAPANESE]
!
! ※※※ 注意!!! ※※※
!
! このファイルは "gtool_historyauto_internal.rb2f90" から Ruby 3.3.8
! によって自動生成されたファイルです.
! このファイルを直接編集しませんようお願い致します.
!
!
!++
 
module gtool_historyauto_internal
 
  ! gtool_history モジュール
  ! "gtool_history" module
  !
  use gtool_history, only: gt_history_axis, gt_history_varinfo, gt_history
 
  ! NAMELIST の使用を想定したヒストリデータ出力情報管理用ユーティリティ
  ! Utilities for history data output information management assuming use of NAMELIST
  !
  use gtool_history_nmlinfo, only: gthst_nmlinfo
 
  ! NetCDF ライブラリで規定される最大の次元の数
  ! Maximum number of dimensions prescribed by the NetCDF library
  !
  use netcdf, only: nf90_max_dims, nf90_max_vars
 
  ! 暦と日付の取り扱い
  ! Calendar and date handler
  !
  use dc_calendar, only: dc_cal, dc_cal_date
                              ! 暦と日時を表現するデータ型.
                              ! Derived data type for calendar and date
 
  ! 日付および時刻の取り扱い (旧版)
  ! Date and time handler (Old version)
  !
  use dc_date_types, only: dc_datetime, dc_difftime
                              ! 日時の差を表現するデータ型.
                              ! Data type for difference about date and time
  ! 種別型パラメタ
  ! Kind type parameter
  !
  use dc_types, only: dp, &      ! 倍精度実数型. Double precision.
    &                 string, &  ! 文字列.       Strings.
    &                 token, &   ! キーワード.   Keywords.
    &                 stderr     ! 標準エラー出力.  Standard error output
 
  implicit none
  private
  public:: gt_history_axis_data, gt_history_multi
  public:: slice_info, axes_weight, space_avr_info
 
  public:: hstvarsoutputcheck, hstfilecreate, averagereduce
 
  ! 次元数
  ! Number of dimensions
  !
  integer, save, public:: numdims
 
  ! 座標重み変数の数
  ! Number of variables of axes weight
  !
  integer, save, public:: numwgts = 0
 
  ! 座標重み変数の接尾詞
  ! Suffix of variables of axes weight
  !
  character(*), parameter, public:: wgtsuf = '_weight'
 
  ! 変数の数
  ! Number of variables
  !
  integer, save, public:: numvars = 0
 
  ! 1 つの変数が依存可能な次元の数
  ! Number of dimensions on which one variable can depend
  !
  integer, parameter, public:: max_dims_depended_by_var = 7
 
  ! 出力間隔をΔtで割った際に, 余りとして許容される範囲
  ! Allowable range of remainder of output interval divided by delta t
  !
  real(dp), parameter, public:: max_remainder_range = 1.0e-3_dp
 
  ! 出力ファイルの基本メタデータ
  ! Basic meta data for output file
  !
  character(STRING), save, public:: title_save
  character(STRING), save, public:: source_save
  character(STRING), save, public:: institution_save
  character(STRING), save, public:: conventions_save
  character(TOKEN),  save, public:: gt_version_save
  character(TOKEN),  save, public:: rank_save
 
  ! 時刻データ
  ! Time data
  !
!  integer, parameter, public:: save_tstepnum = 3
  integer, parameter, public:: save_tstepnum = 1
                              ! 保存する時間ステップの数.
                              ! Number of saved time step
  real(dp), save, public:: saved_time(1:save_tstepnum)
  integer, save, public:: checked_tstepnum = 0
                              ! チェックされた時間ステップの数.
                              ! Number of checked time step
  integer, save, public:: checked_tstep_varnum = 0
                              ! チェックされた変数の数.
                              ! Number of checked variables
  integer, save, public:: saved_tstep = 1
                              ! 前回チェックされた時間ステップ.
                              ! (HstVarsOutputCheck で使用する).
                              !
                              ! Time step checked at previous time
                              ! (Used in "HstVarsOutputCheck").
 
  ! 時刻の単位 (HistoryAutoCreate の dims によって指定されたもの)
  ! Unit of time (specified by "dims" of "HistoryAutoCreate")
  !
  character(TOKEN), save, public:: time_unit_bycreate = ''
  character(STRING), save, public:: time_unit_suffix = ''
 
  ! NAMELIST の使用を想定したヒストリデータ出力情報管理用ユーティリティ
  ! Utilities for history data output information management assuming use of NAMELIST
  !
  type(gthst_nmlinfo), save, public:: gthstnml
 
  type gt_history_axis_data
    real(dp), pointer:: a_axis(:) =>null()
  end type gt_history_axis_data
 
  ! 座標軸情報
  ! Information of axes
  !
  type(gt_history_axis), save, target, public:: gthst_axes(1:nf90_max_dims)
  type(gt_history_axis_data), save, target, public:: data_axes(1:nf90_max_dims)
  type(gt_history_axis_data), save, target, public:: data_axes_whole(1:nf90_max_dims)
  type(gt_history_varinfo), save, public:: gthst_weights(1:nf90_max_dims)
  type(gt_history_axis_data), save, target, public:: data_weights(1:nf90_max_dims)
 
  ! MPI 関連の情報
  ! Information about MPI
  !
  logical, save, public:: save_mpi_split  = .false.
  logical, save, public:: save_mpi_gather = .false.
 
  type gt_history_multi
    type(gt_history), pointer:: gthist =>null()
  end type gt_history_multi
 
  ! 変数情報
  ! Information of variables
  !
  integer, parameter, public:: max_vars = nf90_max_vars
                              ! 出力可能な変数の最大値
                              ! Maximum value of output variables
  type(gt_history_varinfo), save, public:: gthst_vars(1:max_vars)
  character(TOKEN), save, public:: varname_vars(1:max_vars) = ''
  type(gt_history_multi), save, public:: gthst_history_vars(1:max_vars)
  logical, save, public:: output_valid_vars(1:max_vars) = .false.
                              ! 変数出力が有効か否か.
                              ! Whether output of variables is valid or not.
  logical, save, public:: create_timing_vars(1:max_vars, 1:save_tstepnum) = .false.
                              ! 各時間ステップではファイルを作成するか
                              ! Whether file is created or not at eath time step.
  logical, save, public:: close_timing_vars(1:max_vars, 1:save_tstepnum) = .false.
                              ! 各時間ステップではファイルをクローズするか
                              ! Whether file is closed or not at eath time step.
  logical, save, public:: renew_timing_vars(1:max_vars, 1:save_tstepnum) = .false.
                              ! 各時間ステップではファイルを再オープンするか
                              ! Whether file is closed and opened or not at eath time step.
 
  logical, save, public:: output_timing_vars(1:max_vars, 1:save_tstepnum) = .false.
                              ! 各時間ステップでは出力を行うか否か.
                              ! Whether output is done or not at eath time step.
  logical, save, public:: output_timing_avr_vars(1:max_vars, 1:save_tstepnum) = .false.
                              ! 各時間ステップでは平均値出力を行うか否か.
                              ! Whether output of averaged values is done or not at eath time step.
  real(dp), save, public:: interval_time_vars(1:max_vars)
                              ! 出力時間間隔.
                              ! Interval time of output.
  integer, save, public:: interval_unitsym_vars(1:max_vars)
                              ! 出力時間間隔の単位 (シンボル).
                              ! Units (symbols) of interval time of output.
  real(dp), save, public:: prev_outtime_vars(1:max_vars)
                              ! 前回に出力した時間.
                              ! Time of previous output
  logical, save, public:: tavr_vars(1:max_vars) = .false.
                              ! 時間平均フラグ.
                              ! Flag for time average
  real(dp), save, public:: origin_time_vars(1:max_vars)
                              ! 出力開始時刻.
                              ! Start time of output
  real(dp), save, public:: terminus_time_vars(1:max_vars)
                              ! ファイルをクローズする時刻.
                              ! time of closure of file
  logical, save, public:: histaddvar_vars(1:max_vars) = .false.
                              ! HistoryAddVariable 済みかどうか
                              ! Whether "HistoryAddVariable" is done or not.
  real(dp), save, public:: newfile_inttime_vars(1:max_vars)
                              ! ファイルを新規に作り直す時間間隔.
                              ! Interval time of remake of file
  real(dp), save, public:: newfile_createtime_vars(1:max_vars)
                              ! ファイルを新規に作り直した時間.
                              ! Time of remake of file
  logical, save, public:: flag_output_prev_vars(1:max_vars) = .false.
                              ! ファイル出力を一度でも行ったかどうかのフラグ
                              ! Flag implying that file is output previously
  real(dp), save, public:: zero_time
                              ! ゼロ秒. Zero second
 
  type slice_info
    integer, pointer:: st(:) =>null()
                              ! 空間方向の開始点.
                              ! Start points of spaces.
    integer, pointer:: ed(:) =>null()
                              ! 空間方向の終了点.
                              ! End points of spaces.
    integer, pointer:: sd(:) =>null()
                              ! 空間方向の刻み幅.
                              ! Strides of spaces
  end type slice_info
 
  ! データの切り出し情報
  ! Information of slices of data
  !
  type(slice_info), save, target, public:: slice_vars(1:max_vars)
 
  type axes_weight
                        real(dp), pointer:: wgt1(:) =>null()
                    
    real(dp), pointer:: wgt2(:) =>null()
                    
    real(dp), pointer:: wgt3(:) =>null()
                    
    real(dp), pointer:: wgt4(:) =>null()
                    
    real(dp), pointer:: wgt5(:) =>null()
                    
    real(dp), pointer:: wgt6(:) =>null()
                    
    real(dp), pointer:: wgt7(:) =>null()
                    
  end type axes_weight
 
  ! 座標重み情報
  ! Information of axes weight
  !
  type(axes_weight), save, target, public:: weight_vars(1:max_vars)
 
  type space_avr_info
    logical, pointer:: avr(:) =>null()
                              ! 平均化のフラグ.
                              ! Flag of average.
  end type space_avr_info
 
  ! データの切り出し情報
  ! Information of slices of data
  !
  type(space_avr_info), save, target, public:: space_avr_vars(1:max_vars)
 
  ! 登録変数を全て出力するためのフラグ.
  ! Flag for output all registered variables.
  !
  logical, save, public:: all_output_save = .false.
 
  ! 変数登録は確定されているか.
  ! * HistoryAutoAllVarFix が呼ばれると .true. になる.
  ! * 一度 .true. になると, HistoryAutoAddVariable を呼ぶことはできない.
  !
  ! Whether register of variables is fixed.
  ! * When "HistoryAutoAllVarFix" is called, this argument becomes .true.
  ! * Once this argument becomes .true., "HistoryAutoAddVariable" can not be called.
  !
  logical, save, public:: flag_allvarfixed = .false.
 
  ! 暦情報
  ! Calendar
  !
  type(dc_cal), save, public:: cal_save
 
!!$  ! モデルの開始日時
!!$  ! Start date and time of a model
!!$  !
!!$  type(DC_CAL_DATE), save, public:: start_date_save
 
  ! 初期設定フラグ
  ! Initialization flag
  !
  logical, save, public:: initialized = .false.
 
  character(*), parameter, public:: sub_sname = "HistAuto"
 
  character(*), parameter, public:: version = &
    & '$Name:  $' // &
    & '$Id: gtool_historyauto_internal.rb2f90,v 1.6 2010-07-04 22:01:51 morikawa Exp $'
 
  interface hstvarsoutputcheck
    module procedure hstvarsoutputcheck
  end interface
 
  interface hstfilecreate
    module procedure hstfilecreate
  end interface
 
  interface averagereduce
 
 
    module procedure averagereducereal1
 
 
    module procedure averagereducereal2
 
 
    module procedure averagereducereal3
 
 
    module procedure averagereducereal4
 
 
    module procedure averagereducereal5
 
 
    module procedure averagereducereal6
 
 
    module procedure averagereducereal7
 
 
    module procedure averagereducedouble1
 
 
    module procedure averagereducedouble2
 
 
    module procedure averagereducedouble3
 
 
    module procedure averagereducedouble4
 
 
    module procedure averagereducedouble5
 
 
    module procedure averagereducedouble6
 
 
    module procedure averagereducedouble7
 
 
    module procedure averagereduceint1
 
 
    module procedure averagereduceint2
 
 
    module procedure averagereduceint3
 
 
    module procedure averagereduceint4
 
 
    module procedure averagereduceint5
 
 
    module procedure averagereduceint6
 
 
    module procedure averagereduceint7
 
 
  end interface
 
contains
 
  !-------------------------------------------------------------------
  !------------- 内部サブルーチン ; Internal Subroutines -------------
  !-------------------------------------------------------------------
 
  subroutine hstvarsoutputcheck ( time, stime_index )
    use dc_trace, only: beginsub, endsub, dbgmessage
    use dc_error, only: storeerror, dc_noerr
    use gtool_history, only: historyinitialized, historyclose
    use dc_date_types, only: dc_difftime
    use dc_date, only: operator(==), operator(>), operator(<), &
      & operator(>=), operator(<=), operator(-), dcdifftimeputline, &
      & evalsec
    implicit none
    real(DP), intent(in):: time
                              ! 現在時刻. Current time
    integer, intent(out):: stime_index
 
    integer:: tstep
    integer:: stat, i, startnum, endnum
    character(STRING):: cause_c
    character(*), parameter:: subname = "HstVarsOutputCheck"
    real(DP), parameter:: time_eps = epsilon(1.0_dp)
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = ""
 
    ! 与えられた時刻がチェック済みかどうかを調べる
    ! Examine whether given time is already checked or not
    !
    timestepsearch: do
      do i = saved_tstep, checked_tstepnum
        if ( abs(saved_time(i) - time) <= &
          & time_eps * max(1.0_dp, abs(time), abs(saved_time(i))) ) then
          tstep = i
          exit timestepsearch
        end if
      end do
      do i = 1, saved_tstep - 1
        if ( abs(saved_time(i) - time) <= &
          & time_eps * max(1.0_dp, abs(time), abs(saved_time(i))) ) then
          tstep = i
          exit timestepsearch
        end if
      end do
 
      tstep = 0
      exit timestepsearch
    end do timestepsearch
 
    saved_tstep = tstep
 
    if ( saved_tstep /= 0 .and. checked_tstep_varnum == numvars ) then
      ! * output_timing_vars(:,saved_tstep) を使う.
      ! * saved_tstep を stime_index として返す.
 
      stime_index = saved_tstep
      call dbgmessage( 'saved_tstep=<%d> is already checked.', &
        & i =(/ saved_tstep /) )
      goto 999
    end if
 
    ! チェックする時間ステップと, 変数 ID の設定
    ! Configure checked time step, and variable ID
    !
    if ( saved_tstep /= 0 ) then
      startnum = checked_tstep_varnum + 1
      endnum   = numvars
 
      stime_index = saved_tstep
    else
      startnum = 1
      endnum   = numvars
 
      if ( save_tstepnum < 2 ) then
        checked_tstepnum = 1
        saved_time(checked_tstepnum) = time
 
        saved_tstep = checked_tstepnum
        stime_index = saved_tstep
 
      elseif ( .not. checked_tstepnum < save_tstepnum ) then
        create_timing_vars(:,1:checked_tstepnum-1) = create_timing_vars(:,2:checked_tstepnum)
        close_timing_vars(:,1:checked_tstepnum-1) = close_timing_vars(:,2:checked_tstepnum)
        renew_timing_vars(:,1:checked_tstepnum-1) = renew_timing_vars(:,2:checked_tstepnum)
        output_timing_vars(:,1:checked_tstepnum-1) = output_timing_vars(:,2:checked_tstepnum)
        output_timing_avr_vars(:,1:checked_tstepnum-1) = output_timing_avr_vars(:,2:checked_tstepnum)
 
        saved_time(1:checked_tstepnum-1) = saved_time(2:checked_tstepnum)
        saved_time(checked_tstepnum) = time
 
        saved_tstep = checked_tstepnum
        stime_index = saved_tstep
 
      else
        checked_tstepnum = checked_tstepnum + 1
        saved_time(checked_tstepnum) = time
 
        saved_tstep = checked_tstepnum
        stime_index = saved_tstep
      end if
    end if
 
    call dbgmessage( 'numvar=<%d:%d> in saved_tstep=<%d> will be checked from now.', &
      & i =(/ startnum, endnum, saved_tstep /) )
 
 
    ! それぞれのタイミングをチェックして各変数に格納
    !
    ! * ファイルオープン:      create_timing_vars
    ! * ファイルクローズ:      close_timing_vars
    ! * ファイルクローズ/作成: renew_timing_vars
    ! * データ出力:            output_timing_vars
    ! * データ平均化:          output_avr_timing_vars
 
    create_timing_vars(startnum:endnum, checked_tstepnum:checked_tstepnum) = .false.
    close_timing_vars(startnum:endnum, checked_tstepnum:checked_tstepnum) = .false.
    renew_timing_vars(startnum:endnum, checked_tstepnum:checked_tstepnum) = .false.
    output_timing_vars(startnum:endnum, checked_tstepnum:checked_tstepnum) = .false.
    output_timing_avr_vars(startnum:endnum, checked_tstepnum:checked_tstepnum) = .false.
 
    do i = startnum, endnum
 
      if ( .not. output_valid_vars(i) ) cycle
 
      if ( origin_time_vars(i) > time ) cycle
 
      if (             origin_time_vars(i) <= time &
        &  .and.       (      terminus_time_vars(i) < zero_time &
        &                .or. terminus_time_vars(i) >= time      ) &
        &  .and. .not. histaddvar_vars(i)            ) then
 
        create_timing_vars(i,checked_tstepnum) = .true.
 
        if ( newfile_inttime_vars(i) > zero_time ) then
          newfile_createtime_vars(i) = time
        end if
 
        output_timing_vars(i,checked_tstepnum) = .true.
        output_timing_avr_vars(i,checked_tstepnum) = tavr_vars(i)
        cycle
      end if
 
      if ( terminus_time_vars(i) > zero_time .and. terminus_time_vars(i) < time ) then
        close_timing_vars(i,checked_tstepnum) = .true.
        output_timing_vars(i,checked_tstepnum) = .false.
        output_timing_avr_vars(i,checked_tstepnum) = .false.
        cycle
      end if
 
      ! * newfile_inttime_vars だけでは, ぴったり一致しないとうまく行かない.
      !   * そこで...
      !     * 前回に出力した時刻を記憶しておく.
      !     * 前回の時刻と今回の時刻の差が newfile_inttime_vars
      !       よりも大きい場合には現ファイルを閉じ, 新ファイルを作成する.
 
      if ( newfile_inttime_vars(i) > zero_time ) then
        if ( time - newfile_createtime_vars(i) >= newfile_inttime_vars(i) ) then
          renew_timing_vars(i,checked_tstepnum) = .true.
 
          output_timing_vars(i,checked_tstepnum) = .true.
          output_timing_avr_vars(i,checked_tstepnum) = tavr_vars(i)
 
          cycle
        end if
      end if
 
      if ( time - prev_outtime_vars(i) >= interval_time_vars(i) ) then
        output_timing_vars(i,checked_tstepnum) = .true.
        output_timing_avr_vars(i,checked_tstepnum) = tavr_vars(i)
        cycle
      end if
 
      output_timing_vars(i,checked_tstepnum) = .false.
      output_timing_avr_vars(i,checked_tstepnum) = tavr_vars(i)
 
    end do
 
    checked_tstep_varnum = numvars
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c)
    call endsub(subname)
  end subroutine hstvarsoutputcheck
 
  subroutine hstfilecreate( &
    & gthist,  & ! (inout)
    & varname, & ! (in)
    & time &     ! (in)
    & )
    use dc_trace, only: beginsub, endsub
    use dc_error, only: storeerror, dc_noerr, hst_ebadvarname, &
      & hst_ebadslice, hst_empinoaxisdata
    use dc_calendar, only: dccalconvertbyunit
    use dc_date_types, only: dc_difftime
    use dc_date, only: dcdifftimecreate, evalbyunit
    use dc_string, only: cprintf, strinclude, tochar, joinchar
    use dc_message, only: messagenotify
    use gtool_history_nmlinfo_generic, only: &
      & hstnmlinfooutputvalid, hstnmlinfoinquire, hstnmlinfoputline
    use gtool_history, only: gt_history, &
      & historycreate, historyaddvariable, historyaddattr, &
      & historyinitialized, historyput, historyputaxismpi, &
      & historyaxiscreate, historyaxisinquire, historyaxiscopy, &
      & historyvarinfoinquire, historyvarinfocreate, &
      & historyvarinfocopy, historyvarinfoinitialized, &
      & historyvarinfoclear
 
    implicit none
    type(gt_history), intent(inout):: gthist
                              ! gtool_history モジュール用構造体.
                              ! Derived type for "gtool_history" module
    character(*), intent(in):: varname
                              ! 変数の名前.
                              ! Variable name
    real(DP), intent(in):: time
                              ! 現在時刻. Current time
 
    character(TOKEN):: interval_unit
                              ! データの出力間隔の単位.
                              ! Unit for interval of history data output
    real(DP):: origin_value
                              ! データの出力開始時刻の数値.
                              ! Numerical value for start time of history data output
    character(TOKEN):: origin_unit
                              ! データの出力開始時刻の単位.
                              ! Unit for start time of history data output
 
    real(DP):: origin_sec
    integer:: newfile_intvalue
    real(DP):: newfile_intvalued
                              ! ファイル分割時間間隔.
                              ! Interval of time of separation of a file.
    character(TOKEN):: newfile_intunit
                              ! ファイル分割時間間隔の単位.
                              ! Unit of interval of time of separation of a file.
 
    character(STRING):: file, file_base, file_suffix, file_newfile_time, file_rank
                              ! 出力ファイル名.
                              ! Output file name.
    integer:: stat, i, j, k, vnum, numdims_slice, dim_size, cnt
    character(STRING):: name, units, longname, cause_c, wgt_name
    character(TOKEN):: xtype
    type(gt_history_axis):: gthst_axes_time
    type(gt_history_axis), pointer:: gthst_axes_slices(:) =>null()
    type(gt_history_axis_data), pointer:: data_axes_slices(:) =>null()
    type(gt_history_axis_data), pointer:: data_weights_slices(:) =>null()
    real(DP):: wgt_sum, wgt_sum_s
    logical:: slice_valid
    integer:: slice_start(1:numdims-1)
                              ! 空間方向の開始点.
                              ! Start points of spaces.
    integer:: slice_end(1:numdims-1)
                              ! 空間方向の終了点.
                              ! End points of spaces.
    integer:: slice_stride(1:numdims-1)
                              ! 空間方向の刻み幅.
                              ! Strides of spaces
 
    character(*), parameter:: subname = "HstFileCreate"
  continue
    call beginsub(subname, 'varname=%c', c1 = trim(varname) )
    stat = dc_noerr
    cause_c = ""
 
    ! varname から変数情報の探査
    ! Search information of a variable from "varname"
    !
    vnum = 0
    do i = 1, numvars
      call historyvarinfoinquire( &
        & varinfo = gthst_vars(i), &  ! (in)
        & name = name )               ! (out)
      if ( trim(varname) == trim(name) ) vnum = i
    end do
 
    if ( vnum == 0 ) then
      stat = hst_ebadvarname
      cause_c = varname
      goto 999
    end if
 
    ! 出力が有効かどうかを確認する
    ! Confirm whether the output is effective
    !
    if ( .not. hstnmlinfooutputvalid( gthstnml, varname ) ) then
      goto 999
    end if
 
    ! 出力間隔の単位に応じて時間座標情報の作り直し
    ! Remake time axis information correspond to units of output interval
    !
    call hstnmlinfoinquire( &
      & gthstnml = gthstnml, &                 ! (in)
      & name = varname, &                      ! (in)
      & file = file, &                         ! (out)
      & interval_unit  = interval_unit )       ! (out)
 
    call historyaxiscopy( &
      & gthst_axes_time, &     ! (out)
      & gthst_axes(numdims), & ! (in)
      & units = trim(interval_unit) // ' ' // &
      &         trim(time_unit_suffix) ) ! (in)
 
    ! 空間方向のスライスに対応して, 座標および座標重み情報の作り直し
    ! Remake axes and weights information correspond to spatial slices
    !
    call hstnmlinfoinquire( &
      & gthstnml = gthstnml, &         ! (in)
      & name = varname, &              ! (in)
      & slice_start = slice_start, &   ! (out)
      & slice_end = slice_end, &       ! (out)
      & slice_stride = slice_stride )  ! (out)
 
    ! ファイルが未作成の場合は, まずファイル作成
    ! At first, the file is created if the file is not created yet
    !
    if ( .not. historyinitialized( gthist ) ) then
 
      if (       all( slice_start  == (/ ( 1, i = 1, numdims -1 ) /) ) &
        &  .and. all( slice_end    <  (/ ( 1, i = 1, numdims -1 ) /) ) &
        &  .and. all( slice_stride == (/ ( 1, i = 1, numdims -1 ) /) )  ) then
 
        allocate( gthst_axes_slices(1:numdims) )
        gthst_axes_slices(1:numdims-1)     = gthst_axes(1:numdims-1)
        gthst_axes_slices(numdims:numdims) = gthst_axes_time
 
        data_axes_slices               => data_axes
        data_weights_slices            => data_weights
        slice_valid = .false.
 
      else
        allocate( gthst_axes_slices(1:numdims) )
        allocate( data_axes_slices(1:numdims) )
        allocate( data_weights_slices(1:numdims) )
 
        do i = 1, numdims-1
 
          ! スライス値の有効性をチェック
          ! Check validity of slices
          !
          if ( slice_start(i) < 1 ) then
            stat = hst_ebadslice
            cause_c = cprintf('slice_start=%d', &
              &         i = (/ slice_start(i) /) )
            goto 999
          end if
 
          if ( slice_stride(i) < 1 ) then
            stat = hst_ebadslice
            cause_c = cprintf('slice_stride=%d', &
              &         i = (/ slice_stride(i) /) )
            goto 999
          end if
 
          ! 再生成の必要性をチェック
          ! Check necessity of remaking
          !
          if (       ( slice_start(i)  == 1 ) &
            &  .and. ( slice_end(i)    <  1 ) &
            &  .and. ( slice_stride(i) == 1 )  ) then
 
            call historyaxiscopy( &
              & axis_dest = gthst_axes_slices(i) , &  ! (out)
              & axis_src  = gthst_axes(i) )           ! (in)
 
            data_axes_slices(i) = data_axes(i)
 
            cycle
          end if
 
          ! 座標情報の再生成
          ! Remake information of axis
          !
          call historyaxisinquire( &
            & axis = gthst_axes(i), &  ! (in)
            & name = name, &           ! (out)
            & size = dim_size, &       ! (out)
            & longname = longname, &   ! (out)
            & units = units, &         ! (out)
            & xtype = xtype )          ! (out)
 
          ! 終点のスライス値の補正 ; Correct end points of slices
          if ( slice_end(i) < 1 ) slice_end(i) = dim_size
          if ( slice_end(i) > dim_size ) then
            call messagenotify( 'W', subname, &
              & 'slice options to (%c) are undesirable ' // &
              & '(@slice_end=%d). @slice_end is corrected forcibly to (%d) ', &
              & c1 = trim(name), &
              & i = (/ slice_end(i), dim_size /) )
 
            slice_end(i) = dim_size
          end if
 
          ! スライス値の有効性をチェック ; Check validity of slices
          if ( slice_start(i) > slice_end(i) ) then
            stat = hst_ebadslice
            cause_c = cprintf('slice_start=%d, slice_end=%d', &
              &         i = (/ slice_start(i), slice_end(i) /) )
            goto 999
          end if
 
          numdims_slice = int( ( slice_end(i) - slice_start(i) + 1 ) / slice_stride(i) )
 
          ! スライス値の有効性をチェック ; Check validity of slices
          if ( numdims_slice < 1 ) then
            call messagenotify( 'W', subname, &
              & 'slice options to (%c) are invalid. ' // &
              & '(@slice_start=%d @slice_end=%d @slice_stride=%d)', &
              & c1 = trim(name), &
              & i = (/ slice_start(i), slice_end(i), slice_stride(i) /) )
            stat = hst_ebadslice
            cause_c = cprintf('slice_start=%d, slice_end=%d, slice_stride=%d', &
              &         i = (/ slice_start(i), slice_end(i), slice_stride(i) /) )
            goto 999
          end if
 
          call historyaxiscreate( &
            & axis = gthst_axes_slices(i), &  ! (out)
            & name = name, &                  ! (in)
            & size = numdims_slice, &         ! (in)
            & longname = longname, &          ! (in)
            & units = units, &                ! (in)
            & xtype = xtype )                 ! (in)
 
 
          ! 座標データの再生成
          ! Regenerate data of axis
          !
          allocate( data_axes_slices(i) % a_axis( numdims_slice ) )
          cnt = 1
          do j = slice_start(i), slice_end(i), slice_stride(i)
            data_axes_slices(i) % a_axis( cnt ) = data_axes(i) % a_axis( j )
            cnt = cnt + 1
          end do
 
          ! 座標重みデータの再生成
          ! Remake information of axis data
          !
          do j = 1, numwgts
            call historyvarinfoinquire( &
              & varinfo = gthst_weights(j), & ! (in)
              & name = wgt_name )             ! (out) optional
 
            if ( trim(name) // wgtsuf == trim(wgt_name) ) then
 
              ! 座標重みの計算は結構いい加減...
              ! Calculation about axis weight is irresponsible...
              !
              wgt_sum = sum( data_weights(j) % a_axis )
 
              allocate( data_weights_slices(j) % a_axis( numdims_slice ) )
              cnt = 1
              do k = slice_start(i), slice_end(i), slice_stride(i)
                data_weights_slices(j) % a_axis( cnt ) = data_weights(j) % a_axis( k )
                cnt = cnt + 1
              end do
 
              wgt_sum_s = sum( data_weights_slices(j) % a_axis )
              data_weights_slices(j) % a_axis = data_weights_slices(j) % a_axis * ( wgt_sum / wgt_sum_s )
 
            end if
 
          end do
 
        end do
 
        ! 空間切り出しされていない座標に関する座標重みデータを作成
        ! Make data of axis weight not sliced
        !
        do i = 1, numwgts
          if ( .not. associated( data_weights_slices(i) % a_axis ) ) then
            allocate( data_weights_slices(i) % a_axis( size(data_weights(i) % a_axis ) ) )
            data_weights_slices(i) % a_axis = data_weights(i) % a_axis
          end if
        end do
 
        ! 時刻次元のコピー
        ! Copy time dimension
        !
        gthst_axes_slices(numdims) = gthst_axes_time
 
        slice_valid = .true.
      end if
 
      ! HistoryCreate のための設定値の取得
      ! Get the settings for "HistoryCreate"
      !
      call hstnmlinfoinquire( &
        & gthstnml = gthstnml, &                 ! (in)
        & name = varname, &                      ! (in)
        & file = file, &                         ! (out)
        & origin_value   = origin_value, &       ! (out)
        & origin_unit    = origin_unit, &        ! (out)
        & interval_unit  = interval_unit, &      ! (out)
        & newfile_intvalue = newfile_intvalue, & ! (out)
        & newfile_intunit = newfile_intunit )    ! (out)
 
      ! データ出力時刻の設定
      ! Configure data output time
      !
      origin_sec = &
        & dccalconvertbyunit( &
        &   real( origin_value, dp ), origin_unit, 'sec', cal_save )
 
!!$      ! dc_date モジュール使用時
!!$      !
!!$      call DCDiffTimeCreate( &
!!$        & origin_sec, &           ! (out)
!!$        & origin_value, origin_unit )  ! (in)
 
      if ( newfile_intvalue < 1 ) then
 
        origin_value = dccalconvertbyunit( &
          & origin_sec, 'sec', interval_unit, cal_save )
 
!        origin_value = EvalbyUnit( origin_sec, interval_unit )
      else
 
        origin_value = &
          & dccalconvertbyunit( time, 'sec', interval_unit, cal_save )
 
!        origin_value = EvalbyUnit( time, interval_unit )
      end if
 
      ! ファイル名の設定
      ! Configure file name
      !
      if ( len_trim( file ) - index(file, '.nc', .true.) == 2 ) then
        file_base = file(1:len_trim( file ) - 3)
        file_suffix = '.nc'
      else
        file_base = file
        file_suffix = ''
      end if
      if ( trim(rank_save) == '' ) then
        file_rank = ''
      else
        file_rank = '_rank' // trim( adjustl(rank_save) )
      end if
      if ( newfile_intvalue > 0 ) then
        newfile_intvalued = &
          & dccalconvertbyunit( time, 'sec', newfile_intunit, cal_save )
 
        file_newfile_time = &
          & cprintf( '_time%08d', i = (/ int( newfile_intvalued ) /) )
!          &   i = (/ int( EvalbyUnit( time, newfile_intunit ) ) /) )
      else
        file_newfile_time = ''
      end if
 
      file = trim(file_base) // trim(file_rank) // trim(file_newfile_time) // trim(file_suffix)
 
      ! HistoryCreate によるファイル作成
      ! Files are created by "HistoryCreate"
      !
      call historycreate( &
        & history = gthist, &                                     ! (inout)
        & file = file, title = title_save, &                      ! (in)
        & source = source_save, institution = institution_save, & ! (in)
        & axes = gthst_axes_slices(1:numdims), &                  ! (in)
        & origind = origin_value, &                               ! (in)
        & conventions = conventions_save, &                       ! (in)
        & gt_version = gt_version_save, &                         ! (in)
        & flag_mpi_split = save_mpi_split, &                      ! (in)
        & flag_mpi_gather = save_mpi_gather )                     ! (in)
 
      ! 座標データを出力
      ! Output axes data
      !
      do i = 1, numdims - 1
        call historyaxisinquire( &
          & axis = gthst_axes_slices(i), &  ! (in)
          & name = name )                   ! (out)
        call historyput( &
          & history = gthist, &                    ! (inout) optional
          & varname = name, &                      ! (in)
          & array = data_axes_slices(i) % a_axis ) ! (in)
      end do
 
      ! MPI 用に領域全体の座標データを出力
      ! Output axes data in whole area for MPI
      !
      if ( save_mpi_gather ) then
        do i = 1, numdims - 1
          call historyaxisinquire( &
            & axis = gthst_axes_slices(i), &  ! (in)
            & name = name )                   ! (out)
 
          if ( .not. associated( data_axes_whole(i) % a_axis ) ) then
            call messagenotify('W', subname, &
              & 'data of axis (%c) in whole area is lack. ' // &
              & 'Specify the data by "HistoryAutoPutAxisMPI" explicitly.', &
              & c1 = trim(name) )
            stat = hst_empinoaxisdata
            cause_c = name
          end if
 
          call historyputaxismpi( &
            & history = gthist, &                   ! (inout) optional
            & varname = name, &                     ! (in)
            & array = data_axes_whole(i) % a_axis ) ! (in)
        end do
      end if
 
      ! 割付解除
      ! Deallocation
      !
      if ( slice_valid ) then
        deallocate( gthst_axes_slices )
        deallocate( data_axes_slices )
      else
        deallocate( gthst_axes_slices )
        nullify( data_axes_slices )
      end if
 
      ! 座標重みデータを追加
      ! Add axes weights data
      !
      do i = 1, numwgts
        call historyaddvariable( &
          & history = gthist, &           ! (inout)
          & varinfo = gthst_weights(i) )  ! (in)
        call historyvarinfoinquire( &
          & varinfo = gthst_weights(i), &  ! (in)
          & name = name )                  ! (out)
        call historyput( &
          & history = gthist, &                       ! (inout) optional
          & varname = name, &                         ! (in)
          & array = data_weights_slices(i) % a_axis ) ! (in)
      end do
 
      if ( slice_valid ) then
        deallocate( data_weights_slices )
      else
        nullify( data_weights_slices )
      end if
 
    ! ファイル作成おしまい; Creation of file is finished
    end if
 
 
    ! 変数情報を追加
    ! Add information of variables
    !
    call historyaddvariable( &
      & varinfo = gthst_vars(vnum), &  ! (in)
      & history = gthist )             ! (inout) optional
 
999 continue
    call storeerror(stat, subname, cause_c = cause_c)
    call endsub(subname)
  end subroutine hstfilecreate
 
 
  subroutine averagereducereal1( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real, intent(in), target:: array(:)
    logical, intent(in):: space_average(1)
                        real(DP), intent(in):: weight1(:)
                    
    real, pointer:: array_avr(:) ! (out)
 
    real, pointer:: array_avr_work(:)
 
                        real, pointer:: array_avr_work1(:)
                    
 
    integer:: array_shape(1)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                          
        &  ) )
      array_avr_work1 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1) = array_avr_work1(1) + &
        & array_avr_work(i) * real(weight1(i), kind=kind(array_avr_work1))
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / &
        & real(weight_sum, kind=kind(array_avr_work1))
      array_avr_work  => array_avr_work1
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
 
  end subroutine averagereducereal1
 
 
  subroutine averagereducereal2( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real, intent(in), target:: array(:,:)
    logical, intent(in):: space_average(2)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real, pointer:: array_avr(:,:) ! (out)
 
    real, pointer:: array_avr_work(:,:)
 
                        real, pointer:: array_avr_work1(:,:)
                    
    real, pointer:: array_avr_work2(:,:)
                    
 
    integer:: array_shape(2)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        &  ) )
      array_avr_work1 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:) = array_avr_work1(1,:) + &
        & array_avr_work(i,:) * real(weight1(i), kind=kind(array_avr_work1))
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / &
        & real(weight_sum, kind=kind(array_avr_work1))
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        &  ) )
      array_avr_work2 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1) = array_avr_work2(:,1) + &
        & array_avr_work(:,i) * real(weight2(i), kind=kind(array_avr_work2))
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / &
        & real(weight_sum, kind=kind(array_avr_work2))
      array_avr_work  => array_avr_work2
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
 
  end subroutine averagereducereal2
 
 
  subroutine averagereducereal3( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real, intent(in), target:: array(:,:,:)
    logical, intent(in):: space_average(3)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real, pointer:: array_avr(:,:,:) ! (out)
 
    real, pointer:: array_avr_work(:,:,:)
 
                        real, pointer:: array_avr_work1(:,:,:)
                    
    real, pointer:: array_avr_work2(:,:,:)
                    
    real, pointer:: array_avr_work3(:,:,:)
                    
 
    integer:: array_shape(3)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work1 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:) = array_avr_work1(1,:,:) + &
        & array_avr_work(i,:,:) * real(weight1(i), kind=kind(array_avr_work1))
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / &
        & real(weight_sum, kind=kind(array_avr_work1))
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work2 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:) = array_avr_work2(:,1,:) + &
        & array_avr_work(:,i,:) * real(weight2(i), kind=kind(array_avr_work2))
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / &
        & real(weight_sum, kind=kind(array_avr_work2))
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work3 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1) = array_avr_work3(:,:,1) + &
        & array_avr_work(:,:,i) * real(weight3(i), kind=kind(array_avr_work3))
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / &
        & real(weight_sum, kind=kind(array_avr_work3))
      array_avr_work  => array_avr_work3
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
 
  end subroutine averagereducereal3
 
 
  subroutine averagereducereal4( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real, intent(in), target:: array(:,:,:,:)
    logical, intent(in):: space_average(4)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real, pointer:: array_avr(:,:,:,:) ! (out)
 
    real, pointer:: array_avr_work(:,:,:,:)
 
                        real, pointer:: array_avr_work1(:,:,:,:)
                    
    real, pointer:: array_avr_work2(:,:,:,:)
                    
    real, pointer:: array_avr_work3(:,:,:,:)
                    
    real, pointer:: array_avr_work4(:,:,:,:)
                    
 
    integer:: array_shape(4)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work1 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:) = array_avr_work1(1,:,:,:) + &
        & array_avr_work(i,:,:,:) * real(weight1(i), kind=kind(array_avr_work1))
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / &
        & real(weight_sum, kind=kind(array_avr_work1))
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work2 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:) = array_avr_work2(:,1,:,:) + &
        & array_avr_work(:,i,:,:) * real(weight2(i), kind=kind(array_avr_work2))
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / &
        & real(weight_sum, kind=kind(array_avr_work2))
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work3 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:) = array_avr_work3(:,:,1,:) + &
        & array_avr_work(:,:,i,:) * real(weight3(i), kind=kind(array_avr_work3))
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / &
        & real(weight_sum, kind=kind(array_avr_work3))
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work4 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1) = array_avr_work4(:,:,:,1) + &
        & array_avr_work(:,:,:,i) * real(weight4(i), kind=kind(array_avr_work4))
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = array_avr_work4 / &
        & real(weight_sum, kind=kind(array_avr_work4))
      array_avr_work  => array_avr_work4
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
 
  end subroutine averagereducereal4
 
 
  subroutine averagereducereal5( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real, intent(in), target:: array(:,:,:,:,:)
    logical, intent(in):: space_average(5)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    real, pointer:: array_avr(:,:,:,:,:) ! (out)
 
    real, pointer:: array_avr_work(:,:,:,:,:)
 
                        real, pointer:: array_avr_work1(:,:,:,:,:)
                    
    real, pointer:: array_avr_work2(:,:,:,:,:)
                    
    real, pointer:: array_avr_work3(:,:,:,:,:)
                    
    real, pointer:: array_avr_work4(:,:,:,:,:)
                    
    real, pointer:: array_avr_work5(:,:,:,:,:)
                    
 
    integer:: array_shape(5)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work1 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:) = array_avr_work1(1,:,:,:,:) + &
        & array_avr_work(i,:,:,:,:) * real(weight1(i), kind=kind(array_avr_work1))
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / &
        & real(weight_sum, kind=kind(array_avr_work1))
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work2 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:) = array_avr_work2(:,1,:,:,:) + &
        & array_avr_work(:,i,:,:,:) * real(weight2(i), kind=kind(array_avr_work2))
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / &
        & real(weight_sum, kind=kind(array_avr_work2))
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work3 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:) = array_avr_work3(:,:,1,:,:) + &
        & array_avr_work(:,:,i,:,:) * real(weight3(i), kind=kind(array_avr_work3))
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / &
        & real(weight_sum, kind=kind(array_avr_work3))
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work4 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:) = array_avr_work4(:,:,:,1,:) + &
        & array_avr_work(:,:,:,i,:) * real(weight4(i), kind=kind(array_avr_work4))
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = array_avr_work4 / &
        & real(weight_sum, kind=kind(array_avr_work4))
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work5 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1) = array_avr_work5(:,:,:,:,1) + &
        & array_avr_work(:,:,:,:,i) * real(weight5(i), kind=kind(array_avr_work5))
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = array_avr_work5 / &
        & real(weight_sum, kind=kind(array_avr_work5))
      array_avr_work  => array_avr_work5
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
 
  end subroutine averagereducereal5
 
 
  subroutine averagereducereal6( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & weight6, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real, intent(in), target:: array(:,:,:,:,:,:)
    logical, intent(in):: space_average(6)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    real(DP), intent(in):: weight6(:)
                    
    real, pointer:: array_avr(:,:,:,:,:,:) ! (out)
 
    real, pointer:: array_avr_work(:,:,:,:,:,:)
 
                        real, pointer:: array_avr_work1(:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work2(:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work3(:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work4(:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work5(:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work6(:,:,:,:,:,:)
                    
 
    integer:: array_shape(6)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work1 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:,:) = array_avr_work1(1,:,:,:,:,:) + &
        & array_avr_work(i,:,:,:,:,:) * real(weight1(i), kind=kind(array_avr_work1))
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / &
        & real(weight_sum, kind=kind(array_avr_work1))
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work2 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:,:) = array_avr_work2(:,1,:,:,:,:) + &
        & array_avr_work(:,i,:,:,:,:) * real(weight2(i), kind=kind(array_avr_work2))
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / &
        & real(weight_sum, kind=kind(array_avr_work2))
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work3 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:,:) = array_avr_work3(:,:,1,:,:,:) + &
        & array_avr_work(:,:,i,:,:,:) * real(weight3(i), kind=kind(array_avr_work3))
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / &
        & real(weight_sum, kind=kind(array_avr_work3))
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work4 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:,:) = array_avr_work4(:,:,:,1,:,:) + &
        & array_avr_work(:,:,:,i,:,:) * real(weight4(i), kind=kind(array_avr_work4))
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = array_avr_work4 / &
        & real(weight_sum, kind=kind(array_avr_work4))
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work5 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1,:) = array_avr_work5(:,:,:,:,1,:) + &
        & array_avr_work(:,:,:,:,i,:) * real(weight5(i), kind=kind(array_avr_work5))
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = array_avr_work5 / &
        & real(weight_sum, kind=kind(array_avr_work5))
      array_avr_work  => array_avr_work5
    end if
 
                        
 
    if ( space_average(6) ) then
      dim_size = array_shape(6)
      array_shape(6) = 1
      allocate( array_avr_work6( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work6 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work6(:,:,:,:,:,1) = array_avr_work6(:,:,:,:,:,1) + &
        & array_avr_work(:,:,:,:,:,i) * real(weight6(i), kind=kind(array_avr_work6))
        weight_sum = weight_sum + weight6(i)
      end do
      array_avr_work6 = array_avr_work6 / &
        & real(weight_sum, kind=kind(array_avr_work6))
      array_avr_work  => array_avr_work6
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & , array_shape(6) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
      if ( associated( array_avr_work6 ) ) deallocate( array_avr_work6 )
                    
 
  end subroutine averagereducereal6
 
 
  subroutine averagereducereal7( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & weight6, &  ! (in)
                    
    & weight7, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real, intent(in), target:: array(:,:,:,:,:,:,:)
    logical, intent(in):: space_average(7)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    real(DP), intent(in):: weight6(:)
                    
    real(DP), intent(in):: weight7(:)
                    
    real, pointer:: array_avr(:,:,:,:,:,:,:) ! (out)
 
    real, pointer:: array_avr_work(:,:,:,:,:,:,:)
 
                        real, pointer:: array_avr_work1(:,:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work2(:,:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work3(:,:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work4(:,:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work5(:,:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work6(:,:,:,:,:,:,:)
                    
    real, pointer:: array_avr_work7(:,:,:,:,:,:,:)
                    
 
    integer:: array_shape(7)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work1 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:,:,:) = array_avr_work1(1,:,:,:,:,:,:) + &
        & array_avr_work(i,:,:,:,:,:,:) * real(weight1(i), kind=kind(array_avr_work1))
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / &
        & real(weight_sum, kind=kind(array_avr_work1))
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work2 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:,:,:) = array_avr_work2(:,1,:,:,:,:,:) + &
        & array_avr_work(:,i,:,:,:,:,:) * real(weight2(i), kind=kind(array_avr_work2))
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / &
        & real(weight_sum, kind=kind(array_avr_work2))
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work3 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:,:,:) = array_avr_work3(:,:,1,:,:,:,:) + &
        & array_avr_work(:,:,i,:,:,:,:) * real(weight3(i), kind=kind(array_avr_work3))
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / &
        & real(weight_sum, kind=kind(array_avr_work3))
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work4 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:,:,:) = array_avr_work4(:,:,:,1,:,:,:) + &
        & array_avr_work(:,:,:,i,:,:,:) * real(weight4(i), kind=kind(array_avr_work4))
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = array_avr_work4 / &
        & real(weight_sum, kind=kind(array_avr_work4))
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work5 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1,:,:) = array_avr_work5(:,:,:,:,1,:,:) + &
        & array_avr_work(:,:,:,:,i,:,:) * real(weight5(i), kind=kind(array_avr_work5))
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = array_avr_work5 / &
        & real(weight_sum, kind=kind(array_avr_work5))
      array_avr_work  => array_avr_work5
    end if
 
                        
 
    if ( space_average(6) ) then
      dim_size = array_shape(6)
      array_shape(6) = 1
      allocate( array_avr_work6( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work6 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work6(:,:,:,:,:,1,:) = array_avr_work6(:,:,:,:,:,1,:) + &
        & array_avr_work(:,:,:,:,:,i,:) * real(weight6(i), kind=kind(array_avr_work6))
        weight_sum = weight_sum + weight6(i)
      end do
      array_avr_work6 = array_avr_work6 / &
        & real(weight_sum, kind=kind(array_avr_work6))
      array_avr_work  => array_avr_work6
    end if
 
                        
 
    if ( space_average(7) ) then
      dim_size = array_shape(7)
      array_shape(7) = 1
      allocate( array_avr_work7( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work7 = 0.0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work7(:,:,:,:,:,:,1) = array_avr_work7(:,:,:,:,:,:,1) + &
        & array_avr_work(:,:,:,:,:,:,i) * real(weight7(i), kind=kind(array_avr_work7))
        weight_sum = weight_sum + weight7(i)
      end do
      array_avr_work7 = array_avr_work7 / &
        & real(weight_sum, kind=kind(array_avr_work7))
      array_avr_work  => array_avr_work7
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & , array_shape(6) &
                    
        & , array_shape(7) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
      if ( associated( array_avr_work6 ) ) deallocate( array_avr_work6 )
                    
      if ( associated( array_avr_work7 ) ) deallocate( array_avr_work7 )
                    
 
  end subroutine averagereducereal7
 
 
  subroutine averagereducedouble1( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real(DP), intent(in), target:: array(:)
    logical, intent(in):: space_average(1)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), pointer:: array_avr(:) ! (out)
 
    real(DP), pointer:: array_avr_work(:)
 
                        real(DP), pointer:: array_avr_work1(:)
                    
 
    integer:: array_shape(1)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                          
        &  ) )
      array_avr_work1 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1) = array_avr_work1(1) + array_avr_work(i) * weight1(i)
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / weight_sum
      array_avr_work  => array_avr_work1
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
 
  end subroutine averagereducedouble1
 
 
  subroutine averagereducedouble2( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real(DP), intent(in), target:: array(:,:)
    logical, intent(in):: space_average(2)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), pointer:: array_avr(:,:) ! (out)
 
    real(DP), pointer:: array_avr_work(:,:)
 
                        real(DP), pointer:: array_avr_work1(:,:)
                    
    real(DP), pointer:: array_avr_work2(:,:)
                    
 
    integer:: array_shape(2)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        &  ) )
      array_avr_work1 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:) = array_avr_work1(1,:) + array_avr_work(i,:) * weight1(i)
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / weight_sum
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        &  ) )
      array_avr_work2 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1) = array_avr_work2(:,1) + array_avr_work(:,i) * weight2(i)
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / weight_sum
      array_avr_work  => array_avr_work2
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
 
  end subroutine averagereducedouble2
 
 
  subroutine averagereducedouble3( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real(DP), intent(in), target:: array(:,:,:)
    logical, intent(in):: space_average(3)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), pointer:: array_avr(:,:,:) ! (out)
 
    real(DP), pointer:: array_avr_work(:,:,:)
 
                        real(DP), pointer:: array_avr_work1(:,:,:)
                    
    real(DP), pointer:: array_avr_work2(:,:,:)
                    
    real(DP), pointer:: array_avr_work3(:,:,:)
                    
 
    integer:: array_shape(3)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work1 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:) = array_avr_work1(1,:,:) + array_avr_work(i,:,:) * weight1(i)
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / weight_sum
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work2 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:) = array_avr_work2(:,1,:) + array_avr_work(:,i,:) * weight2(i)
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / weight_sum
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work3 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1) = array_avr_work3(:,:,1) + array_avr_work(:,:,i) * weight3(i)
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / weight_sum
      array_avr_work  => array_avr_work3
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
 
  end subroutine averagereducedouble3
 
 
  subroutine averagereducedouble4( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real(DP), intent(in), target:: array(:,:,:,:)
    logical, intent(in):: space_average(4)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), pointer:: array_avr(:,:,:,:) ! (out)
 
    real(DP), pointer:: array_avr_work(:,:,:,:)
 
                        real(DP), pointer:: array_avr_work1(:,:,:,:)
                    
    real(DP), pointer:: array_avr_work2(:,:,:,:)
                    
    real(DP), pointer:: array_avr_work3(:,:,:,:)
                    
    real(DP), pointer:: array_avr_work4(:,:,:,:)
                    
 
    integer:: array_shape(4)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work1 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:) = array_avr_work1(1,:,:,:) + array_avr_work(i,:,:,:) * weight1(i)
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / weight_sum
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work2 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:) = array_avr_work2(:,1,:,:) + array_avr_work(:,i,:,:) * weight2(i)
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / weight_sum
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work3 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:) = array_avr_work3(:,:,1,:) + array_avr_work(:,:,i,:) * weight3(i)
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / weight_sum
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work4 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1) = array_avr_work4(:,:,:,1) + array_avr_work(:,:,:,i) * weight4(i)
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = array_avr_work4 / weight_sum
      array_avr_work  => array_avr_work4
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
 
  end subroutine averagereducedouble4
 
 
  subroutine averagereducedouble5( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real(DP), intent(in), target:: array(:,:,:,:,:)
    logical, intent(in):: space_average(5)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    real(DP), pointer:: array_avr(:,:,:,:,:) ! (out)
 
    real(DP), pointer:: array_avr_work(:,:,:,:,:)
 
                        real(DP), pointer:: array_avr_work1(:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work2(:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work3(:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work4(:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work5(:,:,:,:,:)
                    
 
    integer:: array_shape(5)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work1 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:) = array_avr_work1(1,:,:,:,:) + array_avr_work(i,:,:,:,:) * weight1(i)
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / weight_sum
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work2 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:) = array_avr_work2(:,1,:,:,:) + array_avr_work(:,i,:,:,:) * weight2(i)
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / weight_sum
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work3 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:) = array_avr_work3(:,:,1,:,:) + array_avr_work(:,:,i,:,:) * weight3(i)
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / weight_sum
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work4 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:) = array_avr_work4(:,:,:,1,:) + array_avr_work(:,:,:,i,:) * weight4(i)
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = array_avr_work4 / weight_sum
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work5 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1) = array_avr_work5(:,:,:,:,1) + array_avr_work(:,:,:,:,i) * weight5(i)
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = array_avr_work5 / weight_sum
      array_avr_work  => array_avr_work5
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
 
  end subroutine averagereducedouble5
 
 
  subroutine averagereducedouble6( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & weight6, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real(DP), intent(in), target:: array(:,:,:,:,:,:)
    logical, intent(in):: space_average(6)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    real(DP), intent(in):: weight6(:)
                    
    real(DP), pointer:: array_avr(:,:,:,:,:,:) ! (out)
 
    real(DP), pointer:: array_avr_work(:,:,:,:,:,:)
 
                        real(DP), pointer:: array_avr_work1(:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work2(:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work3(:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work4(:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work5(:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work6(:,:,:,:,:,:)
                    
 
    integer:: array_shape(6)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work1 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:,:) = array_avr_work1(1,:,:,:,:,:) + array_avr_work(i,:,:,:,:,:) * weight1(i)
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / weight_sum
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work2 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:,:) = array_avr_work2(:,1,:,:,:,:) + array_avr_work(:,i,:,:,:,:) * weight2(i)
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / weight_sum
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work3 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:,:) = array_avr_work3(:,:,1,:,:,:) + array_avr_work(:,:,i,:,:,:) * weight3(i)
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / weight_sum
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work4 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:,:) = array_avr_work4(:,:,:,1,:,:) + array_avr_work(:,:,:,i,:,:) * weight4(i)
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = array_avr_work4 / weight_sum
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work5 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1,:) = array_avr_work5(:,:,:,:,1,:) + array_avr_work(:,:,:,:,i,:) * weight5(i)
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = array_avr_work5 / weight_sum
      array_avr_work  => array_avr_work5
    end if
 
                        
 
    if ( space_average(6) ) then
      dim_size = array_shape(6)
      array_shape(6) = 1
      allocate( array_avr_work6( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work6 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work6(:,:,:,:,:,1) = array_avr_work6(:,:,:,:,:,1) + array_avr_work(:,:,:,:,:,i) * weight6(i)
        weight_sum = weight_sum + weight6(i)
      end do
      array_avr_work6 = array_avr_work6 / weight_sum
      array_avr_work  => array_avr_work6
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & , array_shape(6) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
      if ( associated( array_avr_work6 ) ) deallocate( array_avr_work6 )
                    
 
  end subroutine averagereducedouble6
 
 
  subroutine averagereducedouble7( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & weight6, &  ! (in)
                    
    & weight7, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    real(DP), intent(in), target:: array(:,:,:,:,:,:,:)
    logical, intent(in):: space_average(7)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    real(DP), intent(in):: weight6(:)
                    
    real(DP), intent(in):: weight7(:)
                    
    real(DP), pointer:: array_avr(:,:,:,:,:,:,:) ! (out)
 
    real(DP), pointer:: array_avr_work(:,:,:,:,:,:,:)
 
                        real(DP), pointer:: array_avr_work1(:,:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work2(:,:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work3(:,:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work4(:,:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work5(:,:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work6(:,:,:,:,:,:,:)
                    
    real(DP), pointer:: array_avr_work7(:,:,:,:,:,:,:)
                    
 
    integer:: array_shape(7)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work1 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:,:,:) = array_avr_work1(1,:,:,:,:,:,:) + array_avr_work(i,:,:,:,:,:,:) * weight1(i)
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = array_avr_work1 / weight_sum
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work2 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:,:,:) = array_avr_work2(:,1,:,:,:,:,:) + array_avr_work(:,i,:,:,:,:,:) * weight2(i)
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = array_avr_work2 / weight_sum
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work3 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:,:,:) = array_avr_work3(:,:,1,:,:,:,:) + array_avr_work(:,:,i,:,:,:,:) * weight3(i)
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = array_avr_work3 / weight_sum
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work4 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:,:,:) = array_avr_work4(:,:,:,1,:,:,:) + array_avr_work(:,:,:,i,:,:,:) * weight4(i)
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = array_avr_work4 / weight_sum
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work5 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1,:,:) = array_avr_work5(:,:,:,:,1,:,:) + array_avr_work(:,:,:,:,i,:,:) * weight5(i)
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = array_avr_work5 / weight_sum
      array_avr_work  => array_avr_work5
    end if
 
                        
 
    if ( space_average(6) ) then
      dim_size = array_shape(6)
      array_shape(6) = 1
      allocate( array_avr_work6( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work6 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work6(:,:,:,:,:,1,:) = array_avr_work6(:,:,:,:,:,1,:) + array_avr_work(:,:,:,:,:,i,:) * weight6(i)
        weight_sum = weight_sum + weight6(i)
      end do
      array_avr_work6 = array_avr_work6 / weight_sum
      array_avr_work  => array_avr_work6
    end if
 
                        
 
    if ( space_average(7) ) then
      dim_size = array_shape(7)
      array_shape(7) = 1
      allocate( array_avr_work7( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work7 = 0.0_dp
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work7(:,:,:,:,:,:,1) = array_avr_work7(:,:,:,:,:,:,1) + array_avr_work(:,:,:,:,:,:,i) * weight7(i)
        weight_sum = weight_sum + weight7(i)
      end do
      array_avr_work7 = array_avr_work7 / weight_sum
      array_avr_work  => array_avr_work7
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & , array_shape(6) &
                    
        & , array_shape(7) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
      if ( associated( array_avr_work6 ) ) deallocate( array_avr_work6 )
                    
      if ( associated( array_avr_work7 ) ) deallocate( array_avr_work7 )
                    
 
  end subroutine averagereducedouble7
 
 
  subroutine averagereduceint1( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    integer, intent(in), target:: array(:)
    logical, intent(in):: space_average(1)
                        real(DP), intent(in):: weight1(:)
                    
    integer, pointer:: array_avr(:) ! (out)
 
    integer, pointer:: array_avr_work(:)
 
                        integer, pointer:: array_avr_work1(:)
                    
 
    integer:: array_shape(1)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                          
        &  ) )
      array_avr_work1 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1) = int( &
        & real(array_avr_work1(1), kind=dp) + &
        & real(array_avr_work(i), kind=dp) * weight1(i), &
        & kind=kind(array_avr_work1) )
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = int( &
        & real(array_avr_work1, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work1) )
      array_avr_work  => array_avr_work1
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
 
  end subroutine averagereduceint1
 
 
  subroutine averagereduceint2( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    integer, intent(in), target:: array(:,:)
    logical, intent(in):: space_average(2)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    integer, pointer:: array_avr(:,:) ! (out)
 
    integer, pointer:: array_avr_work(:,:)
 
                        integer, pointer:: array_avr_work1(:,:)
                    
    integer, pointer:: array_avr_work2(:,:)
                    
 
    integer:: array_shape(2)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        &  ) )
      array_avr_work1 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:) = int( &
        & real(array_avr_work1(1,:), kind=dp) + &
        & real(array_avr_work(i,:), kind=dp) * weight1(i), &
        & kind=kind(array_avr_work1) )
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = int( &
        & real(array_avr_work1, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work1) )
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        &  ) )
      array_avr_work2 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1) = int( &
        & real(array_avr_work2(:,1), kind=dp) + &
        & real(array_avr_work(:,i), kind=dp) * weight2(i), &
        & kind=kind(array_avr_work2) )
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = int( &
        & real(array_avr_work2, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work2) )
      array_avr_work  => array_avr_work2
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
 
  end subroutine averagereduceint2
 
 
  subroutine averagereduceint3( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    integer, intent(in), target:: array(:,:,:)
    logical, intent(in):: space_average(3)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    integer, pointer:: array_avr(:,:,:) ! (out)
 
    integer, pointer:: array_avr_work(:,:,:)
 
                        integer, pointer:: array_avr_work1(:,:,:)
                    
    integer, pointer:: array_avr_work2(:,:,:)
                    
    integer, pointer:: array_avr_work3(:,:,:)
                    
 
    integer:: array_shape(3)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work1 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:) = int( &
        & real(array_avr_work1(1,:,:), kind=dp) + &
        & real(array_avr_work(i,:,:), kind=dp) * weight1(i), &
        & kind=kind(array_avr_work1) )
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = int( &
        & real(array_avr_work1, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work1) )
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work2 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:) = int( &
        & real(array_avr_work2(:,1,:), kind=dp) + &
        & real(array_avr_work(:,i,:), kind=dp) * weight2(i), &
        & kind=kind(array_avr_work2) )
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = int( &
        & real(array_avr_work2, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work2) )
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        &  ) )
      array_avr_work3 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1) = int( &
        & real(array_avr_work3(:,:,1), kind=dp) + &
        & real(array_avr_work(:,:,i), kind=dp) * weight3(i), &
        & kind=kind(array_avr_work3) )
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = int( &
        & real(array_avr_work3, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work3) )
      array_avr_work  => array_avr_work3
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
 
  end subroutine averagereduceint3
 
 
  subroutine averagereduceint4( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    integer, intent(in), target:: array(:,:,:,:)
    logical, intent(in):: space_average(4)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    integer, pointer:: array_avr(:,:,:,:) ! (out)
 
    integer, pointer:: array_avr_work(:,:,:,:)
 
                        integer, pointer:: array_avr_work1(:,:,:,:)
                    
    integer, pointer:: array_avr_work2(:,:,:,:)
                    
    integer, pointer:: array_avr_work3(:,:,:,:)
                    
    integer, pointer:: array_avr_work4(:,:,:,:)
                    
 
    integer:: array_shape(4)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work1 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:) = int( &
        & real(array_avr_work1(1,:,:,:), kind=dp) + &
        & real(array_avr_work(i,:,:,:), kind=dp) * weight1(i), &
        & kind=kind(array_avr_work1) )
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = int( &
        & real(array_avr_work1, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work1) )
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work2 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:) = int( &
        & real(array_avr_work2(:,1,:,:), kind=dp) + &
        & real(array_avr_work(:,i,:,:), kind=dp) * weight2(i), &
        & kind=kind(array_avr_work2) )
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = int( &
        & real(array_avr_work2, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work2) )
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work3 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:) = int( &
        & real(array_avr_work3(:,:,1,:), kind=dp) + &
        & real(array_avr_work(:,:,i,:), kind=dp) * weight3(i), &
        & kind=kind(array_avr_work3) )
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = int( &
        & real(array_avr_work3, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work3) )
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        &  ) )
      array_avr_work4 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1) = int( &
        & real(array_avr_work4(:,:,:,1), kind=dp) + &
        & real(array_avr_work(:,:,:,i), kind=dp) * weight4(i), &
        & kind=kind(array_avr_work4) )
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = int( &
        & real(array_avr_work4, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work4) )
      array_avr_work  => array_avr_work4
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
 
  end subroutine averagereduceint4
 
 
  subroutine averagereduceint5( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    integer, intent(in), target:: array(:,:,:,:,:)
    logical, intent(in):: space_average(5)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    integer, pointer:: array_avr(:,:,:,:,:) ! (out)
 
    integer, pointer:: array_avr_work(:,:,:,:,:)
 
                        integer, pointer:: array_avr_work1(:,:,:,:,:)
                    
    integer, pointer:: array_avr_work2(:,:,:,:,:)
                    
    integer, pointer:: array_avr_work3(:,:,:,:,:)
                    
    integer, pointer:: array_avr_work4(:,:,:,:,:)
                    
    integer, pointer:: array_avr_work5(:,:,:,:,:)
                    
 
    integer:: array_shape(5)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work1 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:) = int( &
        & real(array_avr_work1(1,:,:,:,:), kind=dp) + &
        & real(array_avr_work(i,:,:,:,:), kind=dp) * weight1(i), &
        & kind=kind(array_avr_work1) )
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = int( &
        & real(array_avr_work1, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work1) )
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work2 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:) = int( &
        & real(array_avr_work2(:,1,:,:,:), kind=dp) + &
        & real(array_avr_work(:,i,:,:,:), kind=dp) * weight2(i), &
        & kind=kind(array_avr_work2) )
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = int( &
        & real(array_avr_work2, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work2) )
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work3 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:) = int( &
        & real(array_avr_work3(:,:,1,:,:), kind=dp) + &
        & real(array_avr_work(:,:,i,:,:), kind=dp) * weight3(i), &
        & kind=kind(array_avr_work3) )
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = int( &
        & real(array_avr_work3, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work3) )
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work4 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:) = int( &
        & real(array_avr_work4(:,:,:,1,:), kind=dp) + &
        & real(array_avr_work(:,:,:,i,:), kind=dp) * weight4(i), &
        & kind=kind(array_avr_work4) )
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = int( &
        & real(array_avr_work4, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work4) )
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        &  ) )
      array_avr_work5 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1) = int( &
        & real(array_avr_work5(:,:,:,:,1), kind=dp) + &
        & real(array_avr_work(:,:,:,:,i), kind=dp) * weight5(i), &
        & kind=kind(array_avr_work5) )
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = int( &
        & real(array_avr_work5, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work5) )
      array_avr_work  => array_avr_work5
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
 
  end subroutine averagereduceint5
 
 
  subroutine averagereduceint6( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & weight6, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    integer, intent(in), target:: array(:,:,:,:,:,:)
    logical, intent(in):: space_average(6)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    real(DP), intent(in):: weight6(:)
                    
    integer, pointer:: array_avr(:,:,:,:,:,:) ! (out)
 
    integer, pointer:: array_avr_work(:,:,:,:,:,:)
 
                        integer, pointer:: array_avr_work1(:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work2(:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work3(:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work4(:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work5(:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work6(:,:,:,:,:,:)
                    
 
    integer:: array_shape(6)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work1 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:,:) = int( &
        & real(array_avr_work1(1,:,:,:,:,:), kind=dp) + &
        & real(array_avr_work(i,:,:,:,:,:), kind=dp) * weight1(i), &
        & kind=kind(array_avr_work1) )
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = int( &
        & real(array_avr_work1, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work1) )
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work2 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:,:) = int( &
        & real(array_avr_work2(:,1,:,:,:,:), kind=dp) + &
        & real(array_avr_work(:,i,:,:,:,:), kind=dp) * weight2(i), &
        & kind=kind(array_avr_work2) )
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = int( &
        & real(array_avr_work2, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work2) )
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work3 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:,:) = int( &
        & real(array_avr_work3(:,:,1,:,:,:), kind=dp) + &
        & real(array_avr_work(:,:,i,:,:,:), kind=dp) * weight3(i), &
        & kind=kind(array_avr_work3) )
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = int( &
        & real(array_avr_work3, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work3) )
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work4 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:,:) = int( &
        & real(array_avr_work4(:,:,:,1,:,:), kind=dp) + &
        & real(array_avr_work(:,:,:,i,:,:), kind=dp) * weight4(i), &
        & kind=kind(array_avr_work4) )
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = int( &
        & real(array_avr_work4, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work4) )
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work5 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1,:) = int( &
        & real(array_avr_work5(:,:,:,:,1,:), kind=dp) + &
        & real(array_avr_work(:,:,:,:,i,:), kind=dp) * weight5(i), &
        & kind=kind(array_avr_work5) )
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = int( &
        & real(array_avr_work5, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work5) )
      array_avr_work  => array_avr_work5
    end if
 
                        
 
    if ( space_average(6) ) then
      dim_size = array_shape(6)
      array_shape(6) = 1
      allocate( array_avr_work6( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        &  ) )
      array_avr_work6 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work6(:,:,:,:,:,1) = int( &
        & real(array_avr_work6(:,:,:,:,:,1), kind=dp) + &
        & real(array_avr_work(:,:,:,:,:,i), kind=dp) * weight6(i), &
        & kind=kind(array_avr_work6) )
        weight_sum = weight_sum + weight6(i)
      end do
      array_avr_work6 = int( &
        & real(array_avr_work6, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work6) )
      array_avr_work  => array_avr_work6
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & , array_shape(6) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
      if ( associated( array_avr_work6 ) ) deallocate( array_avr_work6 )
                    
 
  end subroutine averagereduceint6
 
 
  subroutine averagereduceint7( &
    & array, space_average, &     ! (in)
                        & weight1, &  ! (in)
                    
    & weight2, &  ! (in)
                    
    & weight3, &  ! (in)
                    
    & weight4, &  ! (in)
                    
    & weight5, &  ! (in)
                    
    & weight6, &  ! (in)
                    
    & weight7, &  ! (in)
                    
    & array_avr &                 ! (out)
    )
    implicit none
    integer, intent(in), target:: array(:,:,:,:,:,:,:)
    logical, intent(in):: space_average(7)
                        real(DP), intent(in):: weight1(:)
                    
    real(DP), intent(in):: weight2(:)
                    
    real(DP), intent(in):: weight3(:)
                    
    real(DP), intent(in):: weight4(:)
                    
    real(DP), intent(in):: weight5(:)
                    
    real(DP), intent(in):: weight6(:)
                    
    real(DP), intent(in):: weight7(:)
                    
    integer, pointer:: array_avr(:,:,:,:,:,:,:) ! (out)
 
    integer, pointer:: array_avr_work(:,:,:,:,:,:,:)
 
                        integer, pointer:: array_avr_work1(:,:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work2(:,:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work3(:,:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work4(:,:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work5(:,:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work6(:,:,:,:,:,:,:)
                    
    integer, pointer:: array_avr_work7(:,:,:,:,:,:,:)
                    
 
    integer:: array_shape(7)
    integer:: i, dim_size
    real(DP):: weight_sum
  continue
 
    array_shape = shape( array )
    array_avr_work => array
 
                    
                      
                        
    if ( space_average(1) ) then
      dim_size = array_shape(1)
      array_shape(1) = 1
      allocate( array_avr_work1( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work1 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work1(1,:,:,:,:,:,:) = int( &
        & real(array_avr_work1(1,:,:,:,:,:,:), kind=dp) + &
        & real(array_avr_work(i,:,:,:,:,:,:), kind=dp) * weight1(i), &
        & kind=kind(array_avr_work1) )
        weight_sum = weight_sum + weight1(i)
      end do
      array_avr_work1 = int( &
        & real(array_avr_work1, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work1) )
      array_avr_work  => array_avr_work1
    end if
 
                        
 
    if ( space_average(2) ) then
      dim_size = array_shape(2)
      array_shape(2) = 1
      allocate( array_avr_work2( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work2 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work2(:,1,:,:,:,:,:) = int( &
        & real(array_avr_work2(:,1,:,:,:,:,:), kind=dp) + &
        & real(array_avr_work(:,i,:,:,:,:,:), kind=dp) * weight2(i), &
        & kind=kind(array_avr_work2) )
        weight_sum = weight_sum + weight2(i)
      end do
      array_avr_work2 = int( &
        & real(array_avr_work2, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work2) )
      array_avr_work  => array_avr_work2
    end if
 
                        
 
    if ( space_average(3) ) then
      dim_size = array_shape(3)
      array_shape(3) = 1
      allocate( array_avr_work3( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work3 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work3(:,:,1,:,:,:,:) = int( &
        & real(array_avr_work3(:,:,1,:,:,:,:), kind=dp) + &
        & real(array_avr_work(:,:,i,:,:,:,:), kind=dp) * weight3(i), &
        & kind=kind(array_avr_work3) )
        weight_sum = weight_sum + weight3(i)
      end do
      array_avr_work3 = int( &
        & real(array_avr_work3, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work3) )
      array_avr_work  => array_avr_work3
    end if
 
                        
 
    if ( space_average(4) ) then
      dim_size = array_shape(4)
      array_shape(4) = 1
      allocate( array_avr_work4( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work4 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work4(:,:,:,1,:,:,:) = int( &
        & real(array_avr_work4(:,:,:,1,:,:,:), kind=dp) + &
        & real(array_avr_work(:,:,:,i,:,:,:), kind=dp) * weight4(i), &
        & kind=kind(array_avr_work4) )
        weight_sum = weight_sum + weight4(i)
      end do
      array_avr_work4 = int( &
        & real(array_avr_work4, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work4) )
      array_avr_work  => array_avr_work4
    end if
 
                        
 
    if ( space_average(5) ) then
      dim_size = array_shape(5)
      array_shape(5) = 1
      allocate( array_avr_work5( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work5 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work5(:,:,:,:,1,:,:) = int( &
        & real(array_avr_work5(:,:,:,:,1,:,:), kind=dp) + &
        & real(array_avr_work(:,:,:,:,i,:,:), kind=dp) * weight5(i), &
        & kind=kind(array_avr_work5) )
        weight_sum = weight_sum + weight5(i)
      end do
      array_avr_work5 = int( &
        & real(array_avr_work5, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work5) )
      array_avr_work  => array_avr_work5
    end if
 
                        
 
    if ( space_average(6) ) then
      dim_size = array_shape(6)
      array_shape(6) = 1
      allocate( array_avr_work6( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work6 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work6(:,:,:,:,:,1,:) = int( &
        & real(array_avr_work6(:,:,:,:,:,1,:), kind=dp) + &
        & real(array_avr_work(:,:,:,:,:,i,:), kind=dp) * weight6(i), &
        & kind=kind(array_avr_work6) )
        weight_sum = weight_sum + weight6(i)
      end do
      array_avr_work6 = int( &
        & real(array_avr_work6, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work6) )
      array_avr_work  => array_avr_work6
    end if
 
                        
 
    if ( space_average(7) ) then
      dim_size = array_shape(7)
      array_shape(7) = 1
      allocate( array_avr_work7( array_shape(1) &
                                  & , array_shape(2) &
                          
        & , array_shape(3) &
                          
        & , array_shape(4) &
                          
        & , array_shape(5) &
                          
        & , array_shape(6) &
                          
        & , array_shape(7) &
                          
        &  ) )
      array_avr_work7 = 0
      weight_sum = 0.0_dp
      do i = 1, dim_size
        array_avr_work7(:,:,:,:,:,:,1) = int( &
        & real(array_avr_work7(:,:,:,:,:,:,1), kind=dp) + &
        & real(array_avr_work(:,:,:,:,:,:,i), kind=dp) * weight7(i), &
        & kind=kind(array_avr_work7) )
        weight_sum = weight_sum + weight7(i)
      end do
      array_avr_work7 = int( &
        & real(array_avr_work7, kind=dp) / weight_sum, &
        & kind=kind(array_avr_work7) )
      array_avr_work  => array_avr_work7
    end if
 
                        
                      
 
                    
 
 
      allocate( array_avr( array_shape(1) &
                            & , array_shape(2) &
                    
        & , array_shape(3) &
                    
        & , array_shape(4) &
                    
        & , array_shape(5) &
                    
        & , array_shape(6) &
                    
        & , array_shape(7) &
                    
        & ) )
 
      array_avr = array_avr_work
 
      nullify( array_avr_work )
 
                          if ( associated( array_avr_work1 ) ) deallocate( array_avr_work1 )
                    
      if ( associated( array_avr_work2 ) ) deallocate( array_avr_work2 )
                    
      if ( associated( array_avr_work3 ) ) deallocate( array_avr_work3 )
                    
      if ( associated( array_avr_work4 ) ) deallocate( array_avr_work4 )
                    
      if ( associated( array_avr_work5 ) ) deallocate( array_avr_work5 )
                    
      if ( associated( array_avr_work6 ) ) deallocate( array_avr_work6 )
                    
      if ( associated( array_avr_work7 ) ) deallocate( array_avr_work7 )
                    
 
  end subroutine averagereduceint7
 
 
end module gtool_historyauto_internal
 
!--
! vi:set readonly sw=4 ts=8:
!
!Local Variables:
!mode: f90
!buffer-read-only: t
!End:
!
!++