Subroutine : |
|
xy_SurfAlbedo( 0:imax-1, 1:jmax ) : | real(DP), intent(in )
|
xy_SurfTemp( 0:imax-1, 1:jmax ) : | real(DP), intent(in )
|
xyr_Press( 0:imax-1, 1:jmax, 0:kmax ) : | real(DP), intent(in )
|
xyz_Press( 0:imax-1, 1:jmax, 1:kmax ) : | real(DP), intent(in )
|
xyz_Temp( 0:imax-1, 1:jmax, 1:kmax ) : | real(DP), intent(in )
|
xyz_QH2OVap( 0:imax-1, 1:jmax, 1:kmax ) : | 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 RadSimpleFlux( xy_SurfAlbedo, xy_SurfTemp, xyr_Press, xyz_Press, xyz_Temp, xyz_QH2OVap, xyr_RadSUwFlux, xyr_RadSDwFlux, xyr_RadLUwFlux, xyr_RadLDwFlux, xyra_DelRadLUwFlux, xyra_DelRadLDwFlux )
! USE statements
!
! 物理・数学定数設定
! Physical and mathematical constants settings
!
use constants0, only: PI ! $ \pi $ .
! 円周率. Circular constant
!
! Physical constants settings
!
use constants, only: Grav ! $ g $ [m s-2].
!
! Gravitational acceleration
! 座標データ設定
! Axes data settings
!
use axesset, only : y_Lat
! 短波入射 (太陽入射)
! Short wave (insolation) incoming
!
use rad_short_income, only : RadShortIncome
!
! Solve radiative transfer equation in two stream approximation
!
use rad_rte_two_stream_app, only: RadRTETwoStreamAppSW
! 長波フラックスの計算
! Calculate long wave flux
!
use rad_simple_LW, only : RadSimpleLWFlux
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(in ) :: xyr_Press ( 0:imax-1, 1:jmax, 0:kmax )
real(DP), intent(in ) :: xyz_Press ( 0:imax-1, 1:jmax, 1:kmax )
real(DP), intent(in ) :: xyz_Temp ( 0:imax-1, 1:jmax, 1:kmax )
real(DP), intent(in ) :: xyz_QH2OVap ( 0:imax-1, 1:jmax, 1:kmax )
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) :: DistFromStarScld
! Distance between the central star and the planet
real(DP) :: DiurnalMeanFactor
real(DP) :: SolarFluxTOA
!!$ real(DP) :: QeRatio
real(DP) :: xyz_SSA (0:imax-1, 1:jmax, 1:kmax)
real(DP) :: xyz_AF (0:imax-1, 1:jmax, 1:kmax)
real(DP) :: xy_InAngle (0:imax-1, 1:jmax)
real(DP) :: xyz_DelAtmMass(0:imax-1, 1:jmax, 1:kmax)
real(DP) :: xyr_OptDep (0:imax-1, 1:jmax, 0:kmax)
real(DP) :: xy_SurfEmis (0:imax-1, 1:jmax)
integer :: k
! 初期化
! Initialization
!
if ( .not. rad_Simple_inited ) call RadSimpleInit
! Short wave radiation
!
! 短波入射の計算
! Calculate short wave (insolation) incoming radiation
!
call RadShortIncome( xy_InAngle = xy_InAngle, DistFromStarScld = DistFromStarScld, DiurnalMeanFactor = DiurnalMeanFactor )
SolarFluxTOA = SolarConst / DistFromStarScld**2 * DiurnalMeanFactor
xyz_SSA = 1.0_DP - 1.0d-10
xyz_AF = 0.0_DP
xyr_OptDep = 0.0_DP
call RadRTETwoStreamAppSW( xyz_SSA, xyz_AF, xyr_OptDep, xy_SurfAlbedo, SolarFluxTOA, xy_InAngle, xyr_RadSUwFlux, xyr_RadSDwFlux )
! Long wave radiation
!
do k = 1, kmax
xyz_DelAtmMass(:,:,k) = ( xyr_Press(:,:,k-1) - xyr_Press(:,:,k ) ) / Grav
end do
xy_SurfEmis = 1.0_DP
call RadSimpleLWFlux( xy_SurfEmis, xyz_Temp, xyz_Press, xy_SurfTemp, xyz_DelAtmMass, xyz_QH2OVap, xyr_RadLUwFlux, xyr_RadLDwFlux, xyra_DelRadLUwFlux, xyra_DelRadLDwFlux )
end subroutine RadSimpleFlux