| Class | radiation_short_income_sr |
| In: |
radiation/radiation_short_income_sr.f90
|
Note that Japanese and English are described in parallel.
同期回転惑星を想定した短波入射 (太陽入射) を計算します.
Calculate short wave (insolation) incoming radiation assuming a planet with synchronized rotation.
| ShortIncoming : | 短波入射 (太陽入射) の計算 |
| ———— : | ———— |
| ShortIncoming : | Calculate short wave (insolation) incoming radiation. |
| Subroutine : | |||
| xy_IncomRadSFlux(0:imax-1, 1:jmax) : | real(DP), intent(out)
| ||
| xy_InAngle(0:imax-1, 1:jmax) : | real(DP), intent(out)
|
短波入射 (太陽入射) を計算します. 同期回転惑星を想定しており, 常に経度 90 度で最大入射, 経度 270 度で最小入射となっています.
Calculate short wave (insolation) incoming radiation. A planet with synchronized rotation are assumed. Incoming is max at latitude 90 deg., and min at 270 deg. constantly.
subroutine ShortIncoming( xy_IncomRadSFlux, xy_InAngle )
!
! 短波入射 (太陽入射) を計算します.
! 同期回転惑星を想定しており,
! 常に経度 90 度で最大入射, 経度 270 度で最小入射となっています.
!
! Calculate short wave (insolation) incoming radiation.
! A planet with synchronized rotation are assumed.
! Incoming is max at latitude 90 deg., and min at 270 deg. constantly.
!
! モジュール引用 ; USE statements
!
! 座標データ設定
! Axes data settings
!
use axesset, only: x_Lon, y_Lat ! $ \varphi $ [rad.] . 緯度. Latitude
! 物理定数設定
! Physical constants settings
!
use constants, only: PI ! $ \pi $ .
! 円周率. Circular constant
! 宣言文 ; Declaration statements
!
implicit none
real(DP), intent(out):: xy_IncomRadSFlux (0:imax-1, 1:jmax)
! 短波 (日射) フラックス.
! Short wave (insolation) flux
real(DP), intent(out):: xy_InAngle (0:imax-1, 1:jmax)
! sec (入射角).
! sec (angle of incidence)
! 作業変数
! Work variables
!
real(DP):: CosZet
real(DP):: AngleMaxLon
integer:: i ! 経度方向に回る DO ループ用作業変数
! Work variables for DO loop in longitude
integer:: j ! 緯度方向に回る DO ループ用作業変数
! Work variables for DO loop in latitude
! 実行文 ; Executable statement
!
! 初期化
! Initialization
!
if ( .not. radiation_short_income_inited ) call ShtIncomeInit
! 経度 90 度が最大入射, 経度 270 度が最小入射
! Incoming is max at latitude 90 deg., and min at 270 deg.
!
do i = 0, imax - 1
do j = 1, jmax
AngleMaxLon = - PI / 2.0_DP + x_Lon( i )
CosZet = cos( y_Lat(j) ) * cos( AngleMaxLon )
if ( CosZet > 0.0_DP ) then
xy_IncomRadSFlux( i,j ) = - SolarConst * CosZet
xy_InAngle( i,j ) = 1.0_DP / CosZet
else
xy_IncomRadSFlux( i,j ) = 0.
xy_InAngle( i,j ) = 0.
end if
end do
end do
end subroutine ShortIncoming
| Variable : | |||
| radiation_short_income_inited = .false. : | logical, save, public
|
| 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
| Subroutine : |
radiation_short_income_sr モジュールの初期化を行います. NAMELIST#radiation_short_income_nml の読み込みはこの手続きで行われます.
"radiation_short_income_sr" module is initialized. "NAMELIST#radiation_short_income_nml" is loaded in this procedure.
This procedure input/output NAMELIST#radiation_short_income_nml .
subroutine ShtIncomeInit
!
! radiation_short_income_sr モジュールの初期化を行います.
! NAMELIST#radiation_short_income_nml の読み込みはこの手続きで行われます.
!
! "radiation_short_income_sr" module is initialized.
! "NAMELIST#radiation_short_income_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
! ヒストリデータ出力
! 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 /radiation_short_income_nml/ SolarConst
!
! デフォルト値については初期化手続 "radiation_short_income#ShtIncomeInit"
! のソースコードを参照のこと.
!
! Refer to source codes in the initialization procedure
! "radiation_short_income#ShtIncomeInit" for the default values.
!
! 実行文 ; Executable statement
!
if ( radiation_short_income_inited ) return
call InitCheck
! デフォルト値の設定
! Default values settings
!
SolarConst = 1380.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 = radiation_short_income_nml, iostat = iostat_nml ) ! (out)
close( unit_nml )
call NmlutilMsg( iostat_nml, module_name ) ! (in)
end if
! 保存用の変数の割り付け
! Allocate variables for saving
!
! ヒストリデータ出力のためのへの変数登録
! Register of variables for history data output
!
!!$ call HistoryAutoAddVariable( 'xxxxx' , &
!!$ & (/ 'lon ', 'lat ', 'sig', 'time'/), &
!!$ & 'xxxx', 'W m-2' )
! 印字 ; Print
!
call MessageNotify( 'M', module_name, '----- Initialization Messages -----' )
call MessageNotify( 'M', module_name, 'ShortIncomming:' )
call MessageNotify( 'M', module_name, ' SolarConst = %f', d = (/ SolarConst /) )
call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) )
radiation_short_income_inited = .true.
end subroutine ShtIncomeInit
| Constant : | |||
| module_name = ‘radiation_short_income_sr‘ : | character(*), parameter
|
| Constant : | |||
| version = ’$Name: dcpam5-20100224 $’ // ’$Id: radiation_short_income_sr.f90,v 1.4 2009-08-09 12:30:52 yot Exp $’ : | character(*), parameter
|
| Variable : | |||
| xy_InAngle(:,:) : | real(DP), allocatable, save
|
| Variable : | |||
| xy_IncomRadSFlux(:,:) : | real(DP), allocatable, save
|