Class ChemCalc_3d
In: moist/chemcalc_3d.f90

化学関連の諸量を計算するためのモジュール. AMP と Antoine の飽和蒸気圧式を用いて以下を求める. デフォルトでは AMP 式を使うようにしてある.

 * 飽和蒸気圧
 * 飽和蒸気圧の温度微分
 * 潜熱

Methods

Included Modules

dc_types dc_message ChemData gridset_3d

Public Instance methods

Subroutine :

初期化ルーチン

[Source]

  subroutine ChemCalc_Init( )
    !
    !初期化ルーチン
    !

    !暗黙の型宣言禁止
    implicit none

    !入出力変数
    character(20)      :: Name
    integer            :: id

    !-----------------------------------------------------------
    ! 初期化
    !-----------------------------------------------------------
    call MessageNotify( "M", "ChemCalcSpc_Init", "SpcNum = %d", i=(/SpcNum/) )

    !Antoine の飽和蒸気圧式の係数
    a_antA = ChemData_SvapPress_AntoineA
    a_antB = ChemData_SvapPress_AntoineB
    a_antC = ChemData_SvapPress_AntoineC
    a_antU = ChemData_SvapPress_AntoineUnit

    !AMP 式の飽和蒸気圧式の係数
    a_ampA = ChemData_SvapPress_AMPA
    a_ampB = ChemData_SvapPress_AMPB
    a_ampC = ChemData_SvapPress_AMPC
    a_ampD = ChemData_SvapPress_AMPD
    a_ampE = ChemData_SvapPress_AMPE

    !分子量
    a_MolWt = ChemData_MolWt
    
    !NH4SH の反応熱の初期化
    !  NH4SH 1kg に対する反応熱にする.
    Name = 'NH4SH-s'
    id   = ChemData_OneSpcID( Name )  
    
    ReactHeatNH4SHPerMol  = GasRUniv * 10834.0d0
    ReactHeatNH4SH = GasRUniv * 10834.0d0 / MolWt( id )

    call MessageNotify( "M", "ChemCalc_Init", "ReactHeatNH4SH = %f", d=(/ReactHeatNH4SH/) )
    call MessageNotify( "M", "ChemCalc_Init", "NH4SH MolWt = %f", d=(/MolWt(id)/) )

  end subroutine ChemCalc_Init
Function :
CpPerMolRef :real(DP)
: 標準状態での単位モル当たりの比熱
ID :integer, intent(in)
: 化学種の ID

引数で与えられた化学種に対して, 標準状態での単位モル当たりの定圧比熱を計算

[Source]

  function CpPerMolRef(ID)
    !
    !引数で与えられた化学種に対して, 標準状態での単位モル当たりの定圧比熱を計算
    !

    !暗黙の型宣言禁止
    implicit none
    
    !入出力変数  
    real(DP)            :: CpPerMolRef  !標準状態での単位モル当たりの比熱
    integer, intent(in) :: ID           !化学種の ID

    
    !データベースから情報取得
    CpPerMolRef = ChemData_CpPerMolRef(ID)

  end function CpPerMolRef
Function :
CpRef :real(DP)
: 標準状態での単位質量当たりの比熱
ID :integer, intent(in)
: 化学種の ID

引数で与えられた化学種に対して, 標準状態での単位質量当たりの定圧比熱を計算

[Source]

  function CpRef(ID)
    !
    !引数で与えられた化学種に対して, 標準状態での単位質量当たりの定圧比熱を計算
    !

    !暗黙の型宣言禁止
    implicit none
    
    !入出力変数  
    real(DP)            :: CpRef        !標準状態での単位質量当たりの比熱
    integer, intent(in) :: ID           !化学種の ID

    
    !データベースから情報取得
    CpRef = ChemData_CpRef(ID)

  end function CpRef
Function :
CvRef :real(DP)
: 標準状態での単位質量当たりの比熱
ID :integer, intent(in)
: 化学種の ID

引数で与えられた化学種に対して, 標準状態での単位質量当たりの定圧比熱を計算

[Source]

  function CvRef(ID)
    !
    !引数で与えられた化学種に対して, 標準状態での単位質量当たりの定圧比熱を計算
    !

    !暗黙の型宣言禁止
    implicit none
    
    !入出力変数  
    real(DP)            :: CvRef       !標準状態での単位質量当たりの比熱
    integer, intent(in) :: ID           !化学種の ID

    
    !データベースから情報取得
    CvRef = ChemData_CvRef(ID)

  end function CvRef
Function :
GasR :real(DP)
: 分子量
ID :integer, intent(in)
: 化学種の ID

