historyautochkoutput.f90

Go to the documentation of this file.

!= 変数の出力設定の確認手続
!= Check output setting of a variable
!
! Authors::   Yasuhiro MORIKAWA
! Version::   $Id: historyautochkoutput.f90,v 1.2 2010-12-28 09:53:09 morikawa Exp $
! Tag Name::  $Name:  $
! Copyright:: Copyright (C) GFD Dennou Club, 2010-. All rights reserved. 
! License::   See COPYRIGHT[link:../../../COPYRIGHT]
!
  function historyautochkoutput( varname ) result(result)
    !
    ! 変数の出力設定の確認を行います. 
    !
    ! *varname* に指定された変数名が出力されるよう設定されている
    ! 場合には *true* が, 出力されないよう設定されている場合は
    ! *false* が返ります. 
    ! 
    ! *HistoryAutoCreate* による初期設定がなされていない場合や, 
    ! *varname* に指定された変数名が HistoryAutoAddVariable によって
    ! 登録されていない場合, 返り値に *false* が返ります. 
    !
    ! Check output setting of a variable. 
    !
    ! If output setting of *varname* is valid, *true* is returned. 
    ! If output setting of *varname* is invalid, *false* is returned. 
    !
    ! If initialization with "HistoryAutoCreate" is not done yet, *false* is returned. 
    ! If *varname* is invalid, *false* is returned. 
    !
    use gtool_historyauto_internal, only: initialized, numdims, numvars, &
      & varname_vars, output_valid_vars
    use dc_types, only: dp, string, token
 
    implicit none
    logical:: result
    character(*), intent(in):: varname
                              ! 変数の名前. 
                              ! Variable name
 
    integer:: i, vnum
    integer, save:: svnum = 1
    character(*), parameter:: subname = "HistoryAutoChkOutput"
  continue
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      result = .false.
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      result = .false.
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 出力設定の確認
    ! Check output setting
    !
    result = output_valid_vars( vnum )
 
999 continue
  end function historyautochkoutput
 
 
 
  function historyautochkoutputtiming( time, varname ) result(result)
    !
    ! 変数の出力設定の確認を行います. 
    !
    ! *varname* に指定された変数名が、*time* のタイミングで出力されるよう
    ! 設定されている場合には *true* が, 出力されないよう設定されている場合は
    ! *false* が返ります. 
    !
    ! なお, gtool_historyauto は時間ステップ可変に対応しているため, 
    ! *time* のタイミングで出力されるかどうかについては, 
    ! その前回に出力された時間に依存します. 従って, とある時間の
    ! 出力設定の確認は下記のように HistoryAutoPut の前で使用して下さい. 
    !
    !   if ( HistoryAutoChkOutputTiming( time, var ) ) then
    !     <some operation ...>
    !   end if 
    !
    !   call HistoryAutoPut( time, var, data )
    !
    ! 以下のように使用した場合には, 期待するような返り値が
    ! 得られないことにご留意下さい. 
    ! 
    !   do i = 1, 10 
    !     write(*,*) HistoryAutoChkOutputTiming( i * timestep, var )
    !   end do
    !
    ! *HistoryAutoCreate* による初期設定がなされていない場合や, 
    ! *varname* に指定された変数名が HistoryAutoAddVariable によって
    ! 登録されていない場合, 返り値に *false* が返ります. 
    !
    ! Check output setting of a variable. 
    !
    ! If *varname* is output on *time*, *true* is returned. 
    !
    ! Please use this function as follows. 
    ! 
    !   if ( HistoryAutoChkOutputTiming( time, var ) ) then
    !     <some operation ...>
    !   end if 
    !
    !   call HistoryAutoPut( time, var, data )
    !
    ! Following usage does not return correct values.
    ! 
    !   do i = 1, 10 
    !     write(*,*) HistoryAutoChkOutputTiming( i * timestep, var )
    !   end do
    !
    ! If initialization with "HistoryAutoCreate" is not done yet, *false* is returned. 
    ! If *varname* is invalid, *false* is returned. 
    !
    use gtool_historyauto_internal, only: initialized, numdims, numvars, &
      & varname_vars, output_valid_vars, origin_time_vars, histaddvar_vars, &
      & newfile_inttime_vars, newfile_createtime_vars, prev_outtime_vars, &
      & interval_time_vars, terminus_time_vars
    use dc_types, only: dp, string, token
 
    implicit none
    logical:: result
    real(dp), intent(in):: time
                              ! データの時刻. 
                              ! Time of data
    character(*), intent(in):: varname
                              ! 変数の名前. 
                              ! Variable name
 
    real(dp), parameter:: zero_time = 0.0_dp
    integer:: i, vnum
    integer, save:: svnum = 1
    character(*), parameter:: subname = "HistoryAutoChkOutputTiming"
  continue
 
    ! 初期設定チェック
    ! Check initialization
    !
    if ( .not. initialized ) then
      result = .false.
      goto 999
    end if
 
    ! 変数 ID のサーチ
    ! Search variable ID
    !
    varsearch: do
      do i = svnum, numvars
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
      do i = 1, svnum - 1
        if ( trim( varname_vars(i) ) == trim(varname) ) then
          vnum = i
          exit varsearch
        end if
      end do
 
      result = .false.
      goto 999
    end do varsearch
 
    svnum = vnum
 
    ! 出力設定の確認
    ! Check output setting
    !
    result = output_valid_vars( vnum )
 
    if ( .not. result ) goto 999
 
    if ( origin_time_vars(vnum) > time ) then
      result = .false. 
      goto 999
    end if
 
    if (             origin_time_vars(vnum) <= time &
      &  .and.       (      terminus_time_vars(vnum) < zero_time &
      &                .or. terminus_time_vars(vnum) >= time      ) &
      &  .and. .not. histaddvar_vars(vnum)            ) then
 
      result = .true.
      goto 999
    end if
 
    if ( terminus_time_vars(vnum) > zero_time .and. terminus_time_vars(vnum) < time ) then
      result = .false.
      goto 999
    end if
 
    if ( newfile_inttime_vars(vnum) > zero_time ) then
      if ( time - newfile_createtime_vars(vnum) >= newfile_inttime_vars(vnum) ) then
        result = .true.
        goto 999
      end if
    end if
 
    if ( time - prev_outtime_vars(vnum) >= interval_time_vars(vnum) ) then
      result = .true.
      goto 999
    end if
 
    result = .false.
 
999 continue
  end function historyautochkoutputtiming