!= 物理定数設定 ! != Physical constants settings ! ! Authors:: Yasuhiro MORIKAWA ! Version:: $Id: gridset.f90,v 1.5 2007/06/09 13:12:30 sugiyama Exp $ ! Tag Name:: $Name: $ ! Copyright:: Copyright (C) GFD Dennou Club, 2008. All rights reserved. ! License:: See COPYRIGHT[link:../../../COPYRIGHT] ! module constants ! != 物理定数設定 ! != Physical constants settings ! ! Note that Japanese and English are described in parallel. ! ! 物理定数の設定および保管を行います. ! デフォルト値は地球大気を想定した値が設定されています. ! これらの値は NAMELIST 変数群 constants_nml によって変更することが ! 可能です. ! ! Physical constants are set and stored. ! By default, values on atmosphere of earth are set. ! These values can be changed by NAMELIST group name ! "constants_nml". ! !== Procedures List ! ! ConstantsInit :: 物理定数の設定 ! ------------ :: ------------ ! ConstantsInit :: Settings of physical constants ! !== NAMELIST ! ! NAMELIST#constants_nml ! ! モジュール引用 ; USE statements ! ! 種別型パラメタ ! Kind type parameter ! use dc_types, only: DP ! 倍精度実数型. Double precision. ! 宣言文 ; Declaration statements ! implicit none private ! 公開手続き ! Public procedure ! public:: ConstantsInit ! 公開変数 ! Public variables ! logical, save, public:: constants_inited = .false. ! 初期設定フラグ. ! Initialization flag real(DP), parameter, public:: PI = 3.1415926535897932_DP ! $ \pi $ . ! 円周率. Circular constant real(DP), parameter, public:: GasRUniv = 8.314_DP ! $ R^{*} $ [J K-1 mol-1]. ! 普遍気体定数. Universal gas constant real(DP), parameter, public:: StB = 5.67e-8_DP ! $ \sigma_{SB} $ . ! ステファンボルツマン定数. ! Stefan-Boltzmann constant real(DP), parameter, public:: FKarm = 0.4_DP ! $ k $ . ! カルマン定数. ! Karman constant real(DP), save, public:: RPlanet ! $ a $ [m]. ! 惑星半径. ! Radius of planet real(DP), save, public:: Omega ! $ \Omega $ [s-1]. ! 回転角速度. ! Angular velocity real(DP), save, public:: Grav ! $ g $ [m s-2]. ! 重力加速度. ! Gravitational acceleration real(DP), save, public:: MolWtDry ! $ ?? $ [kg mol-1]. ! 乾燥大気の平均分子量. ! Mean molecular weight of dry air real(DP), save, public:: CpDry ! $ C_p $ [J kg-1 K-1]. ! 乾燥大気の定圧比熱. ! Specific heat of air at constant pressure real(DP), save, public:: GasRDry ! $ R $ [J kg-1 K-1]. ! 乾燥大気の気体定数. ! Gas constant of air real(DP), save, public:: MolWtWet ! $ ?? $ [kg mol-1]. ! 凝結成分の平均分子量. ! Mean molecular weight of condensible elements real(DP), save, public:: CpWet ! $ C_v $ [J kg-1 K-1] . ! 凝結成分の定圧比熱. ! Specific heat of condensible elements at constant pressure real(DP), save, public:: GasRWet ! $ R_v $ [J kg-1 K-1]. ! 凝結成分の気体定数. ! Gas constant of condensible elements real(DP), save, public:: LatentHeat ! $ L $ [J kg-1] . ! 凝結の潜熱. ! Latent heat of condensation real(DP), save, public:: EpsV ! $ \epsilon_v $ . ! 水蒸気分子量比. ! Molecular weight of water vapor ! 非公開変数 ! Private variables ! character(*), parameter:: module_name = 'constants' ! モジュールの名称. ! Module name character(*), parameter:: version = & & '$Name: $' // & & '$Id: $' ! モジュールのバージョン ! Module version ! INTERFACE 文 ; INTERFACE statements ! interface ConstantsInit module procedure ConstantsInit end interface contains subroutine ConstantsInit ! ! constants モジュールの初期化を行います. ! NAMELIST#constants_nml の読み込みはこの手続きで行われます. ! ! "constants" module is initialized. ! NAMELIST#constants_nml is loaded in this procedure. ! ! モジュール引用 ; USE statements ! ! NAMELIST ファイル入力に関するユーティリティ ! Utilities for NAMELIST file input ! use namelist_util, only: namelist_filename, NmlutilMsg ! ファイル入出力補助 ! File I/O support ! use dc_iounit, only: FileOpen ! 種別型パラメタ ! Kind type parameter ! use dc_types, only: STDOUT ! 標準出力の装置番号. Unit number of standard output ! メッセージ出力 ! Message output ! use dc_message, only: MessageNotify ! 宣言文 ; Declaration statements ! implicit none integer:: unit_nml ! NAMELIST ファイルオープン用装置番号. ! Unit number for NAMELIST file open integer:: iostat_nml ! NAMELIST 読み込み時の IOSTAT. ! IOSTAT of NAMELIST read ! NAMELIST 変数群 ! NAMELIST group name ! namelist /constants_nml/ & & RPlanet, Omega, Grav, & & MolWtDry, CpDry, GasRDry, & & MolWtWet, CpWet, GasRWet, & & LatentHeat, EpsV ! ! デフォルト値については初期化手続 "constants#ConstantsInit" ! のソースコードを参照のこと. ! ! Refer to source codes in the initialization procedure ! "constants#ConstantsInit" for the default values. ! character(*), parameter:: subname = 'ConstantsInit' ! 実行文 ; Executable statement ! if ( constants_inited ) return call InitCheck ! デフォルト値の設定 ! Default values settings ! RPlanet = 6.371e6_DP Omega = 2.0_DP * PI / ( 60.0_DP * 60.0_DP * 23.9345_DP ) Grav = 9.8_DP CpDry = 1004.6_DP MolWtDry = 28.964e-3_DP GasRDry = GasRUniv / MolWtDry CpWet = 1810.0_DP MolWtWet = 18.01528e-3_DP GasRWet = GasRUniv / MolWtWet LatentHeat = 2.5e6_DP EpsV = MolWtWet / MolWtDry ! NAMELIST からの入力 ! Input from NAMELIST ! if ( trim(namelist_filename) /= '' ) then call FileOpen( unit_nml, & ! (out) & namelist_filename, mode = 'r' ) ! (in) rewind( unit_nml ) read( unit_nml, & ! (in) & nml = constants_nml, & ! (out) & iostat = iostat_nml ) ! (out) close( unit_nml ) call NmlutilMsg( iostat_nml, module_name ) ! (in) if ( iostat_nml == 0 ) write( STDOUT, nml = constants_nml ) end if ! 印字 ; Print ! call MessageNotify( 'M', module_name, '----- Initialization Messages -----' ) call MessageNotify( 'M', module_name, ' PI = %f', d = (/ PI /) ) call MessageNotify( 'M', module_name, ' GasRUniv = %f', d = (/ GasRUniv /) ) call MessageNotify( 'M', module_name, ' StB = %f', d = (/ StB /) ) call MessageNotify( 'M', module_name, ' FKarm = %f', d = (/ FKarm /) ) call MessageNotify( 'M', module_name, ' RPlanet = %f', d = (/ RPlanet /) ) call MessageNotify( 'M', module_name, ' Omega = %f', d = (/ Omega /) ) call MessageNotify( 'M', module_name, ' Grav = %f', d = (/ Grav /) ) call MessageNotify( 'M', module_name, ' CpDry = %f', d = (/ CpDry /) ) call MessageNotify( 'M', module_name, ' MolWtDry = %f', d = (/ MolWtDry /) ) call MessageNotify( 'M', module_name, ' GasRDry = %f', d = (/ GasRDry /) ) call MessageNotify( 'M', module_name, ' CpWet = %f', d = (/ CpWet /) ) call MessageNotify( 'M', module_name, ' MolWtWet = %f', d = (/ MolWtWet /) ) call MessageNotify( 'M', module_name, ' GasRWet = %f', d = (/ GasRWet /) ) call MessageNotify( 'M', module_name, ' LatentHeat = %f', d = (/ LatentHeat /) ) call MessageNotify( 'M', module_name, ' EpsV = %f', d = (/ EpsV /) ) call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) ) constants_inited = .true. end subroutine ConstantsInit subroutine InitCheck ! ! 依存モジュールの初期化チェック ! ! Check initialization of dependency modules ! モジュール引用 ; USE statements ! ! NAMELIST ファイル入力に関するユーティリティ ! Utilities for NAMELIST file input ! use namelist_util, only: namelist_util_inited ! メッセージ出力 ! Message output ! use dc_message, only: MessageNotify ! 実行文 ; Executable statement ! if ( .not. namelist_util_inited ) & & call MessageNotify( 'E', module_name, '"namelist_util" module is not initialized.' ) end subroutine InitCheck end module constants