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