| Class | Radiation |
| In: |
physics/radiation.f90
|
モデルの放射過程を計算するためのパッケージ型モジュール 具体的には以下の項を計算するための関数を格納する.
* 一様冷却
| Subroutine : | |
| cfgfile : | character(*), intent(in) |
NAMELIST から放射強制の設定を取得
This procedure input/output NAMELIST#radiation .
subroutine Radiation_init(cfgfile)
!
!NAMELIST から放射強制の設定を取得
!
!暗黙の型宣言禁止
implicit none
!入力変数
character(*), intent(in) :: cfgfile
real(8) :: RadHeightUp !放射強制を与える鉛直領域の上限
real(8) :: RadHeightDown !放射強制を与える鉛直領域の下限
integer :: k !ループ変数
! NAMELIST から情報を取得
NAMELIST /radiation/ RadHeatRate, RadHeightUp, RadHeightDown
open (10, FILE=cfgfile)
read(10, NML=radiation)
close(10)
allocate( xz_RadHeight(DimXMin:DimXMax, DimZMin:DimZMax) )
! 放射強制が存在する領域の設定
! RadHeightDown < s_Z < RadHeightUp の範囲では値が 1 となる
! 係数を用意する.
!
xz_RadHeight = 1.0d0
do k = DimZMin, DimZMax
if ( s_Z(k) <= RadHeightDown ) then
xz_RadHeight(:,k) = 0.0d0
elseif( s_Z(k) >= RadHeightUp ) then
xz_RadHeight(:,k) = 0.0d0
end if
end do
call MessageNotify( "M", "Radiation", "RadHeatRate = %f", d=(/RadHeatRate/))
call MessageNotify( "M", "Radiation", "RadHeightUp = %f", d=(/RadHeightUP/))
call MessageNotify( "M", "Radiation", "RadHeightDown= %f", d=(/RadHeightDown/))
end subroutine Radiation_init
| Function : | |||
| xz_NewtonCool(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8)
| ||
| xz_PotTemp(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in)
|
温位の放射強制項. 地表面から RadHeight で指定された高度までの間で一様放射冷却を与える.
function xz_NewtonCool(xz_PotTemp)
!
! 温位の放射強制項.
! 地表面から RadHeight で指定された高度までの間で一様放射冷却を与える.
!
!暗黙の型宣言を禁止
implicit none
!変数定義
real(8), intent(in) :: xz_PotTemp(DimXMin:DimXMax, DimZMin:DimZMax)
!温位の擾乱場
real(8) :: xz_NewtonCool(DimXMin:DimXMax, DimZMin:DimZMax)
!放射強制
real(8) :: xz_PotTempMean(DimXMin:DimXMax, DimZMin:DimZMax)
!温位擾乱の東西平均
real(8) :: EFTime !放射緩和時間
integer :: k
!放射緩和時間は 5 日
EFTime = 5.0d0 * DayTime
do k = DimZMin, DimZMax
xz_PotTempMean(:,k) = sum( xz_PotTemp(:,k) ) / real(DimXMax - DimXMin + 1)
end do
xz_NewtonCool = - xz_PotTempMean / EFTime
call StorePotTempDamp( xz_NewtonCool )
end function xz_NewtonCool
| Function : | |||
| xz_RadHeatConst(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8)
| ||
| xz_Exner(DimXMin:DimXMax, DimZMin:DimZMax) : | real(8), intent(in)
|
温位の放射強制項. 地表面から RadHeight で指定された高度までの間で一様放射冷却を与える.
function xz_RadHeatConst(xz_Exner)
!
! 温位の放射強制項.
! 地表面から RadHeight で指定された高度までの間で一様放射冷却を与える.
!
!暗黙の型宣言を禁止
implicit none
!変数定義
real(8), intent(in) :: xz_Exner(DimXMin:DimXMax, DimZMin:DimZMax)
!エクスナー関数の擾乱場
real(8) :: xz_RadHeatConst(DimXMin:DimXMax, DimZMin:DimZMax)
!放射強制
xz_RadHeatConst = xz_RadHeight * RadHeatRate / ( ( xz_ExnerBasicZ + xz_Exner ) * DayTime )
call StorePotTempRad( xz_RadHeatConst )
end function xz_RadHeatConst