scale_atmos_phy_ch_rn222 Module Reference

module atmosphere / physics / chemistry / RN222 More...

Functions/Subroutines
subroutine, public	atmos_phy_ch_rn222_setup
	Setup. More...
subroutine, public	atmos_phy_ch_rn222_tendency (KA, KS, KE, IA, IS, IE, JA, JS, JE, QA_CH, DENS, QTRC, RHOQ_t)
	Chemistry Microphysics. More...

Variables
integer, public	atmos_phy_ch_rn222_ntracers = QA_CH
character(len=h_short), dimension(qa_ch), public	atmos_phy_ch_rn222_name = (/ "RN222" /)
character(len=h_mid), dimension(qa_ch), public	atmos_phy_ch_rn222_desc = (/ "Ratio of Rn222 to total mass" /)
character(len=h_short), dimension(qa_ch), public	atmos_phy_ch_rn222_unit = (/ "Bq/kg" /)

Detailed Description

module atmosphere / physics / chemistry / RN222

Description

Component for rn222 tracer

Author

Team SCALE

NAMELIST

• PARAM_ATMOS_PHY_CH_RN222

  name	type	default value	comment
  ATMOS_PHY_CH_RN222_HALF_LIFE	real(RP)	3.30048E+5_RP

  
  

History Output

No history output

Function/Subroutine Documentation

atmos_phy_ch_rn222_setup()

subroutine, public scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_setup

Setup.

Definition at line 59 of file scale_atmos_phy_ch_rn222.F90.

    use scale_prc, only: &
       prc_abort
    implicit none
 
    real(RP) :: ATMOS_PHY_CH_Rn222_half_life = 3.30048e+5_rp
 
    namelist / param_atmos_phy_ch_rn222 / &
       atmos_phy_ch_rn222_half_life
 
    integer :: ierr
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_CH_rn222_setup",*) 'Setup'
    log_info("ATMOS_PHY_CH_rn222_setup",*) 'rn222 process'
 
    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_ch_rn222,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_CH_rn222_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_CH_rn222_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_CH_RN222. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_phy_ch_rn222)
 
    atmos_phy_ch_rn222_decay_ratio = log(2.0_rp) / atmos_phy_ch_rn222_half_life
 
    log_newline
    log_info("ATMOS_PHY_CH_rn222_setup",*) 'Characteristics of Rn222'
    log_info_cont('(A,E16.6)')             'Half life   [s]   : ', atmos_phy_ch_rn222_half_life
    log_info_cont('(A,E16.6)')             'Decay ratio [1/s] : ', atmos_phy_ch_rn222_decay_ratio
 
    return

References scale_io::io_fid_conf, and scale_prc::prc_abort().

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ch_rn222_tendency()

subroutine, public scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_tendency	(	integer, intent(in)	KA,
		integer, intent(in)	KS,
		integer, intent(in)	KE,
		integer, intent(in)	IA,
		integer, intent(in)	IS,
		integer, intent(in)	IE,
		integer, intent(in)	JA,
		integer, intent(in)	JS,
		integer, intent(in)	JE,
		integer, intent(in)	QA_CH,
		real(rp), dimension (ka,ia,ja), intent(in)	DENS,
		real(rp), dimension (ka,ia,ja,qa_ch), intent(in)	QTRC,
		real(rp), dimension(ka,ia,ja,qa_ch), intent(inout)	RHOQ_t
	)

Chemistry Microphysics.

Definition at line 106 of file scale_atmos_phy_ch_rn222.F90.

    implicit none
 
    integer,  intent(in)    :: KA, KS, KE
    integer,  intent(in)    :: IA, IS, IE
    integer,  intent(in)    :: JA, JS, JE
    integer,  intent(in)    :: QA_CH
    real(RP), intent(in)    :: DENS  (KA,IA,JA)
    real(RP), intent(in)    :: QTRC  (KA,IA,JA,QA_CH)
    real(RP), intent(inout) :: RHOQ_t(KA,IA,JA,QA_CH)
 
    integer  :: k, i, j
    !---------------------------------------------------------------------------
 
    log_progress(*) 'atmosphere / physics / chemistry / Rn222'
 
    !--- Decay based on half life
    do j  = js, je
    do i  = is, ie
    do k  = ks, ke
       rhoq_t(k,i,j,i_ch_rn222) = rhoq_t(k,i,j,i_ch_rn222) &
                                - dens(k,i,j) * qtrc(k,i,j,1) * atmos_phy_ch_rn222_decay_ratio ! [Bq/m3/s]
    enddo
    enddo
    enddo
 
    return

References scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_tracer::k, scale_atmos_grid_cartesc_index::ke, and scale_atmos_grid_cartesc_index::ks.

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency().

Here is the caller graph for this function:

Variable Documentation

atmos_phy_ch_rn222_ntracers

integer, public scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_ntracers = QA_CH

Definition at line 36 of file scale_atmos_phy_ch_rn222.F90.

  integer,                public :: ATMOS_PHY_CH_rn222_ntracers = qa_ch

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup().

atmos_phy_ch_rn222_name

character(len=h_short), dimension(qa_ch), public scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_name = (/ "RN222" /)

Definition at line 38 of file scale_atmos_phy_ch_rn222.F90.

  character(len=H_SHORT), public :: ATMOS_PHY_CH_rn222_NAME(QA_CH) = (/ "RN222" /)

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup().

atmos_phy_ch_rn222_desc

character(len=h_mid), dimension(qa_ch), public scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_desc = (/ "Ratio of Rn222 to total mass" /)

Definition at line 39 of file scale_atmos_phy_ch_rn222.F90.

  character(len=H_MID)  , public :: ATMOS_PHY_CH_rn222_DESC(QA_CH) = (/ "Ratio of Rn222 to total mass" /)

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup().

atmos_phy_ch_rn222_unit

character(len=h_short), dimension(qa_ch), public scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_unit = (/ "Bq/kg" /)

Definition at line 40 of file scale_atmos_phy_ch_rn222.F90.

  character(len=H_SHORT), public :: ATMOS_PHY_CH_rn222_UNIT(QA_CH) = (/ "Bq/kg" /)

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup().