dcdatetimediv.f90

Go to the documentation of this file.

! -*- mode: f90; coding: utf-8 -*-
!-----------------------------------------------------------------------
! Copyright (c) 2000-2026 Gtool Development Group. All rights reserved.
!-----------------------------------------------------------------------
 
function dcdatetime_div_fi(diff, denominator) result(result)
  use dc_date_types, only: cyclic_mdays
  use dc_scaledsec, only: dc_scaled_sec, &
    & operator(<), operator(>), operator(<=), operator(>=), &
    & operator(+), operator(-), operator(*), operator(/), &
    & modscl => mod, modulo, int, abs, sign
  use dc_date_types, only: dc_datetime, dc_difftime
  implicit none
  type(dc_difftime):: result
  type(dc_difftime), intent(in):: diff
  integer, intent(in):: denominator
continue
  result % mon = int( diff % mon / denominator )
  ! 月からの近似的繰り下がりは日単位でしか行わない
  result % day = &
    &   int( diff % day / denominator ) &
    & + int( (cyclic_mdays * modscl(diff % mon, denominator)) / denominator )
  result % sec = diff % sec / denominator + &
    & (diff % day_seconds * modscl(diff % day, denominator)) / &
    & denominator
  result % nondim_flag = diff % nondim_flag
end function dcdatetime_div_fi
 
function dcdatetime_div_fr(diff, denominator) result(result)
  use dc_date_generic, only: operator(/)
  use dc_date_types, only: dc_datetime, dc_difftime
  use dc_types, only: dp
  implicit none
  type(dc_difftime):: result
  type(dc_difftime), intent(in):: diff
  real, intent(in):: denominator
continue
  result = diff / real(denominator, dp)
end function dcdatetime_div_fr
 
function dcdatetime_div_fd(diff, denominator) result(result)
  use dc_date_internal, only: dcdate_normalize
  use dc_date_types, only: cyclic_mdays, dc_datetime, dc_difftime
  use dc_scaledsec, only: dc_scaled_sec, &
    & operator(<), operator(>), operator(<=), operator(>=), &
    & operator(+), operator(-), operator(*), operator(/), &
    & modulo, int, abs, sign
  use dc_types, only: dp
  implicit none
  type(dc_difftime):: result
  type(dc_difftime), intent(in):: diff
  real(dp), intent(in):: denominator
  type(dc_scaled_sec):: month, day
continue
  month = int( diff % mon / denominator )
  result % mon = int(month)
  day =   int( diff % day / denominator ) &
    &   + int(cyclic_mdays * (month - result % mon))
  result % day = int(day)
  result % sec = &
    & diff % sec / denominator + (day - result % day) * diff % day_seconds
  result % day_seconds = diff % day_seconds
  result % nondim_flag = diff % nondim_flag
  call dcdate_normalize(result % day, result % sec, result % day_seconds, result % nondim_flag)
end function dcdatetime_div_fd
 
function dcdatetime_div_ff(diff1, diff2) result(result)
  use dc_date_types, only: cyclic_mdays
  use dc_scaledsec, only: dc_scaled_sec, assignment(=), &
    & operator(<), operator(>), operator(<=), operator(>=), &
    & operator(+), operator(-), operator(*), operator(/), &
    & modulo, int, abs, sign
  use dc_date_types, only: dc_datetime, dc_difftime
  use dc_types, only: dp
  implicit none
  real(dp):: result
  type(dc_difftime), intent(in):: diff1, diff2
continue
  ! ゼロ割対応コードが必要か?
  result = &
    & (diff1 % day_seconds * (cyclic_mdays * diff1 % mon + diff1 % day) &
    & + diff1 % sec) / &
    & (diff2 % day_seconds * (cyclic_mdays * diff2 % mon + diff2 % day) &
    & + diff2 % sec)
end function dcdatetime_div_ff