| Class | rad_none |
| In: |
radiation/rad_none.f90
|
Note that Japanese and English are described in parallel.
| RadEarthV2Flux : | 放射フラックスの計算 |
| ———— : | ———— |
| RadEarthV2Flux : | Calculate radiation flux |
| Subroutine : | |
| xy_SurfAlbedo(0:imax-1, 1:jmax) : | real(DP), intent(in ) |
| xy_SurfTemp(0:imax-1, 1:jmax) : | real(DP), intent(in ) |
| xyr_RadSUwFlux(0:imax-1, 1:jmax, 0:kmax) : | real(DP), intent(out) |
| xyr_RadSDwFlux(0:imax-1, 1:jmax, 0:kmax) : | real(DP), intent(out) |
| xyr_RadLUwFlux(0:imax-1, 1:jmax, 0:kmax) : | real(DP), intent(out) |
| xyr_RadLDwFlux(0:imax-1, 1:jmax, 0:kmax) : | real(DP), intent(out) |
| xyra_DelRadLUwFlux(0:imax-1, 1:jmax, 0:kmax, 0:1) : | real(DP), intent(out) |
| xyra_DelRadLDwFlux(0:imax-1, 1:jmax, 0:kmax, 0:1) : | real(DP), intent(out) |
subroutine RadNoneFlux( xy_SurfAlbedo, xy_SurfTemp, xyr_RadSUwFlux, xyr_RadSDwFlux, xyr_RadLUwFlux, xyr_RadLDwFlux, xyra_DelRadLUwFlux, xyra_DelRadLDwFlux )
! USE statements
!
! 物理・数学定数設定
! Physical and mathematical constants settings
!
use constants0, only: StB ! $ \sigma_{SB} $ .
! ステファンボルツマン定数.
! Stefan-Boltzmann constant
! 太陽放射フラックスの設定
! Set solar constant
!
use set_solarconst, only : SetSolarConst
! 短波入射 (太陽入射)
! Short wave (insolation) incoming
!
use rad_short_income, only : RadShortIncome
real(DP), intent(in ) :: xy_SurfAlbedo (0:imax-1, 1:jmax)
real(DP), intent(in ) :: xy_SurfTemp (0:imax-1, 1:jmax)
real(DP), intent(out) :: xyr_RadSUwFlux (0:imax-1, 1:jmax, 0:kmax)
real(DP), intent(out) :: xyr_RadSDwFlux (0:imax-1, 1:jmax, 0:kmax)
real(DP), intent(out) :: xyr_RadLUwFlux (0:imax-1, 1:jmax, 0:kmax)
real(DP), intent(out) :: xyr_RadLDwFlux (0:imax-1, 1:jmax, 0:kmax)
real(DP), intent(out) :: xyra_DelRadLUwFlux(0:imax-1, 1:jmax, 0:kmax, 0:1)
real(DP), intent(out) :: xyra_DelRadLDwFlux(0:imax-1, 1:jmax, 0:kmax, 0:1)
! Work variables
!
real(DP) :: SolarConst
real(DP) :: DistFromStarScld
real(DP) :: DiurnalMeanFactor
real(DP) :: SolarFluxTOA
real(DP) :: xy_CosSZA(0:imax-1, 1:jmax)
integer :: i
integer :: j
integer :: k
! 初期化確認
! Initialization check
!
if ( .not. rad_none_inited ) then
call MessageNotify( 'E', module_name, 'This module has not been initialized.' )
end if
! 太陽放射フラックスの設定
! Set solar constant
!
call SetSolarConst( SolarConst )
! 短波入射の計算
! Calculate short wave (insolation) incoming radiation
!
call RadShortIncome( xy_CosZet = xy_CosSZA, DistFromStarScld = DistFromStarScld, DiurnalMeanFactor = DiurnalMeanFactor )
SolarFluxTOA = SolarConst / DistFromStarScld**2 * DiurnalMeanFactor
do k = 0, kmax
do j = 1, jmax
do i = 0, imax-1
if ( xy_CosSZA(i,j) >= 0.0_DP ) then
xyr_RadSDwFlux(i,j,k) = SolarFluxTOA * xy_CosSZA(i,j)
else
xyr_RadSDwFlux(i,j,k) = 0.0_DP
end if
end do
end do
end do
do k = 0, kmax
xyr_RadSUwFlux(:,:,k) = xy_SurfAlbedo * xyr_RadSDwFlux(:,:,k)
end do
do k = 0, kmax
xyr_RadLUwFlux(:,:,k) = StB * xy_SurfTemp**4
end do
xyr_RadLDwFlux = 0.0_DP
xyra_DelRadLUwFlux = 0.0_DP
xyra_DelRadLDwFlux = 0.0_DP
! Output variables
!
end subroutine RadNoneFlux
| Subroutine : |
subroutine RadNoneInit
! ファイル入出力補助
! File I/O support
!
!!$ use dc_iounit, only: FileOpen
! NAMELIST ファイル入力に関するユーティリティ
! Utilities for NAMELIST file input
!
!!$ use namelist_util, only: namelist_filename, NmlutilMsg, NmlutilAryValid
! 太陽放射フラックスの設定
! Set solar constant
!
use set_solarconst, only : SetSolarConstInit
! 短波入射 (太陽入射)
! Short wave (insolation) incoming
!
use rad_short_income, only : RadShortIncomeInit
! 宣言文 ; Declaration statements
!
!!$ integer:: unit_nml ! NAMELIST ファイルオープン用装置番号.
!!$ ! Unit number for NAMELIST file open
!!$ integer:: iostat_nml ! NAMELIST 読み込み時の IOSTAT.
!!$ ! IOSTAT of NAMELIST read
! NAMELIST 変数群
! NAMELIST group name
!
!!$ namelist /rad_Earth_V2_nml/ &
!!$ & CloudIceREffMethod, &
!!$ & CloudWatREff, &
!!$ & CloudIceREff
!!$ & SWVer, LWVer
!
! デフォルト値については初期化手続 "rad_Earth_V2#RadEarthV2Init"
! のソースコードを参照のこと.
!
! Refer to source codes in the initialization procedure
! "rad_Earth_V2#RadEarthV2Init" for the default values.
!
if ( rad_none_inited ) return
! デフォルト値の設定
! Default values settings
!
!!$ SWVer = 1
!!$ LWVer = 3
! 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 = rad_Earth_V2_nml, & ! (out)
!!$ & iostat = iostat_nml ) ! (out)
!!$ close( unit_nml )
!!$
!!$ call NmlutilMsg( iostat_nml, module_name ) ! (in)
!!$ end if
! 太陽放射フラックスの設定
! Set solar constant
!
call SetSolarConstInit
! 短波入射 (太陽入射)
! Short wave (insolation) incoming
!
call RadShortIncomeInit
! 印字 ; Print
!
call MessageNotify( 'M', module_name, '----- Initialization Messages -----' )
!!$ call MessageNotify( 'M', module_name, 'SWVer = %d', i = (/ SWVer /) )
!!$ call MessageNotify( 'M', module_name, 'LWVer = %d', i = (/ LWVer /) )
call MessageNotify( 'M', module_name, '-- version = %c', c1 = trim(version) )
rad_none_inited = .true.
end subroutine RadNoneInit
| Constant : | |||
| version = ’$Name: $’ // ’$Id: rad_none.f90,v 1.1 2015/01/31 06:17:21 yot Exp $’ : | character(*), parameter
|