Class | rad_Earth_V2 |
In: |
radiation/rad_Earth_V2.f90
|
Note that Japanese and English are described in parallel.
地球大気向け放射モデル.
This is a radiation model for the Earth‘s atmospehre.
Radiation in the wavenumber range from 0 to 3000 cm-1 is calculated in the routine of long wave radiation. Radiation in the wavenumber range from 1000 to 57143 cm-1 (0.175 to 10 micron) is calculated in the routine of shortwave radiation.
RadEarthV2Flux : | 放射フラックスの計算 |
———— : | ———— |
RadEarthV2Flux : | Calculate radiation flux |
Subroutine : | |||
xy_SurfAlbedo(0:imax-1, 1:jmax) : | real(DP), intent(in ) | ||
xyz_Press(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(in ) | ||
xyr_Press(0:imax-1, 1:jmax, 0: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 )
| ||
xyz_QH2OLiq(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(in ) | ||
xyz_QH2OSol(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(in ) | ||
xyz_CloudCover(0:imax-1, 1:jmax, 1:kmax) : | 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 RadEarthV2Flux( xy_SurfAlbedo, xyz_Press, xyr_Press, xyz_Temp, xyz_QH2OVap, xyz_QH2OLiq, xyz_QH2OSol, xyz_CloudCover, xy_SurfTemp, xyr_RadSUwFlux, xyr_RadSDwFlux, xyr_RadLUwFlux, xyr_RadLDwFlux, xyra_DelRadLUwFlux, xyra_DelRadLDwFlux ) ! USE statements ! real(DP), intent(in ) :: xy_SurfAlbedo (0:imax-1, 1:jmax) real(DP), intent(in ) :: xyz_Press (0:imax-1, 1:jmax, 1:kmax) real(DP), intent(in ) :: xyr_Press (0:imax-1, 1:jmax, 0: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) ! $ q $ . 混合比. Mass mixing ratio of constituents (1) real(DP), intent(in ) :: xyz_QH2OLiq (0:imax-1, 1:jmax, 1:kmax) real(DP), intent(in ) :: xyz_QH2OSol (0:imax-1, 1:jmax, 1:kmax) real(DP), intent(in ) :: xyz_CloudCover (0:imax-1, 1:jmax, 1:kmax) 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) !!$ real(DP):: xyr_RadLUwFlux (0:imax-1, 1:jmax, 0:kmax) !!$ real(DP):: xyr_RadLDwFlux (0:imax-1, 1:jmax, 0:kmax) !!$ real(DP):: xyra_DelRadLUwFlux(0:imax-1, 1:jmax, 0:kmax, 0:1) !!$ real(DP):: xyra_DelRadLDwFlux(0:imax-1, 1:jmax, 0:kmax, 0:1) ! Work variables ! real(DP) :: xyz_DelAtmMass (0:imax-1, 1:jmax, 1:kmax) integer :: k ! 初期化確認 ! Initialization check ! if ( .not. rad_Earth_V2_inited ) then call MessageNotify( 'E', module_name, 'This module has not been initialized.' ) end if do k = 1, kmax xyz_DelAtmMass(:,:,k) = ( xyr_Press(:,:,k-1) - xyr_Press(:,:,k ) ) / Grav end do call RadEarthV2FluxCore( xy_SurfAlbedo, xyz_DelAtmMass, xyz_Press, xyz_Temp, xyz_QH2OVap, xyz_QH2OLiq, xyz_QH2OSol, xyz_CloudCover, xy_SurfTemp, xyr_RadSUwFlux, xyr_RadSDwFlux, xyr_RadLUwFlux, xyr_RadLDwFlux, xyra_DelRadLUwFlux, xyra_DelRadLDwFlux ) end subroutine RadEarthV2Flux
Subroutine : | |
FlagSnow : | logical, intent(in) |
This procedure input/output NAMELIST#rad_Earth_V2_nml .
subroutine RadEarthV2Init( FlagSnow ) ! ファイル入出力補助 ! File I/O support ! use dc_iounit, only: FileOpen ! ヒストリデータ出力 ! History data output ! use gtool_historyauto, only: HistoryAutoAddVariable ! NAMELIST ファイル入力に関するユーティリティ ! Utilities for NAMELIST file input ! use namelist_util, only: namelist_filename, NmlutilMsg, NmlutilAryValid ! O3 分布の設定 ! Set O3 distribution ! use set_o3, only : SetO3Init ! 地球大気向け短波放射モデル Ver. 2.1 ! short wave radiation model for the Earth's atmosphere Ver. 2.1 ! use rad_Earth_SW_V2_1, only: RadEarthSWV21Init !!$ ! 地球大気向け短波放射モデル Ver. 2.2 !!$ ! short wave radiation model for the Earth's atmosphere Ver. 2.2 !!$ ! !!$ use rad_Earth_SW_V2_2, only: RadEarthSWV22Init ! 地球大気向け長波放射モデル Ver. 2.3 ! long wave radiation model for the Earth's atmosphere Ver. 2.3 ! use rad_Earth_LW_V2_3, only : RadEarthLWV23Init ! 宣言文 ; Declaration statements ! logical, intent(in) :: FlagSnow 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/ VMRCO2 !!$ & SWVer, LWVer ! ! デフォルト値については初期化手続 "rad_Earth_V2#RadEarthV2Init" ! のソースコードを参照のこと. ! ! Refer to source codes in the initialization procedure ! "rad_Earth_V2#RadEarthV2Init" for the default values. ! if ( rad_Earth_V2_inited ) return ! Molecular weights of atmosphere, CO2, and H2O ! MeanMolWeight = 28.0d-3 CO2MolWeight = 44.0d-3 H2OMolWeight = 18.0d-3 ! デフォルト値の設定 ! Default values settings ! VMRCO2 = 382.0d-6 !!$ SWVer = 1 !!$ LWVer = 3 ! 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 = rad_Earth_V2_nml, iostat = iostat_nml ) ! (out) close( unit_nml ) call NmlutilMsg( iostat_nml, module_name ) ! (in) end if ! Allocate a local variable for O3 distribution ! allocate( xyz_QO3(0:imax-1, 1:jmax, 1:kmax) ) ! Initialization of modules used in this module ! ! O3 分布の設定 ! Set O3 distribution ! call SetO3Init ! 地球大気向け短波放射モデル Ver. 2.1 ! short wave radiation model for the Earth's atmosphere Ver. 2.1 ! call RadEarthSWV21Init !!$ ! 地球大気向け短波放射モデル Ver. 2.2 !!$ ! short wave radiation model for the Earth's atmosphere Ver. 2.2 !!$ ! !!$ call RadEarthSWV22Init( & !!$ & FlagSnow & !!$ & ) ! 地球大気向け長波放射モデル Ver. 2.3 ! long wave radiation model for the Earth's atmosphere Ver. 2.3 ! call RadEarthLWV23Init( FlagSnow ) ! ヒストリデータ出力のためのへの変数登録 ! Register of variables for history data output ! call HistoryAutoAddVariable( 'CloudCoverforRad', (/ 'lon ', 'lat ', 'sig ', 'time' /), 'cloud cover', '1' ) ! 印字 ; Print ! call MessageNotify( 'M', module_name, '----- Initialization Messages -----' ) call MessageNotify( 'M', module_name, ' VMRCO2 = %f', d = (/ VMRCO2 /) ) !!$ 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_Earth_V2_inited = .true. end subroutine RadEarthV2Init
Subroutine : | |||
xy_SurfAlbedo(0:imax-1, 1:jmax) : | real(DP), intent(in ) | ||
xyz_DelAtmMass(0:imax-1, 1:jmax, 1: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 )
| ||
xyz_QH2OLiq(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(in ) | ||
xyz_QH2OSol(0:imax-1, 1:jmax, 1:kmax) : | real(DP), intent(in ) | ||
xyz_CloudCover(0:imax-1, 1:jmax, 1:kmax) : | 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 RadEarthV2FluxCore( xy_SurfAlbedo, xyz_DelAtmMass, xyz_Press, xyz_Temp, xyz_QH2OVap, xyz_QH2OLiq, xyz_QH2OSol, xyz_CloudCover, xy_SurfTemp, xyr_RadSUwFlux, xyr_RadSDwFlux, xyr_RadLUwFlux, xyr_RadLDwFlux, xyra_DelRadLUwFlux, xyra_DelRadLDwFlux ) ! USE statements ! ! ヒストリデータ出力 ! History data output ! use gtool_historyauto, only: HistoryAutoPut ! 時刻管理 ! Time control ! use timeset, only: TimeN ! ステップ $ t $ の時刻. ! Time of step $ t $. ! O3 分布の設定 ! Set O3 distribution ! use set_o3, only: SetO3 !!$ ! dcpam 地球大気向け短波放射モデル Ver. 2 !!$ ! dcpam short wave radiation model for the Earth's atmosphere Ver. 2 !!$ ! !!$ use radiation_dcpam_E_SW_V2, only: RadiationDcpamESWV2Flux ! 地球大気向け短波放射モデル Ver. 2.1 ! short wave radiation model for the Earth's atmosphere Ver. 2.1 ! use rad_Earth_SW_V2_1, only: RadEarthSWV21Flux !!$ !!$ ! 地球大気向け短波放射モデル Ver. 2.2 !!$ ! short wave radiation model for the Earth's atmosphere Ver. 2.2 !!$ ! !!$ use rad_Earth_SW_V2_2, only: RadEarthSWV22Flux !!$ ! dcpam 地球大気向け長波放射モデル Ver. 2 !!$ ! dcpam long wave radiation model for the Earth's atmosphere Ver. 2 !!$ ! !!$ use radiation_dcpam_E_LW_V2, only : RadiationDcpamELWV2Flux !!$ !!$ ! dcpam 地球大気向け長波放射モデル Ver. 2.1 !!$ ! dcpam long wave radiation model for the Earth's atmosphere Ver. 2.1 !!$ ! !!$ use radiation_dcpam_E_LW_V2_1, only : RadiationDcpamELWV21Flux !!$ !!$ ! dcpam 地球大気向け長波放射モデル Ver. 2.2 !!$ ! dcpam long wave radiation model for the Earth's atmosphere Ver. 2.2 !!$ ! !!$ use radiation_dcpam_E_LW_V2_2, only : RadiationDcpamELWV22Flux ! 地球大気向け長波放射モデル Ver. 2.3 ! long wave radiation model for the Earth's atmosphere Ver. 2.3 ! use rad_Earth_LW_V2_3, only : RadEarthLWV23Flux real(DP), intent(in ) :: xy_SurfAlbedo (0:imax-1, 1:jmax) real(DP), intent(in ) :: xyz_DelAtmMass (0:imax-1, 1:jmax, 1: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) ! $ q $ . 混合比. Mass mixing ratio of constituents (1) real(DP), intent(in ) :: xyz_QH2OLiq (0:imax-1, 1:jmax, 1:kmax) real(DP), intent(in ) :: xyz_QH2OSol (0:imax-1, 1:jmax, 1:kmax) real(DP), intent(in ) :: xyz_CloudCover (0:imax-1, 1:jmax, 1:kmax) 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) :: xyz_QCO2 (0:imax-1, 1:jmax, 1:kmax) real(DP) :: xyz_DelCO2Mass (0:imax-1, 1:jmax, 1:kmax) real(DP) :: xyz_DelH2OVapMass(0:imax-1, 1:jmax, 1:kmax) real(DP) :: xyz_DelH2OLiqMass(0:imax-1, 1:jmax, 1:kmax) real(DP) :: xyz_DelH2OSolMass(0:imax-1, 1:jmax, 1:kmax) real(DP) :: xyz_DelO3Mass (0:imax-1, 1:jmax, 1:kmax) ! 初期化確認 ! Initialization check ! if ( .not. rad_Earth_V2_inited ) then call MessageNotify( 'E', module_name, 'This module has not been initialized.' ) end if xyz_QCO2 = VMRCO2 * CO2MolWeight / MeanMolWeight ! O3 分布の設定 ! Setting of O3 distribution ! call SetO3( xyz_Press, xyz_QO3 ) xyz_DelCO2Mass = xyz_DelAtmMass * xyz_QCO2 xyz_DelH2OVapMass = xyz_DelAtmMass * xyz_QH2OVap xyz_DelH2OLiqMass = xyz_DelAtmMass * xyz_QH2OLiq xyz_DelH2OSolMass = xyz_DelAtmMass * xyz_QH2OSol xyz_DelO3Mass = xyz_DelAtmMass * xyz_QO3 !!$ select case ( SWVer ) !!$ case ( 0 ) !!$ ! dcpam 地球大気向け短波放射モデル Ver. 2 !!$ ! dcpam short wave radiation model for the Earth's atmosphere Ver. 2 !!$ ! !!$ call RadiationDcpamESWV2Flux( & !!$ & xy_SurfAlbedo, & !!$ & xyz_Press, xyr_Press, xyz_Temp, xyz_QH2OVap, xyz_QO3, xyz_Height, & !!$ & xyr_RadSFlux & !!$ & ) !!$ case ( 1 ) ! 地球大気向け短波放射モデル Ver. 2.1 ! short wave radiation model for the Earth's atmosphere Ver. 2.1 ! call RadEarthSWV21Flux( xy_SurfAlbedo, xyz_DelAtmMass, xyz_DelH2OVapMass, xyz_DelH2OLiqMass, xyz_DelH2OSolMass, xyz_DelO3Mass, xyz_Press, xyz_Temp, xyz_CloudCover, xyr_RadSUwFlux, xyr_RadSDwFlux ) !!$ case ( 2 ) !!$ ! 地球大気向け短波放射モデル Ver. 2.2 !!$ ! short wave radiation model for the Earth's atmosphere Ver. 2.2 !!$ ! !!$ call RadEarthSWV22Flux( & !!$ & xy_SurfAlbedo, & !!$ & xyz_DelAtmMass, & !!$ & xyz_DelH2OVapMass, xyz_DelH2OLiqMass, xyz_DelH2OSolMass, & !!$ & xyz_DelO3Mass, & !!$ & xyz_Press, xyz_Temp, & !!$ & xyz_CloudCover, & !!$ & xyr_RadSUwFlux, xyr_RadSDwFlux & !!$ & ) !!$ case default !!$ call MessageNotify( 'E', module_name, 'SW model version %d is not supported.', i = (/ LWVer /) ) !!$ end select !!$ select case ( LWVer ) !!$ case ( 0 ) !!$ call RadiationDcpamELWV2Flux( & !!$ & xyz_Press, xyr_Press, xyz_Temp, xy_SurfTemp, & ! (in ) !!$ & xyz_QH2OVap, xyz_QO3, & ! (in ) !!$ & xyr_RadLFlux, xyra_DelRadLFlux & ! (out) !!$ & ) !!$ case ( 1 ) !!$ call RadiationDcpamELWV21Flux( & !!$ & xyz_Press, xyr_Press, xyz_Temp, xy_SurfTemp, & ! (in ) !!$ & xyz_QH2OVap, xyz_QO3, & ! (in ) !!$ & xyr_RadLFlux, xyra_DelRadLFlux & ! (out) !!$ & ) !!$ case ( 2 ) !!$ call RadiationDcpamELWV22Flux( & !!$ & xyz_Press, xyr_Press, xyz_Temp, xy_SurfTemp, & ! (in ) !!$ & xyz_QH2OVap, xyz_QO3, & ! (in ) !!$ & xyr_RadLFlux, xyra_DelRadLFlux & ! (out) !!$ & ) !!$ case ( 3 ) call RadEarthLWV23Flux( xyz_DelCO2Mass, xyz_DelH2OVapMass, xyz_DelH2OLiqMass, xyz_DelH2OSolMass, xyz_DelO3Mass, xyz_Press, xyz_Temp, xy_SurfTemp, xyz_QCO2, xyz_QH2OVap, xyz_CloudCover, xyr_RadLUwFlux, xyr_RadLDwFlux, xyra_DelRadLUwFlux, xyra_DelRadLDwFlux ) !!$ case default !!$ call MessageNotify( 'E', module_name, 'LW model version %d is not supported.', i = (/ LWVer /) ) !!$ end select ! Output variables ! call HistoryAutoPut( TimeN, 'CloudCoverforRad', xyz_CloudCover ) end subroutine RadEarthV2FluxCore
Constant : | |||
module_name = ‘rad_Earth_V2‘ : | character(*), parameter
|
Variable : | |||
rad_Earth_V2_inited = .false. : | logical, save
|
Constant : | |||
version = ’$Name: dcpam5-20121110 $’ // ’$Id: rad_Earth_V2.f90,v 1.6 2012-09-08 15:25:05 yot Exp $’ : | character(*), parameter
|