引数で与えられた化学種に対して, 気体定数を計算

[Source]

  function GasR(ID)
    !
    !引数で与えられた化学種に対して, 気体定数を計算
    !

    !暗黙の型宣言禁止
    implicit none
    
    !入出力変数  
    real(DP)            :: GasR          !分子量
    integer, intent(in) :: ID      !化学種の ID
    
    
    !データベースから情報取得
    GasR = ChemData_GasR(ID)

  end function GasR
Function :
LatentHeatPerMol :real(DP)
: 潜熱
ID :integer, intent(in)
: 化学種名
Temp :real(DP),intent(in)
: 温度

引数で与えられた化学種と温度に対して, 潜熱を計算

[Source]

  function LatentHeatPerMol(ID, Temp)
    !
    !引数で与えられた化学種と温度に対して, 潜熱を計算
    !

    !暗黙の型宣言禁止
    implicit none

    !入出力変数
    real(DP)            :: LatentHeatPerMol   !潜熱
    real(DP),intent(in) :: Temp            !温度
    integer, intent(in) :: ID                 !化学種名
    
    !内部変数
    real(DP)            :: DLogSvapPressDTemp
    real(DP),parameter  :: GasRUniv = 8.314d0  !普遍気体定数
    

    !飽和蒸気圧の温度微分
    DLogSvapPressDTemp = - a_ampA(ID) / (Temp ** 2.0d0) + a_ampC(ID) / Temp + a_ampD(ID) + a_ampE(ID) * 2.0d0 * Temp
        
    !潜熱の計算
    LatentHeatPerMol = DLogSvapPressDTemp * GasRUniv * (Temp ** 2.0d0) 

  end function LatentHeatPerMol
Function :
MolWt :real(DP)
: 分子量
ID :integer, intent(in)
: 化学種の ID

引数で与えられた化学種に対して, 分子量を計算

[Source]

  function MolWt(ID)
    !
    !引数で与えられた化学種に対して, 分子量を計算
    !

    !暗黙の型宣言禁止
    implicit none
    
    !入出力変数  
    real(DP)            :: MolWt         !分子量
    integer, intent(in) :: ID      !化学種の ID

    
    !データベースから情報取得
    MolWt = ChemData_MolWt(ID)

  end function MolWt
ReactHeatNH4SH
Variable :
ReactHeatNH4SH :real(DP)
: NH4SH 生成反応熱 [J/K kg]
ReactHeatNH4SHPerMol
Variable :
ReactHeatNH4SHPerMol :real(DP)
: NH4SH 生成反応熱 [J/K mol]
Function :
SvapPress :real(DP)
: 飽和蒸気圧
ID :integer, intent(in)
: 化学種の ID
Temp :real(DP),intent(in)
: 温度 [K]

引数で与えられた化学種と温度に対して, 飽和蒸気圧を計算. AMP 式を利用

[Source]

  function SvapPress(ID, Temp)
    !
    !引数で与えられた化学種と温度に対して, 飽和蒸気圧を計算. AMP 式を利用
    !

    !暗黙の型宣言禁止
    implicit none
    
    !入出力変数  
    real(DP)            :: SvapPress   !飽和蒸気圧
    real(DP),intent(in) :: Temp        !温度 [K]
    integer, intent(in) :: ID          !化学種の ID

    !内部変数
    real(DP)            :: LogSvapPress
    
    !飽和蒸気圧の対数を計算
    !対数が大きくなりすぎないようにする. 
    ! Fujitsu Fortran Compiler では 700 より大きい数の exp を取ると警告が出る.
    LogSvapPress = min( ( a_ampA(ID) / Temp + a_ampB(ID) + a_ampC(ID) * dlog( Temp ) + a_ampD(ID) * Temp + a_ampE(ID) * ( temp ** 2 ) + dlog(1.0d-1) ), 7.0d2 )

    !飽和蒸気圧を計算
    SvapPress =  dexp( LogSvapPress )
    
  end function SvapPress
Function :
xyz_DSvapPressDTemp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) :real(DP)
: 飽和蒸気圧の温度微分 [Pa/K]
ID :integer, intent(in)
: 化学種名
xyz_Temp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) :real(DP),intent(in)
: 温度 [K]

引数で与えられた化学種と温度に対して, 飽和蒸気圧の温度微分を計算

