Class: DExnerDt

モデルの放射過程を計算するためのパッケージ型モジュール具体的には以下の項を計算するための関数を格納する.

 * 一様冷却

Methods

xy_DExnerDt_xy_xyf xyz_DExnerDt_xyz xyz_DExnerDt_xyz_xyzf xyz_DExnerDt_xyzf

Included Modules

dc_types dc_iounit dc_message gridSet basicset constants composition

Public Instance methods

xy_DExnerDt_xy_xyf( xy_DPTempDt, xyf_DQMixDt, kz ) result(xy_DExnerDt_xy_xyf)

*Function* :
`xy_DExnerDt_xy_xyf(imin:imax,jmin:jmax)` :	`real(DP)`
`xy_DPTempDt(imin:imax,jmin:jmax)` :	`real(DP), intent(in)`
`xyf_DQMixDt(imin:imax,jmin:jmax,ncmax)` :	`real(DP), intent(in)`
`kz` :	`integer, intent(in)`

[Source]

  function xy_DExnerDt_xy_xyf( xy_DPTempDt, xyf_DQMixDt, kz )

    !暗黙の型宣言禁止
    implicit none

    !入力変数
    real(DP), intent(in) :: xy_DPTempDt(imin:imax,jmin:jmax)
    real(DP), intent(in) :: xyf_DQMixDt(imin:imax,jmin:jmax,ncmax)
    integer,  intent(in) :: kz
    real(DP)             :: xy_DExnerDt_xy_xyf(imin:imax,jmin:jmax)
    real(DP)             :: xyf_DQMixDtPerMolWt(imin:imax,jmin:jmax,ncmax)
    integer              :: s

    do s = 1, ncmax
      xyf_DQMixDtPerMolWt(:,:,s) = xyf_DQMixDt(:,:,s) / MolWtWet(s)
    end do
    
    xy_DExnerDt_xy_xyf = (xyz_VelSoundBZ(:,:,kz) ** 2.0d0) / (CpDry * xyz_VPTempBZ(:,:,kz)) * ( xy_DPTempDt / xyz_PTempBZ(:,:,kz) + sum(xyf_DQMixDtPerMolWt(:,:,1:GasNum),3) / ( 1.0d0 / MolWtDry + xyz_QMixBZPerMolWt(:,:,kz) ) - sum(xyf_DQMixDt,3) / ( 1.0d0 + xyz_QMixBZ(:,:,kz) ) )
    
  end function xy_DExnerDt_xy_xyf

xyz_DExnerDt_xyz( xyz_DPTempDt ) result(xyz_DExnerDt_xyz)

*Function* :
`xyz_DExnerDt_xyz(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP)`
`xyz_DPTempDt(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP), intent(in)`

[Source]

  function xyz_DExnerDt_xyz( xyz_DPTempDt )

    !暗黙の型宣言禁止
    implicit none
    
    !入力変数
    real(DP), intent(in) :: xyz_DPTempDt(imin:imax,jmin:jmax,kmin:kmax)
    real(DP)             :: xyz_DExnerDt_xyz(imin:imax,jmin:jmax,kmin:kmax)
    
    xyz_DExnerDt_xyz = (xyz_VelSoundBZ ** 2.0d0) / (CpDry * xyz_VPTempBZ) * ( xyz_DPTempDt / xyz_PTempBZ )
    
  end function xyz_DExnerDt_xyz

xyz_DExnerDt_xyz_xyzf( xyz_DPTempDt, xyzf_DQMixDt ) result(xyz_DExnerDt_xyz_xyzf)

*Function* :
`xyz_DExnerDt_xyz_xyzf(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP)`
`xyz_DPTempDt(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP), intent(in)`
`xyzf_DQMixDt(imin:imax,jmin:jmax,kmin:kmax,ncmax)` :	`real(DP), intent(in)`

[Source]

  function xyz_DExnerDt_xyz_xyzf( xyz_DPTempDt, xyzf_DQMixDt )

    !暗黙の型宣言禁止
    implicit none

    !入力変数
    real(DP), intent(in) :: xyz_DPTempDt(imin:imax,jmin:jmax,kmin:kmax)
    real(DP), intent(in) :: xyzf_DQMixDt(imin:imax,jmin:jmax,kmin:kmax,ncmax)
    real(DP)             :: xyz_DExnerDt_xyz_xyzf(imin:imax,jmin:jmax,kmin:kmax)
    real(DP)             :: xyzf_DQMixDtPerMolWt(imin:imax,jmin:jmax,kmin:kmax,ncmax)
    integer              :: s
    
    do s = 1, ncmax
      xyzf_DQMixDtPerMolWt(:,:,:,s) = xyzf_DQMixDt(:,:,:,s) / MolWtWet(s)
    end do
    xyz_DExnerDt_xyz_xyzf = (xyz_VelSoundBZ ** 2.0d0) / (CpDry * xyz_VPTempBZ) * ( xyz_DPTempDt / xyz_PTempBZ + sum(xyzf_DQMixDtPerMolWt(:,:,:,1:GasNum),4) / ( 1.0d0 / MolWtDry + xyz_QMixBZPerMolWt ) - sum(xyzf_DQMixDt,4) / ( 1.0d0 + xyz_QMixBZ ) )

  end function xyz_DExnerDt_xyz_xyzf

xyz_DExnerDt_xyzf( xyzf_DQMixDt ) result(xyz_DExnerDt_xyzf)

*Function* :
`xyz_DExnerDt_xyzf(imin:imax,jmin:jmax,kmin:kmax)` :	`real(DP)`
`xyzf_DQMixDt(imin:imax,jmin:jmax,kmin:kmax,ncmax)` :	`real(DP), intent(in)`

[Source]

  function xyz_DExnerDt_xyzf( xyzf_DQMixDt )

    !暗黙の型宣言禁止
    implicit none

    !入力変数
    real(DP), intent(in) :: xyzf_DQMixDt(imin:imax,jmin:jmax,kmin:kmax,ncmax)
    real(DP)             :: xyz_DExnerDt_xyzf(imin:imax,jmin:jmax,kmin:kmax)
    real(DP)             :: xyzf_DQMixDtPerMolWt(imin:imax,jmin:jmax,kmin:kmax,ncmax)
    integer              :: s
    
    do s = 1, ncmax
      xyzf_DQMixDtPerMolWt(:,:,:,s) = xyzf_DQMixDt(:,:,:,s) / MolWtWet(s)
    end do
    xyz_DExnerDt_xyzf = (xyz_VelSoundBZ ** 2.0d0) / (CpDry * xyz_VPTempBZ) * ( + sum(xyzf_DQMixDtPerMolWt(:,:,:,1:GasNum),4) / ( 1.0d0 / MolWtDry + xyz_QMixBZPerMolWt ) - sum(xyzf_DQMixDt,4) / ( 1.0d0 + xyz_QMixBZ ) )

  end function xyz_DExnerDt_xyzf