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.

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

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_setup().

    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

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.

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency().

    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

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.

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup().

  integer,                public :: atmos_phy_ch_rn222_ntracers = qa_ch

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.

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup().

  character(len=H_SHORT), public :: atmos_phy_ch_rn222_name(qa_ch) = (/ "RN222" /)

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.

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup().

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

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.

Referenced by mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup().

  character(len=H_SHORT), public :: atmos_phy_ch_rn222_unit(qa_ch) = (/ "Bq/kg" /)