historyclose.f90

Go to the documentation of this file.

!> @file historyclose.F
!>
!> @author Yasuhiro MORIKAWA, Eizi TOYODA
!> @copyright Copyright (C) GFD Dennou Club, 2000-2026. All rights reserved. <br/>
!>            License is BSD-2-Clause. See [COPYRIGHT](@ref COPYRIGHT) in detail
!>
!> @en
!> @brief Close gtool4 data output
!> @enden
!>
!> @ja
!> @brief gtool4 データの終了処理
!> @endja
!>
!> @en
!> @brief Close gtool4 data output
!>
!> Performs finalization of data output started by HistoryCreate.
!> When HistoryCreate is used in a program, this subroutine must
!> be called before the program terminates.
!> @param[inout] history History structure (optional)
!> @param[in]    quiet   Suppress messages (optional)
!> @param[out]   err     Error flag (optional)
!> @enden
!>
!> @ja
!> @brief gtool4 データの終了処理
!>
!> HistoryCreate で始まったデータ出力の終了処理をおこなうものです。
!> プログラム内で HistoryCreate を用いた場合、プログラムを終了する
!> 前に必ずこのサブルーチンを呼んで下さい。
!> @param[inout] history 出力ファイル設定構造体 (省略可能)
!> @param[in]    quiet   メッセージ出力抑制 (省略可能)
!> @param[out]   err     エラーフラグ (省略可能)
!> @endja
!>
  subroutine historyclose( history, quiet, err )
    !
    use gtool_history_types, only: gt_history
    use gtool_history_internal, only: default, set_fake_dim_value
    use gtool_history_generic, only: historyvarinfoclear
    use gtdata_generic, only: close, inquire
    use gtdata_types, only: gt_variable
    use dc_message, only: messagenotify
    use dc_url, only: urlsplit
    use dc_present, only: present_and_true
    use dc_error, only: storeerror, dc_noerr, dc_enotinit
    use dc_types, only: string, dp
    use dc_date_types, only: unit_symbol_err
    use dc_trace, only: beginsub, endsub, dbgmessage
    implicit none
    type(gt_history), intent(inout), optional, target:: history
                              ! 出力ファイルの設定に関する情報を
                              ! 格納した構造体
                              ! 
                              ! ここに指定するものは,
                              ! HistoryCreate によって初期設定
                              ! されていなければなりません.
                              ! 
    logical, intent(in), optional:: quiet
                              ! .true. を与えた場合, 
                              ! メッセージ出力が抑制されます. 
                              !
                              ! If ".true." is given, 
                              ! messages are suppressed. 
                              ! 
    logical, intent(out), optional:: err
                              ! 例外処理用フラグ.
                              ! デフォルトでは, この手続き内でエラーが
                              ! 生じた場合, プログラムは強制終了します.
                              ! 引数 *err* が与えられる場合,
                              ! プログラムは強制終了せず, 代わりに
                              ! *err* に .true. が代入されます.
                              !
                              ! Exception handling flag. 
                              ! By default, when error occur in 
                              ! this procedure, the program aborts. 
                              ! If this *err* argument is given, 
                              ! .true. is substituted to *err* and 
                              ! the program does not abort. 
    type(gt_history), pointer:: hst =>null()
    character(STRING):: url, file
    integer:: i, v_size
    integer:: stat
    character(STRING):: cause_c
    character(len = *), parameter:: subname = "HistoryClose"
  continue
    call beginsub(subname)
    stat = dc_noerr
    cause_c = ""
    if (present(history)) then
      hst => history
    else
      hst => default
    endif
    !-----------------------------------------------------------------
    !  初期設定のチェック
    !  Check initialization
    !-----------------------------------------------------------------
    if ( .not. hst % initialized ) then
      stat = dc_enotinit
      cause_c = 'GT_HISTORY'
      goto 999
    end if
    !-----------------------------------------------------------------
    !  メッセージ出力用にファイル名取得
    !  Get filename for output messages
    !-----------------------------------------------------------------
    if ( .not. hst % mpi_gather &
      &  .or. ( hst % mpi_gather .and. hst % mpi_myrank == 0 ) ) then
    if ( .not. present_and_true(quiet) ) then
      call inquire( hst % dimvars(1), & ! (in)
        & url = url )                   ! (out)
      call urlsplit( fullname = url, &  ! (in)
        & file = file )                 ! (out)
    end if
    endif
    !-----------------------------------------------------------------
    !  変数のクローズ
    !  Close variables
    !-----------------------------------------------------------------
    if ( .not. hst % mpi_gather &
      &  .or. ( hst % mpi_gather .and. hst % mpi_myrank == 0 ) ) then
    v_size = size(hst % dimvars)
    do, i = 1, v_size
      if (.not. hst % dim_value_written(i)) &
        call set_fake_dim_value(hst, i)
      call close(hst % dimvars(i))
    enddo
    v_size = size(hst % vars)
    do, i = 1, v_size
      call close(hst % vars(i))
    enddo
    endif
    deallocate(hst % dimvars)
    v_size = size(hst % vars)
    !-----------------------------------------------------------------
    !  設定のクリア
    !  Clear configurations
    !-----------------------------------------------------------------
    hst % unlimited_index = 0
    hst % unlimited_units = ''
    hst % unlimited_units_symbol = unit_symbol_err
    if (associated(hst % dim_value_written)) deallocate(hst % dim_value_written)
    if (associated(hst % vars)) deallocate(hst % vars)
    if (associated(hst % growable_indices)) deallocate(hst % growable_indices)
    if (associated(hst % count)) deallocate(hst % count)
    if (associated(hst % var_avr_count)) deallocate(hst % var_avr_count)
    do, i = 1, v_size
      if (associated(hst % var_avr_data(i) % a_DataAvr)) deallocate(hst % var_avr_data(i) % a_DataAvr)
    enddo
    if (associated(hst % var_avr_data)) deallocate(hst % var_avr_data)
    if (associated(hst % var_avr_firstput)) deallocate(hst % var_avr_firstput)
    if (associated(hst % var_avr_coefsum)) deallocate(hst % var_avr_coefsum)
    if (associated(hst % var_avr_baseint)) deallocate(hst % var_avr_baseint)
    if (associated(hst % var_avr_prevtime)) deallocate(hst % var_avr_prevtime)
    hst % time_bnds = 0.0_dp
    hst % time_bnds_output_count = 0
    hst % time_nv_index = 0
    hst % origin_setting = .false.
    if ( associated( hst % mpi_fileinfo % axes ) ) deallocate( hst % mpi_fileinfo % axes )
    if ( associated( hst % mpi_fileinfo ) ) deallocate( hst % mpi_fileinfo )
    v_size = size(hst % mpi_dimdata_all)
    do, i = 1, v_size
      if ( associated( hst % mpi_dimdata_all(i) % a_Axis ) ) deallocate( hst % mpi_dimdata_all(i) % a_Axis )
      if ( associated( hst % mpi_dimdata_all(i) % attrs ) ) deallocate( hst % mpi_dimdata_all(i) % attrs )
    enddo
    if ( associated( hst % mpi_dimdata_all ) ) deallocate( hst % mpi_dimdata_all )
    v_size = size(hst % mpi_dimdata_each)
    do, i = 1, v_size
      if ( associated( hst % mpi_dimdata_each(i) % a_Axis ) ) deallocate( hst % mpi_dimdata_each(i) % a_Axis )
    enddo
    if ( associated( hst % mpi_dimdata_each ) ) deallocate( hst % mpi_dimdata_each )
    if ( associated( hst % mpi_gthr_info ) ) then
      v_size = size(hst % mpi_gthr_info)
      do, i = 1, v_size
        if ( associated( hst % mpi_gthr_info(i) % index_all ) ) deallocate( hst % mpi_gthr_info(i) % index_all )
        if ( associated( hst % mpi_gthr_info(i) % length ) )    deallocate( hst % mpi_gthr_info(i) % length )
      end do
      deallocate( hst % mpi_gthr_info )
    end if
    if ( associated(hst % mpi_varinfo) ) then
      v_size = size(hst % mpi_varinfo)
      do, i = 1, v_size
        call historyvarinfoclear( hst % mpi_varinfo(i), err )
      end do
      deallocate( hst % mpi_varinfo )
    end if
    if ( associated(hst % mpi_vars_index) ) then
      v_size = size(hst % mpi_vars_index)
      do, i = 1, v_size
        if ( associated( hst % mpi_vars_index(i) % each2all ) ) deallocate( hst % mpi_vars_index(i) % each2all )
        if ( associated( hst % mpi_vars_index(i) % allcount ) ) deallocate( hst % mpi_vars_index(i) % allcount )
      end do
      deallocate( hst % mpi_vars_index )
    end if
    !-----------------------------------------------------------------
    !  メッセージ出力
    !  Output messages
    !-----------------------------------------------------------------
    if ( .not. hst % mpi_gather &
      &  .or. ( hst % mpi_gather .and. hst % mpi_myrank == 0 ) ) then
    if ( .not. present_and_true(quiet) ) then
      call messagenotify('M', subname, &
        & '"%c" is closed', &
        & c1 = trim( file ), rank_mpi = -1 )
    end if
    endif
    hst % mpi_myrank     = -1
    hst % mpi_nprocs = -1
    hst % mpi_gather   = .false.
    hst % mpi_split    = .false.
    !-----------------------------------------------------------------
    !  終了処理, 例外処理
    !  Termination and Exception handling
    !-----------------------------------------------------------------
    hst % initialized = .false.
999 continue
    call storeerror( stat, subname, err, cause_c )
    call endsub( subname )
  end subroutine historyclose