[Source]

  function xyz_DSvapPressDTemp(ID, xyz_Temp)
    !
    !引数で与えられた化学種と温度に対して, 飽和蒸気圧の温度微分を計算
    !
    
    !暗黙の型宣言禁止
    implicit none
  
    !入出力変数    
    real(DP)            :: xyz_DSvapPressDTemp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax)
                                                         !飽和蒸気圧の温度微分 [Pa/K]
    real(DP),intent(in) :: xyz_Temp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) 
                                                         !温度 [K]
    integer, intent(in) :: ID                            !化学種名

    !内部変数
    real(DP)            :: xyz_LogSvapPress(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax)
    real(DP)            :: xyz_DLogSvapPressDTemp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax)
    

    !飽和蒸気圧の対数を計算
    !対数が大きくなりすぎないようにする. 
    ! Fujitsu Fortran Compiler では 700 より大きい数の exp を取ると警告が出る.
    xyz_LogSvapPress = min( ( a_ampA(ID) / xyz_Temp + a_ampB(ID) + a_ampC(ID) * dlog( xyz_Temp ) + a_ampD(ID) * xyz_Temp + a_ampE(ID) * ( xyz_temp ** 2 ) + dlog(1.0d-1) ), 7.0d2 )

    !飽和蒸気圧の温度微分
    xyz_DLogSvapPressDTemp = - a_ampA(ID) / (xyz_Temp ** 2.0d0) + a_ampC(ID) / xyz_Temp + a_ampD(ID) + a_ampE(ID) * 2.0d0 * xyz_Temp
    
    !飽和蒸気圧の温度微分
    xyz_DSvapPressDTemp = xyz_DLogSvapPressDTemp * dexp( xyz_LogSvapPress ) 
    
  end function xyz_DSvapPressDTemp
Function :
xyz_LatentHeat(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) :real(DP)
: 潜熱
ID :integer, intent(in)
: 化学種名
xyz_Temp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) :real(DP),intent(in)
: 温度

引数で与えられた化学種と温度に対して, 潜熱を計算

[Source]

  function xyz_LatentHeat(ID, xyz_Temp)
    !
    !引数で与えられた化学種と温度に対して, 潜熱を計算
    !

    !暗黙の型宣言禁止
    implicit none

    !入出力変数
    real(DP)            :: xyz_LatentHeat(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax)
                                                         !潜熱
    real(DP),intent(in) :: xyz_Temp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) 
                                                         !温度
    integer, intent(in) :: ID                            !化学種名
    
    !内部変数
    real(DP)            :: xyz_DLogSvapPressDTemp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax)
    real(DP),parameter  :: GasRUniv = 8.314d0  !普遍気体定数
    

    !飽和蒸気圧の温度微分
    xyz_DLogSvapPressDTemp = - a_ampA(ID) / (xyz_Temp ** 2.0d0) + a_ampC(ID) / xyz_Temp + a_ampD(ID) + a_ampE(ID) * 2.0d0 * xyz_Temp
        
    !潜熱の計算
    xyz_LatentHeat = xyz_DLogSvapPressDTemp * GasRUniv * (xyz_Temp ** 2.0d0) / a_MolWt(ID)

  end function xyz_LatentHeat
Function :
xyz_SvapPress(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) :real(DP)
: 飽和蒸気圧
ID :integer, intent(in)
: 化学種の ID
xyz_Temp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) :real(DP),intent(in)
: 温度 [K]

引数で与えられた化学種と温度に対して, 飽和蒸気圧を計算. AMP 式を利用

[Source]

  function xyz_SvapPress(ID, xyz_Temp)
    !
    !引数で与えられた化学種と温度に対して, 飽和蒸気圧を計算. AMP 式を利用
    !

    !暗黙の型宣言禁止
    implicit none
    
    !入出力変数  
    real(DP)            :: xyz_SvapPress(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax)
                                                         !飽和蒸気圧
    real(DP),intent(in) :: xyz_Temp(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax) 
                                                         !温度 [K]
    integer, intent(in) :: ID                            !化学種の ID

    !内部変数
    real(DP)            :: xyz_LogSvapPress(DimXMin:DimXMax,DimYMin:DimYMax,DimZMin:DimZMax)
    
    !飽和蒸気圧の対数を計算
    !対数が大きくなりすぎないようにする. 
    ! Fujitsu Fortran Compiler では 700 より大きい数の exp を取ると警告が出る.
    xyz_LogSvapPress = min( ( a_ampA(ID) / xyz_Temp + a_ampB(ID) + a_ampC(ID) * dlog( xyz_Temp ) + a_ampD(ID) * xyz_Temp + a_ampE(ID) * ( xyz_temp ** 2 ) + dlog(1.0d-1) ), 7.0d2 )

    !飽和蒸気圧を計算
    xyz_SvapPress =  dexp( xyz_LogSvapPress )
    
  end function xyz_SvapPress