| Class | held_suarez_1994 |
| In: |
held_suarez_1994/held_suarez_1994.f90
|
Note that Japanese and English are described in parallel.
Held and Suarez (1994) で提案される乾燥大気 GCM ベンチマーク用の 強制と散逸を計算します. 与える強制と散逸として, 温度場の帯状対称場への簡単なニュートン冷却と, 境界層摩擦を表現する下層風のレイリー摩擦を用います. 詳細を以下に記します.
Forcing and dissipation for dry air GCM benchmark suggested by Held and Suarez (1994) are caluclate. We use simple Newtonian relaxation of the temperature field to a zonally symmetric state and Rayleigh damping of low-level winds to represent boundary-layer friction. Their specifications are detailed as follows.
\[
\left( \DP{\Dvect{v}}{t} \right)_{\mathrm{HS94}} =
- k_v (\sigma) \Dvect{v}, \] \[
\left( \DP{T}{t} \right)_{\mathrm{HS94}} =
- k_T (\phi, \sigma) [T - T_{eq} (\phi,p)], \] \[
T_{eq} = \mathrm{max}
\left\{
200 \mathrm{K},
\left[
315 \mathrm{K} - (\Delta T)_y \sin^2\phi
- (\Delta \theta)_z
\log \left(\frac{p}{p_0}\right) \cos^2\phi
\right] \left(\frac{p}{p_0}\right)^\kappa
\right\}, \] \[
k_T = k_a + (k_s - k_a)
\mathrm{max}
\left(0, \frac{\sigma - \sigma_b}{1 - \sigma_b}\right) \cos^4\phi,
\] \[
k_v = k_f
\mathrm{max}
\left(0, \frac{\sigma - \sigma_b}{1 - \sigma_b}\right),
\] \[
\sigma_b = 0.7, \qquad
k_f = 1 \mathrm{day}^{-1}, \qquad
k_a = \Dinv{40} \mathrm{day}^{-1}, \qquad
k_s = \Dinv{4} \mathrm{day}^{-1}, \] \[
(\Delta T)_y = 60 \mathrm{K}, \qquad
(\Delta \theta)_z = 10 \mathrm{K}, \qquad
p_0 = 1000 \mathrm{hPa}, \qquad
\kappa = \frac{R}{c_p}.
\]
Forcing では, 与えられた速度や温度 ( $ t+Delta t$ を想定) に対して以下のように強制と散逸を適用します.
By Forcing, forcing and dissipation are applied to given wind and temperature ($ t+Delta t$ is expected) as follows.
\[
\hat{\Dvect{v}}^{t+\Delta t} =
\Dvect{v}^{t+\Delta t}
+ 2 \Delta t \left( \DP{\Dvect{v}}{t} \right)_{\mathrm{HS94}} \] \[
\hat{T}^{t+\Delta t} =
T^{t+\Delta t}
+ 2 \Delta t \left( \DP{T}{t} \right)_{\mathrm{HS94}}
\]
| Hs94Forcing : | 強制と散逸の計算 |
| Hs94Finalize : | 終了処理 (モジュール内部の変数の割り付け解除) |
| ———— : | ———— |
| Hs94Forcing : | Calculate forcing and dissipation |
| Hs94Finalize : | Termination (deallocate variables in this module) |
| Subroutine : |
モジュール内部の変数の割り付け解除を行います.
Deallocate variables in this module.
subroutine Hs94Finalize
!
! モジュール内部の変数の割り付け解除を行います.
!
! Deallocate variables in this module.
!
! 宣言文 ; Declaration statements
!
implicit none
! 実行文 ; Executable statement
!
if ( .not. held_suarez_1994_inited ) return
! 割り付け解除
! Deallocation
!
if ( allocated( z_kv ) ) deallocate( z_kv )
if ( allocated( yz_kt ) ) deallocate( yz_kt )
held_suarez_1994_inited = .false.
end subroutine Hs94Finalize
| Subroutine : | |||
| 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)
| ||
| xy_Ps(0:imax-1, 1:jmax) : | real(DP), intent(in)
| ||
| xyz_DUDt(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(out)
| ||
| xyz_DVDt(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(out)
| ||
| xyz_DTempDt(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(out)
|
引数として与えられた東西風速 xyz_U, 南北風速 xyz_V, 温度 xyz_Temp から, 温度場の帯状対称場への簡単なニュートン冷却と 境界層摩擦を表現する下層風のレイリー摩擦による 風速と温度の変化率を求め, xyz_DUDt, xyz_DVDt, xyz_DTempDt に返します.
Tendencies by simple Newtonian relaxation of the temperature field to a zonally symmetric state and Rayleigh damping of low-level winds to represent boundary-layer friction are calculated from eastward wind "xyz_U", northward wind "xyz_V", temperature "xyz_Temp". And the tencencies are returned as "xyz_DUDt", "xyz_DVDt", "xyz_DTempDt".
subroutine Hs94Forcing( xyz_U, xyz_V, xyz_Temp, xy_Ps, xyz_DUDt, xyz_DVDt, xyz_DTempDt )
!
! 引数として与えられた東西風速 xyz_U, 南北風速 xyz_V,
! 温度 xyz_Temp から,
! 温度場の帯状対称場への簡単なニュートン冷却と
! 境界層摩擦を表現する下層風のレイリー摩擦による
! 風速と温度の変化率を求め,
! xyz_DUDt, xyz_DVDt, xyz_DTempDt に返します.
!
! Tendencies by simple Newtonian relaxation of the temperature field to a
! zonally symmetric state and Rayleigh damping of low-level winds to
! represent boundary-layer friction are calculated
! from eastward wind "xyz_U", northward wind "xyz_V",
! temperature "xyz_Temp".
! And the tencencies are returned as
! "xyz_DUDt", "xyz_DVDt", "xyz_DTempDt".
!
!
! モジュール引用 ; USE statements
!
! 座標データ設定
! Axes data settings
!
use axesset, only: y_Lat, z_Sigma ! $ \sigma $ レベル (整数).
! Full $ \sigma $ level
! 時刻管理
! Time control
!
use timeset, only: DelTime, TimeN, TimesetClockStart, TimesetClockStop
! ヒストリデータ出力
! History data output
!
use gtool_historyauto, only: HistoryAutoPut
! 宣言文 ; Declaration statements
!
implicit none
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
real(DP), intent(in):: xy_Ps (0:imax-1, 1:jmax)
! $ p_s $ . 地表面気圧.
! Surface pressure
real(DP), intent(out):: xyz_DUDt (0:imax-1, 1:jmax, 1:kmax)
! $ \DP{u}{t} $ . 東西風速変化.
! Eastward wind tendency
real(DP), intent(out):: xyz_DVDt (0:imax-1, 1:jmax, 1:kmax)
! $ \DP{v}{t} $ . 南北風速変化.
! Northward wind tendency
real(DP), intent(out):: xyz_DTempDt (0:imax-1, 1:jmax, 1:kmax)
! $ \DP{T}{t} $ . 温度変化.
! Temperature tendency
! 作業変数
! Work variables
!
real(DP):: xyz_TempEQ (0:imax-1, 1:jmax, 1:kmax)
! $ T_{eq} $ . 平衡温度.
! Equilibrium temperature
real(DP):: xyz_Press (0:imax-1, 1:jmax, 1:kmax)
! $ T $ . 圧力.
! Pressure
integer:: j ! 緯度方向に回る DO ループ用作業変数
! Work variables for DO loop in latitude
integer:: k ! 鉛直方向に回る DO ループ用作業変数
! Work variables for DO loop in vertical direction
! 実行文 ; Executable statement
!
! 計算時間計測開始
! Start measurement of computation time
!
call TimesetClockStart( module_name )
! 初期化
! Initialization
!
if ( .not. held_suarez_1994_inited ) call Hs94Init
! 東西風速 $ u $ と南北風速 $ v $ へレイリー摩擦を適用
! Apply Rayleigh damping to eastward wind $ u $ and northward wind $ v $
!
do k = 1, kmax
xyz_DUDt (:,:,k) = - z_kv (k) * xyz_U (:,:,k)
xyz_DVDt (:,:,k) = - z_kv (k) * xyz_V (:,:,k)
end do
! 温度 $ T $ へニュートン冷却を適用
! Apply Newtonian relaxation to temperature $ T $
!
do k = 1, kmax
xyz_Press(:,:,k) = z_Sigma(k) * xy_Ps
enddo
do j = 1, jmax
xyz_TempEQ(:,j,:) = max( 200.0_DP, ( 315.0_DP - DelTempY * sin( y_Lat(j) ) ** 2 - DelPotTempZ * log( xyz_Press(:,j,:) / P0 ) * cos( y_Lat(j) ) ** 2 ) * ( xyz_Press(:,j,:) / P0 ) ** Kappa )
end do
do k = 1, kmax
do j = 1, jmax
xyz_DTempDt (:,j,k) = - yz_kt (j,k) * ( xyz_Temp (:,j,k) - xyz_TempEQ (:,j,k) )
end do
end do
! ヒストリデータ出力
! History data output
!
call HistoryAutoPut( TimeN, 'DUDtHs94', xyz_DUDt )
call HistoryAutoPut( TimeN, 'DVDtHs94', xyz_DVDt )
call HistoryAutoPut( TimeN, 'DTempDtHs94', xyz_DTempDt )
call HistoryAutoPut( TimeN, 'TempEQHs94', xyz_TempEQ )
! 計算時間計測一時停止
! Pause measurement of computation time
!
call TimesetClockStop( module_name )
end subroutine Hs94Forcing
| Variable : | |||
| held_suarez_1994_inited = .false. : | logical, save, public
|
| Subroutine : |
held_suarez_1994 モジュールの初期化を行います. NAMELIST#held_suarez_1994_nml の読み込みはこの手続きで行われます.
"held_suarez_1994" module is initialized. "NAMELIST#held_suarez_1994_nml" is loaded in this procedure.
subroutine Hs94Init
!
! held_suarez_1994 モジュールの初期化を行います.
! NAMELIST#held_suarez_1994_nml の読み込みはこの手続きで行われます.
!
! "held_suarez_1994" module is initialized.
! "NAMELIST#held_suarez_1994_nml" is loaded in this procedure.
!
! モジュール引用 ; USE statements
!
! 物理定数設定
! Physical constants settings
!
use constants, only: GasRDry, CpDry
! $ C_p $ [J kg-1 K-1].
! 乾燥大気の定圧比熱.
! Specific heat of air at constant pressure
! 座標データ設定
! Axes data settings
!
use axesset, only: y_Lat, z_Sigma ! $ \sigma $ レベル (整数).
! Full $ \sigma $ level
! 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
! ヒストリデータ出力
! History data output
!
use gtool_historyauto, only: HistoryAutoAddVariable
! 宣言文 ; Declaration statements
!
implicit none
real(DP):: day_seconds
! 1 日の秒数.
! Seconds in day.
real(DP):: SigmaB ! $ \sigma_b $ .
real(DP):: kf ! $ k_f $ .
real(DP):: ka ! $ k_a $ .
real(DP):: ks ! $ k_s $ .
integer:: j ! 緯度方向に回る DO ループ用作業変数
! Work variables for DO loop in latitude
integer:: k ! 鉛直方向に回る DO ループ用作業変数
! Work variables for DO loop in vertical direction
!!$ integer:: unit_nml ! NAMELIST ファイルオープン用装置番号.
!!$ ! Unit number for NAMELIST file open
!!$ integer:: iostat_nml ! NAMELIST 読み込み時の IOSTAT.
!!$ ! IOSTAT of NAMELIST read
! NAMELIST 変数群
! NAMELIST group name
!
!!$ namelist /held_suarez_1994_nml/
!!$ !
!!$ ! デフォルト値については初期化手続 "held_suarez_1994#Hs94Init"
!!$ ! のソースコードを参照のこと.
!!$ !
!!$ ! Refer to source codes in the initialization procedure
!!$ ! "held_suarez_1994#Hs94Init" for the default values.
!!$ !
! 実行文 ; Executable statement
!
if ( held_suarez_1994_inited ) return
call InitCheck
! デフォルト値の設定
! 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 = held_suarez_1994_nml, & ! (out)
!!$ & iostat = iostat_nml ) ! (out)
!!$ close( unit_nml )
!!$
!!$ call NmlutilMsg( iostat_nml, module_name ) ! (in)
!!$!$ if ( iostat_nml == 0 ) write( STDOUT, nml = held_suarez_1994_nml )
!!$ end if
! 係数の設定
! Configure coefficients
!
Kappa = GasRDry / CpDry
P0 = 1000.0e2_DP
DelTempY = 60.0_DP
DelPotTempZ = 10.0_DP
day_seconds = 86400.0_DP
SigmaB = 0.7_DP
kf = 1.0_DP / day_seconds
ka = 1.0_DP / ( 40.0_DP * day_seconds )
ks = 1.0_DP / ( 4.0_DP * day_seconds )
allocate( z_kv (1:kmax) )
z_kv = kf * max( 0.0_DP, ( z_Sigma - SigmaB ) / ( 1.0_DP - SigmaB ) )
allocate( yz_kt (1:jmax, 1:kmax) )
do k = 1, kmax
do j = 1, jmax
yz_kt(j,k) = ka + ( ks - ka ) * max( 0.0_DP, ( z_Sigma(k) - SigmaB ) / ( 1.0_DP - SigmaB ) ) * cos( y_Lat(j) ) ** 4
end do
end do
! ヒストリデータ出力のためのへの変数登録
! Register of variables for history data output
!
call HistoryAutoAddVariable( 'DUDtHs94', (/ 'lon ', 'lat ', 'sig ', 'time' /), 'eastward wind tendency', 'm s-2' )
call HistoryAutoAddVariable( 'DVDtHs94', (/ 'lon ', 'lat ', 'sig ', 'time' /), 'northward wind tendency', 'm s-2' )
call HistoryAutoAddVariable( 'DTempDtHs94', (/ 'lon ', 'lat ', 'sig ', 'time' /), 'temperature tendency', 'K s-1' )
call HistoryAutoAddVariable( 'TempEQHs94', (/ 'lon ', 'lat ', 'sig ', 'time' /), 'equilibrium temperature', 'K' )
! 印字 ; Print
!
call MessageNotify( 'M', module_name, '----- Initialization Messages -----' )
call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) )
held_suarez_1994_inited = .true.
end subroutine Hs94Init
| Subroutine : |
依存モジュールの初期化チェック
Check initialization of dependency modules
subroutine InitCheck
!
! 依存モジュールの初期化チェック
!
! Check initialization of dependency modules
! モジュール引用 ; USE statements
!
! NAMELIST ファイル入力に関するユーティリティ
! Utilities for NAMELIST file input
!
use namelist_util, only: namelist_util_inited
! 格子点設定
! Grid points settings
!
use gridset, only: gridset_inited
! 物理定数設定
! Physical constants settings
!
use constants, only: constants_inited
! 座標データ設定
! Axes data settings
!
use axesset, only: axesset_inited
! 時刻管理
! Time control
!
use timeset, only: timeset_inited
! 実行文 ; Executable statement
!
if ( .not. namelist_util_inited ) call MessageNotify( 'E', module_name, '"namelist_util" module is not initialized.' )
if ( .not. gridset_inited ) call MessageNotify( 'E', module_name, '"gridset" module is not initialized.' )
if ( .not. constants_inited ) call MessageNotify( 'E', module_name, '"constants" module is not initialized.' )
if ( .not. axesset_inited ) call MessageNotify( 'E', module_name, '"axesset" module is not initialized.' )
if ( .not. timeset_inited ) call MessageNotify( 'E', module_name, '"timeset" module is not initialized.' )
end subroutine InitCheck
| Variable : | |||
| Kappa : | real(DP)
|
| Constant : | |||
| module_name = ‘held_suarez_1994‘ : | character(*), parameter
|
| Constant : | |||
| version = ’$Name: dcpam5-20100224 $’ // ’$Id: held_suarez_1994.f90,v 1.9 2009-03-16 05:27:28 morikawa Exp $’ : | character(*), parameter
|