dc_units Module Reference

単位系処理用モジュール More...

Data Types
type	units
interface	clear
interface	deallocate
interface	operator(*)
interface	operator(/)
interface	operator(+)

Functions/Subroutines
logical function, public	add_okay (u1, u2)

Detailed Description

単位系処理用モジュール

Author

Youhei SASAKI, Eizi TOYODA, Yasuhiro MORIKAWA

Copyright

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

このモジュールは単位系処理のためのプログラム群を提供します. 物理単位を表現する派生型 UNITS と、単位の算術演算 (乗算、除算、加算) のための演算子を含みます.

演算子・手続一覧

演算子/手続	説明
assignment(=)	文字列から UNITS を構築または文字列へ変換
operator(*)	2つの UNITS を乗算
operator(/)	2つの UNITS を除算
operator(+)	2つの UNITS を加算 (互換性チェック付き)
clear	UNITS 変数の初期化
deallocate	UNITS 変数の解放
add_okay	2つの UNITS が加算互換かチェック

Function/Subroutine Documentation

add_okay()

logical function, public dc_units::add_okay	(	type(units), intent(in)	u1,
		type(units), intent(in)	u2 )

2つの UNITS が加算互換かチェック

2つの単位は同じ次元を持つ場合（すなわち、同じ単位要素と 同じべき乗を持つ場合）に互換となります.

Parameters

[in]	u1	1番目の UNITS 変数
[in]	u2	2番目の UNITS 変数

Returns

互換なら .true., そうでなければ .false.

Definition at line 307 of file dc_units.f90.

    type(UNITS), intent(in):: u1, u2
    type(UNITS):: x
    character(STRING):: debug
    call clear(x)
    x = u1 / u2
    debug = u1
    debug = u2
    debug = x
    if (x%nelems == 0) then
      result = .true.
    else if (all(abs(x%power(1:x%nelems)) < tiny(0.0_dp))) then
      result = .true.
    else
      result = .false.
    endif
    call deallocate(x)