| Class | saturate_DennouAGCMExt |
| In: |
saturate/saturate_DennouAGCMExt.f90
|
Note that Japanese and English are described in parallel.
Dennou AGCM で用いられていた式を用いて飽和蒸気圧を計算します. この式は, 元々 沼口 (1992) に依っています.
Saturation specific humidity and temperature derivative of it are calculated with a formula used in Dennou AGCM. This formula is originaly obtained from Numaguti (1992).
| CalcQVapSat : | 飽和比湿の計算 |
| CalcDQVapSatDTemp : | 飽和比湿の温度微分の計算 |
| ———— : | ———— |
| CalcQVapSat : | Calculate saturation specific humidity |
| CalcDQVapSatDTemp : | Calculate temperature derivative of saturation specific humidity |
| Subroutine : |
saturate_DennouAGCM モジュールの初期化を行います.
"saturate_DennouAGCM" module is initialized.
subroutine SaturateInit
!
! saturate_DennouAGCM モジュールの初期化を行います.
!
!
! "saturate_DennouAGCM" module is initialized.
!
! モジュール引用 ; 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
! 文字列操作
! Character handling
!
use dc_string, only: StoA
! ヒストリデータ出力
! History data output
!
use gtool_historyauto, only: HistoryAutoAddVariable
! 宣言文 ; 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 /saturate_DennouAGCM_nml/ &
!!$ & FlagSublimation
!
! デフォルト値については初期化手続 "saturate_DennouAGCM#SaturateInit"
! のソースコードを参照のこと.
!
! Refer to source codes in the initialization procedure
! "saturate_DennouAGCM#SaturateInit" for the default values.
!
! 実行文 ; Executable statement
!
if ( saturate_DennouAGCMExt_inited ) return
! デフォルト値の設定
! Default values settings
!
! NAMELIST の読み込み
! NAMELIST is input
!
!!$ if ( trim(namelist_filename) /= '' ) then
!!$ call FileOpen( unit_nml, & ! (out)
!!$ & namelist_filename, mode = 'r' ) ! (in)
!!$
!!$ rewind( unit_nml )
!!$ read( unit_nml, & ! (in)
!!$ & nml = saturate_DennouAGCM_nml, & ! (out)
!!$ & iostat = iostat_nml ) ! (out)
!!$ close( unit_nml )
!!$
!!$ call NmlutilMsg( iostat_nml, module_name ) ! (in)
!!$ end if
! Initialization of modules used in this moduled
!
call SaturateInitOnLiq
call SaturateInitOnSol
! 印字 ; Print
!
call MessageNotify( 'M', module_name, '----- Initialization Messages -----' )
!!$ call MessageNotify( 'M', module_name, ' FlagSublimation = %b', l = (/ FlagSublimation /) )
call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) )
saturate_DennouAGCMExt_inited = .true.
end subroutine SaturateInit
| Variable : | |||
| saturate_DennouAGCMExt_inited = .false. : | logical, save, public
|
| Function : | |||
| xyz_DQVapSatDTemp(size(xyz_Temp,1), size(xyz_Temp,2), size(xyz_Temp,3)) : | real(DP)
| ||
| xyz_Temp(:,:,:) : | real(DP), intent(in)
| ||
| xyz_QVapSat(:,:,:) : | real(DP), intent(in)
|
温度 Temp と飽和比湿 QVapSat を用い, 飽和比湿の温度微分 DQVapSatDTemp を求めます.
Calculate temperature derivative of saturation specific humidity DQVapSatDTemp using temperature Temp and saturation specific humidity QVapSat.
function xyz_CalcDQVapSatDTemp( xyz_Temp, xyz_QVapSat ) result( xyz_DQVapSatDTemp )
!
! 温度 *Temp* と飽和比湿 *QVapSat* を用い,
! 飽和比湿の温度微分 *DQVapSatDTemp* を求めます.
!
! Calculate temperature derivative of saturation specific humidity
! *DQVapSatDTemp* using
! temperature *Temp* and saturation specific humidity *QVapSat*.
!
! モジュール引用 ; USE statements
!
! 宣言文 ; Declaration statements
!
implicit none
real(DP), intent(in):: xyz_Temp (:,:,:)
! $ T $ . 温度. Temperature
real(DP), intent(in):: xyz_QVapSat(:,:,:)
! $ q^{*} $ . 飽和比湿. Saturation specific humidity
real(DP):: xyz_DQVapSatDTemp(size(xyz_Temp,1), size(xyz_Temp,2), size(xyz_Temp,3))
! $ \DP{q^{*}}{T} $ . 飽和比湿の温度微分.
! Temperature derivative of saturation specific humidity.
! 作業変数
! Work variables
!
real(DP):: xyz_DQVapSatDTempOnLiq(size(xyz_Temp,1), size(xyz_Temp,2), size(xyz_Temp,3))
real(DP):: xyz_DQVapSatDTempOnSol(size(xyz_Temp,1), size(xyz_Temp,2), size(xyz_Temp,3))
real(DP):: xyz_WatFrac (size(xyz_Temp,1), size(xyz_Temp,2), size(xyz_Temp,3))
! 実行文 ; Executable statement
!
if ( .not. saturate_DennouAGCMExt_inited ) then
call MessageNotify( 'E', module_name, 'This module has not been initialized.' )
end if
xyz_DQVapSatDTempOnLiq = xyz_CalcDQVapSatDTempOnLiq( xyz_Temp, xyz_QVapSat )
xyz_DQVapSatDTempOnSol = xyz_CalcDQVapSatDTempOnSol( xyz_Temp, xyz_QVapSat )
call CloudUtilsWatFraction( xyz_Temp, xyz_WatFrac )
xyz_DQVapSatDTemp = xyz_WatFrac * xyz_DQVapSatDTempOnLiq + ( 1.0_DP - xyz_WatFrac ) * xyz_DQVapSatDTempOnSol
end function xyz_CalcDQVapSatDTemp
| Function : | |||
| xyz_QVapSat(1:size(xyz_Temp,1), 1:size(xyz_Temp,2), 1:size(xyz_Temp,3)) : | real(DP)
| ||
| xyz_Temp(:,:,:) : | real(DP), intent(in)
| ||
| xyz_Press(:,:,:) : | real(DP), intent(in)
|
温度 Temp と気圧 Press を用い, 飽和比湿 QVapSat を求めます.
Calculate saturation specific humidity QVapSat using temperature Temp and air pressure Press.
function xyz_CalcQVapSat( xyz_Temp, xyz_Press ) result( xyz_QVapSat )
!
! 温度 *Temp* と気圧 *Press* を用い,
! 飽和比湿 *QVapSat* を求めます.
!
! Calculate saturation specific humidity *QVapSat* using
! temperature *Temp* and air pressure *Press*.
!
! モジュール引用 ; USE statements
!
! 宣言文 ; Declaration statements
!
implicit none
real(DP), intent(in):: xyz_Temp (:,:,:)
! $ T $ . 温度. Temperature
real(DP), intent(in):: xyz_Press(:,:,:)
! $ p $ . 気圧. Air pressure
real(DP):: xyz_QVapSat(1:size(xyz_Temp,1), 1:size(xyz_Temp,2), 1:size(xyz_Temp,3))
! $ q^{*} $ . 飽和比湿. Saturation specific humidity
! 作業変数
! Work variables
!
real(DP):: xyz_QVapSatOnLiq(1:size(xyz_Temp,1), 1:size(xyz_Temp,2), 1:size(xyz_Temp,3))
real(DP):: xyz_QVapSatOnSol(1:size(xyz_Temp,1), 1:size(xyz_Temp,2), 1:size(xyz_Temp,3))
real(DP):: xyz_WatFrac (1:size(xyz_Temp,1), 1:size(xyz_Temp,2), 1:size(xyz_Temp,3))
! 実行文 ; Executable statement
!
! 初期化確認
! Initialization check
!
if ( .not. saturate_DennouAGCMExt_inited ) then
call MessageNotify( 'E', module_name, 'This module has not been initialized.' )
end if
xyz_QVapSatOnLiq = xyz_CalcQVapSatOnLiq( xyz_Temp, xyz_Press )
xyz_QVapSatOnSol = xyz_CalcQVapSatOnSol( xyz_Temp, xyz_Press )
call CloudUtilsWatFraction( xyz_Temp, xyz_WatFrac )
xyz_QVapSat = xyz_WatFrac * xyz_QVapSatOnLiq + ( 1.0_DP - xyz_WatFrac ) * xyz_QVapSatOnSol
end function xyz_CalcQVapSat
| Constant : | |||
| module_name = ‘saturate_DennouAGCMExt‘ : | character(*), parameter
|
| Constant : | |||
| version = ’$Name: $’ // ’$Id: saturate_DennouAGCMExt.f90,v 1.1 2014/06/29 07:26:29 yot Exp $’ : | character(*), parameter
|