| Class | check_prog_vars |
| In: |
util/check_prog_vars.f90
|
Note that Japanese and English are described in parallel.
予報変数の値の確認
Check values of prognostic variables
| !$ ! SurfaceFlux : | 地表面フラックスの計算 |
| !$ ! SurfaceFluxOutput : | 地表面フラックスの出力 |
| !$ ! ———— : | ———— |
| !$ ! SurfaceFlux : | Calculate surface fluxes |
| !$ ! SurfaceFluxOutput : | Output surface fluxes |
| Subroutine : | |||
| xy_Ps(0:imax-1, 1:jmax) : | real(DP), intent(in)
| ||
| xyz_U(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(in)
| ||
| xyz_V(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(in)
| ||
| xyz_Temp(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(in)
| ||
| xyzf_QMix(0:imax-1, 1:jmax, 1:kmax, 1:ncmax) : | real(DP), intent(in)
|
予報変数の値の確認
Check values of prognostic variables
subroutine CheckProgVars( xy_Ps, xyz_U, xyz_V, xyz_Temp, xyzf_QMix )
!
! 予報変数の値の確認
!
! Check values of prognostic variables
!
! モジュール引用 ; USE statements
!
! 格子点設定
! Grid points settings
!
use gridset, only: imax, jmax, kmax ! 鉛直層数.
! Number of vertical level
! 組成に関わる配列の設定
! Settings of array for atmospheric composition
!
use composition, only: ncmax, IndexH2OVap
! メッセージ出力
! Message output
!
use dc_message, only: MessageNotify
! 時刻管理
! Time control
!
use timeset, only: TimeN, TimesetClockStart, TimesetClockStop
! デバッグ用ユーティリティ
! Utilities for debug
!
use dc_trace, only: DbgMessage, BeginSub, EndSub
! 宣言文 ; Declaration statements
!
implicit none
real(DP), intent(in):: xy_Ps (0:imax-1, 1:jmax)
! $ p_s (t-\Delta t) $ . 地表面気圧. Surface pressure (Pa)
real(DP), intent(in):: xyz_U (0:imax-1, 1:jmax, 1:kmax)
! $ u $ . 東西風速. Eastward wind
real(DP), intent(in):: xyz_V (0:imax-1, 1:jmax, 1:kmax)
! $ v $ . 南北風速. Northward wind
real(DP), intent(in):: xyz_Temp (0:imax-1, 1:jmax, 1:kmax)
! $ T $ . 温度 (整数レベル).
! Temperature (full level)
real(DP), intent(in):: xyzf_QMix(0:imax-1, 1:jmax, 1:kmax, 1:ncmax)
! $ q $ . 比湿. Specific humidity
! 作業変数
! Work variables
!
integer:: i ! 経度方向に回る DO ループ用作業変数
! Work variables for DO loop in longitude
integer:: j ! 緯度方向に回る DO ループ用作業変数
! Work variables for DO loop in latitude
integer:: k ! 鉛直方向に回る DO ループ用作業変数
! Work variables for DO loop in vertical direction
integer:: n ! 組成方向に回る DO ループ用作業変数
! Work variables for DO loop in dimension of constituents
real(DP) :: minval, maxval
integer :: mini, minj, mink
integer :: maxi, maxj, maxk
! 実行文 ; Executable statement
!
! 計算時間計測開始
! Start measurement of computation time
!
call TimesetClockStart( module_name )
! 初期化
! Initialization
!
if ( .not. check_prog_vars_inited ) call CheckProgVarsInit
i = 0
j = 1
k = 1
minval = xyz_U(i,j,k)
maxval = xyz_U(i,j,k)
mini = i
minj = j
mink = k
maxi = i
maxj = j
maxk = k
do k = 1, kmax
do j = 1, jmax
do i = 0, imax-1
if ( xyz_U (i,j,k) < minval ) then
minval = xyz_U (i,j,k)
mini = i
minj = j
mink = k
end if
if ( maxval < xyz_U (i,j,k) ) then
maxval = xyz_U (i,j,k)
maxi = i
maxj = j
maxk = k
end if
end do
end do
end do
if ( minval < VelMin ) then
call MessageNotify( 'M', module_name, 'Extreme value : U(%d,%d,%d) = %f.', i = (/ mini, minj, mink /), d = (/ minval /) )
end if
if ( VelMax < maxval ) then
call MessageNotify( 'M', module_name, 'Extreme value : U(%d,%d,%d) = %f.', i = (/ maxi, maxj, maxk /), d = (/ maxval /) )
end if
i = 0
j = 1
k = 1
minval = xyz_V(i,j,k)
maxval = xyz_V(i,j,k)
mini = i
minj = j
mink = k
maxi = i
maxj = j
maxk = k
do k = 1, kmax
do j = 1, jmax
do i = 0, imax-1
if ( xyz_V (i,j,k) < minval ) then
minval = xyz_V (i,j,k)
mini = i
minj = j
mink = k
end if
if ( maxval < xyz_V (i,j,k) ) then
maxval = xyz_V (i,j,k)
maxi = i
maxj = j
maxk = k
end if
end do
end do
end do
if ( minval < VelMin ) then
call MessageNotify( 'M', module_name, 'Extreme value : V(%d,%d,%d) = %f.', i = (/ mini, minj, mink /), d = (/ minval /) )
end if
if ( VelMax < maxval ) then
call MessageNotify( 'M', module_name, 'Extreme value : V(%d,%d,%d) = %f.', i = (/ maxi, maxj, maxk /), d = (/ maxval /) )
end if
i = 0
j = 1
k = 1
minval = xyz_Temp(i,j,k)
maxval = xyz_Temp(i,j,k)
mini = i
minj = j
mink = k
maxi = i
maxj = j
maxk = k
do k = 1, kmax
do j = 1, jmax
do i = 0, imax-1
if ( xyz_Temp(i,j,k) < minval ) then
minval = xyz_Temp(i,j,k)
mini = i
minj = j
mink = k
end if
if ( maxval < xyz_Temp(i,j,k) ) then
maxval = xyz_Temp(i,j,k)
maxi = i
maxj = j
maxk = k
end if
end do
end do
end do
if ( minval < TempMin ) then
call MessageNotify( 'M', module_name, 'Extreme value : Temp(%d,%d,%d) = %f.', i = (/ mini, minj, mink /), d = (/ minval /) )
end if
if ( TempMax < maxval ) then
call MessageNotify( 'M', module_name, 'Extreme value : Temp(%d,%d,%d) = %f.', i = (/ maxi, maxj, maxk /), d = (/ maxval /) )
end if
! ヒストリデータ出力
! History data output
!
! 計算時間計測一時停止
! Pause measurement of computation time
!
call TimesetClockStop( module_name )
end subroutine CheckProgVars
| Variable : | |||
| check_prog_vars_inited = .false. : | logical, save, public
|
| Subroutine : |
check_prog_vars モジュールの初期化を行います. NAMELIST#check_prog_vars_nml の読み込みはこの手続きで行われます.
"check_prog_vars" module is initialized. "NAMELIST#check_prog_vars_nml" is loaded in this procedure.
This procedure input/output NAMELIST#check_prog_vars_nml .
subroutine CheckProgVarsInit
!
! check_prog_vars モジュールの初期化を行います.
! NAMELIST#check_prog_vars_nml の読み込みはこの手続きで行われます.
!
! "check_prog_vars" module is initialized.
! "NAMELIST#check_prog_vars_nml" is loaded in this procedure.
!
! モジュール引用 ; USE statements
!
! NAMELIST ファイル入力に関するユーティリティ
! Utilities for NAMELIST file input
!
use namelist_util, only: namelist_filename, NmlutilMsg, NmlutilAryValid
! ファイル入出力補助
! 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
! メッセージ出力
! 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 /check_prog_vars_nml/ VelMin, VelMax, TempMin, TempMax
!
! デフォルト値については初期化手続 "check_prog_vars#CheckProgVarsInit"
! のソースコードを参照のこと.
!
! Refer to source codes in the initialization procedure
! "check_prog_vars#CheckProgVarsInit" for the default values.
!
! 実行文 ; Executable statement
!
if ( check_prog_vars_inited ) return
!!$ call InitCheck
! デフォルト値の設定
! Default values settings
!
VelMin = -300.0_DP
VelMax = 300.0_DP
TempMin = 50.0_DP
TempMax = 350.0_DP
! NAMELIST の読み込み
! NAMELIST is input
!
if ( trim(namelist_filename) /= '' ) then
call FileOpen( unit_nml, namelist_filename, mode = 'r' ) ! (in)
rewind( unit_nml )
read( unit_nml, nml = check_prog_vars_nml, iostat = iostat_nml ) ! (out)
close( unit_nml )
call NmlutilMsg( iostat_nml, module_name ) ! (in)
end if
! 印字 ; Print
!
call MessageNotify( 'M', module_name, '----- Initialization Messages -----' )
call MessageNotify( 'M', module_name, ' VelMin = %f', d = (/ VelMin /) )
call MessageNotify( 'M', module_name, ' VelMax = %f', d = (/ VelMax /) )
call MessageNotify( 'M', module_name, ' TempMin = %f', d = (/ TempMin /) )
call MessageNotify( 'M', module_name, ' TempMax = %f', d = (/ TempMax /) )
call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) )
check_prog_vars_inited = .true.
end subroutine CheckProgVarsInit
| Constant : | |||
| module_name = ‘check_prog_vars‘ : | character(*), parameter
|
| Constant : | |||
| version = ’$Name: dcpam5-20100224 $’ // ’$Id: check_prog_vars.f90,v 1.2 2010-01-21 15:41:31 yot Exp $’ : | character(*), parameter
|