mod_atmos_phy_ch_driver Module Reference

module ATMOSPHERE / Physics Chemistry More...

Functions/Subroutines
subroutine, public	atmos_phy_ch_driver_tracer_setup
	Config. More...
subroutine, public	atmos_phy_ch_driver_setup
	Setup. More...
subroutine, public	atmos_phy_ch_driver_finalize
	finalize More...
subroutine, public	atmos_phy_ch_driver_calc_tendency (update_flag)
	Driver. More...
subroutine, public	atmos_phy_ch_driver_ocean_flux (SFLX_QTRC)
	Driver. More...
subroutine, public	atmos_phy_ch_driver_land_flux (SFLX_QTRC)
	Driver. More...
subroutine, public	atmos_phy_ch_driver_urban_flux (SFLX_QTRC)
	Driver. More...

Detailed Description

module ATMOSPHERE / Physics Chemistry

Description

Atmospheric chemistry driver

Author

Team SCALE

NAMELIST

• No namelist group

History Output

name	description	unit	variable
{TRACER_NAME}_t_CH	tendency rho*{TRACER_NAME} in CH; {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	kg/m3/s	RHOQ_t_CH

Function/Subroutine Documentation

atmos_phy_ch_driver_tracer_setup()

subroutine, public mod_atmos_phy_ch_driver::atmos_phy_ch_driver_tracer_setup

Config.

Definition at line 54 of file mod_atmos_phy_ch_driver.F90.

    use mod_atmos_admin, only: &
       atmos_phy_ch_type, &
       atmos_sw_phy_ch
    use scale_tracer, only: &
       tracer_regist
    use scale_atmos_phy_ch_rn222, only: &
       atmos_phy_ch_rn222_ntracers, &
       atmos_phy_ch_rn222_name, &
       atmos_phy_ch_rn222_desc, &
       atmos_phy_ch_rn222_unit
    use mod_atmos_phy_ch_vars, only: &
       qa_ch, &
       qs_ch, &
       qe_ch
    use scale_prc, only: &
       prc_abort
    implicit none
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_CH_driver_tracer_setup",*) 'Setup'
 
    if ( atmos_sw_phy_ch ) then
 
       select case ( atmos_phy_ch_type )
       case ( 'OFF', 'NONE' )
          log_info("ATMOS_PHY_CH_driver_tracer_setup",*) 'this component is never called.'
       case ( 'RN222' )
 
          call tracer_regist( qs_ch,                        & ! [OUT]
                              atmos_phy_ch_rn222_ntracers,  & ! [IN]
                              atmos_phy_ch_rn222_name(:),   & ! [IN]
                              atmos_phy_ch_rn222_desc(:),   & ! [IN]
                              atmos_phy_ch_rn222_unit(:)    ) ! [IN]
          qa_ch = atmos_phy_ch_rn222_ntracers
 
       case default
          log_error("ATMOS_PHY_CH_driver_tracer_setup",*) 'invalid chemistry type(', atmos_phy_ch_type, '). CHECK!'
          call prc_abort
       end select
 
       qe_ch = qs_ch + qa_ch - 1
 
    else
       qa_ch = 0
       qs_ch = -1
       qe_ch = -2
    endif
 
    return

References scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_desc, scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_name, scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_ntracers, scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_unit, mod_atmos_admin::atmos_phy_ch_type, mod_atmos_admin::atmos_sw_phy_ch, scale_prc::prc_abort(), mod_atmos_phy_ch_vars::qa_ch, mod_atmos_phy_ch_vars::qe_ch, mod_atmos_phy_ch_vars::qs_ch, and scale_tracer::tracer_regist().

Referenced by mod_atmos_driver::atmos_driver_tracer_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ch_driver_setup()

subroutine, public mod_atmos_phy_ch_driver::atmos_phy_ch_driver_setup

Setup.

Definition at line 109 of file mod_atmos_phy_ch_driver.F90.

    use scale_atmos_grid_cartesc_real, only: &
       real_lon => atmos_grid_cartesc_real_lon, &
       real_lat => atmos_grid_cartesc_real_lat
    use scale_atmos_phy_ch_rn222, only: &
       atmos_phy_ch_rn222_setup
    use scale_atmos_sfc_ch_rn222, only: &
       atmos_sfc_ch_rn222_setup
    use mod_atmos_admin, only: &
       atmos_phy_ch_type, &
       atmos_sw_phy_ch
    implicit none
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_CH_driver_setup",*) 'Setup'
 
    if ( atmos_sw_phy_ch ) then
 
       select case ( atmos_phy_ch_type )
       case ( 'RN222' )
          call atmos_phy_ch_rn222_setup
          call atmos_sfc_ch_rn222_setup( ia, ja,            & ! [IN]
                                         real_lon, real_lat ) ! [IN]
       end select
 
    else
       log_info("ATMOS_PHY_CH_driver_setup",*) 'this component is never called.'
    endif
 
    return

References scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_lat, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_lon, scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_setup(), mod_atmos_admin::atmos_phy_ch_type, scale_atmos_sfc_ch_rn222::atmos_sfc_ch_rn222_setup(), mod_atmos_admin::atmos_sw_phy_ch, scale_atmos_grid_cartesc_index::ia, and scale_atmos_grid_cartesc_index::ja.

Referenced by mod_atmos_driver::atmos_driver_setup().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ch_driver_finalize()

subroutine, public mod_atmos_phy_ch_driver::atmos_phy_ch_driver_finalize

finalize

Definition at line 144 of file mod_atmos_phy_ch_driver.F90.

    use scale_atmos_sfc_ch_rn222, only: &
       atmos_sfc_ch_rn222_finalize
    use mod_atmos_admin, only: &
       atmos_phy_ch_type, &
       atmos_sw_phy_ch
    implicit none
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_CH_driver_finalize",*) 'Finalize'
 
    if ( atmos_sw_phy_ch ) then
 
       select case ( atmos_phy_ch_type )
       case ( 'RN222' )
          call atmos_sfc_ch_rn222_finalize
       end select
 
    end if
 
    return

References mod_atmos_admin::atmos_phy_ch_type, scale_atmos_sfc_ch_rn222::atmos_sfc_ch_rn222_finalize(), and mod_atmos_admin::atmos_sw_phy_ch.

Referenced by mod_atmos_driver::atmos_driver_finalize().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ch_driver_calc_tendency()

subroutine, public mod_atmos_phy_ch_driver::atmos_phy_ch_driver_calc_tendency

(

logical, intent(in)

update_flag

)

Driver.

Definition at line 170 of file mod_atmos_phy_ch_driver.F90.

    use scale_tracer, only: &
       tracer_name
    use scale_time, only: &
       dt_ch => time_dtsec_atmos_phy_ch
    use scale_statistics, only: &
       statistics_checktotal, &
       statistics_total
    use scale_atmos_grid_cartesc_real, only: &
       atmos_grid_cartesc_real_vol, &
       atmos_grid_cartesc_real_totvol
    use scale_file_history, only: &
       file_history_in
    use mod_atmos_vars, only: &
       dens   => dens_av, &
       qtrc   => qtrc_av, &
       rhoq_t => rhoq_tp
    use mod_atmos_phy_sf_vars, only: &
       sflx_qtrc => atmos_phy_sf_sflx_qtrc
    use mod_atmos_phy_ch_vars, only: &
       rhoq_t_ch => atmos_phy_ch_rhoq_t, &
       qa_ch,        &
       qs_ch,        &
       qe_ch
!       O3        => ATMOS_PHY_CH_O3
    use mod_atmos_admin, only: &
       atmos_phy_ch_type
    use scale_atmos_phy_ch_rn222, only: &
       atmos_phy_ch_rn222_tendency
    implicit none
 
    logical, intent(in) :: update_flag
 
    integer  :: k, i, j, iq
    !---------------------------------------------------------------------------
 
    if ( update_flag ) then
 
!OCL XFILL
       !$acc kernels
       rhoq_t_ch(:,:,:,:) = 0.0_rp
       !$acc end kernels
 
       select case( atmos_phy_ch_type )
       case( 'RN222' )
          !$acc update host(DENS,QTRC(:,:,:,QS_CH:QE_CH),RHOQ_t_CH(:,:,:,QS_CH:QE_CH))
          call atmos_phy_ch_rn222_tendency( ka, ks, ke,                   & ! [IN]
                                            ia, is, ie,                   & ! [IN]
                                            ja, js, je,                   & ! [IN]
                                            qa_ch,                        & ! [IN]
                                            dens(:,:,:),             & ! [IN]
                                            qtrc(:,:,:,qs_ch:qe_ch), & ! [IN]
                                            rhoq_t_ch(:,:,:,qs_ch:qe_ch)  ) ! [INOUT]
          !$acc update device(RHOQ_t_CH(:,:,:,QS_CH:QE_CH))
       end select
 
       do iq = qs_ch, qe_ch
          call file_history_in( rhoq_t_ch(:,:,:,iq), trim(tracer_name(iq))//'_t_CH', &
                                'tendency rho*'//trim(tracer_name(iq))//' in CH',    &
                                'kg/m3/s', fill_halo=.true.                          )
       enddo
    endif
 
    !$omp parallel do private(i,j,k) collapse(3) OMP_SCHEDULE_
    !$acc kernels
    do iq = qs_ch, qe_ch
       do j  = js, je
       do i  = is, ie
       do k  = ks, ke
          rhoq_t(k,i,j,iq) = rhoq_t(k,i,j,iq) + rhoq_t_ch(k,i,j,iq)
       enddo
       enddo
       enddo
    enddo
    !$acc end kernels
 
    if ( statistics_checktotal ) then
       do iq = qs_ch, qe_ch
          call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                 rhoq_t_ch(:,:,:,iq), trim(tracer_name(iq))//'_t_CH', &
                                 atmos_grid_cartesc_real_vol(:,:,:),                  &
                                 atmos_grid_cartesc_real_totvol                       )
       enddo
    endif
 
    return

References scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvol, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_vol, mod_atmos_phy_ch_vars::atmos_phy_ch_rhoq_t, scale_atmos_phy_ch_rn222::atmos_phy_ch_rn222_tendency(), mod_atmos_admin::atmos_phy_ch_type, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc, mod_atmos_vars::dens, mod_atmos_vars::dens_av, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mod_atmos_phy_ch_vars::qa_ch, mod_atmos_phy_ch_vars::qe_ch, mod_atmos_phy_ch_vars::qs_ch, mod_atmos_vars::qtrc, mod_atmos_vars::qtrc_av, mod_atmos_vars::rhoq_tp, scale_statistics::statistics_checktotal, scale_time::time_dtsec_atmos_phy_ch, and scale_tracer::tracer_name.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ch_driver_ocean_flux()

subroutine, public mod_atmos_phy_ch_driver::atmos_phy_ch_driver_ocean_flux

(

real(rp), dimension(oia,oja,qa), intent(inout)

SFLX_QTRC

)

Driver.

Definition at line 261 of file mod_atmos_phy_ch_driver.F90.

    use scale_ocean_grid_cartesc_index
    use scale_atmos_sfc_ch_rn222, only: &
       atmos_sfc_ch_rn222_ocean_flux
    use mod_atmos_admin, only: &
       atmos_phy_ch_type
    use mod_atmos_phy_ch_vars, only: &
       qa_ch, &
       qs_ch, &
       qe_ch
    implicit none
 
    real(RP), intent(inout) :: SFLX_QTRC(OIA,OJA,QA)
 
    !---------------------------------------------------------------------------
 
    select case( atmos_phy_ch_type )
    case( 'RN222' )
       !$acc update host(SFLX_QTRC(:,:,QS_CH:QE_CH))
       call atmos_sfc_ch_rn222_ocean_flux( oia, ois, oie,             & ! [IN]
                                           oja, ojs, oje,             & ! [IN]
                                           qa_ch,                     & ! [IN]
                                           sflx_qtrc(:,:,qs_ch:qe_ch) ) ! [INOUT]
       !$acc update device(SFLX_QTRC(:,:,QS_CH:QE_CH))
    end select
 
    return

References mod_atmos_admin::atmos_phy_ch_type, scale_atmos_sfc_ch_rn222::atmos_sfc_ch_rn222_ocean_flux(), scale_ocean_grid_cartesc_index::oia, scale_ocean_grid_cartesc_index::oie, scale_ocean_grid_cartesc_index::ois, scale_ocean_grid_cartesc_index::oja, scale_ocean_grid_cartesc_index::oje, scale_ocean_grid_cartesc_index::ojs, mod_atmos_phy_ch_vars::qa_ch, mod_atmos_phy_ch_vars::qe_ch, and mod_atmos_phy_ch_vars::qs_ch.

Referenced by mod_ocean_driver::ocean_driver_calc_tendency().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ch_driver_land_flux()

subroutine, public mod_atmos_phy_ch_driver::atmos_phy_ch_driver_land_flux

(

real(rp), dimension(lia,lja,qa), intent(inout)

SFLX_QTRC

)

Driver.

Definition at line 293 of file mod_atmos_phy_ch_driver.F90.

    use scale_land_grid_cartesc_index
    use scale_time, only: &
       time_nowdate
    use scale_atmos_sfc_ch_rn222, only: &
       atmos_sfc_ch_rn222_land_flux
    use mod_atmos_admin, only: &
       atmos_phy_ch_type
    use mod_atmos_phy_ch_vars, only: &
       qa_ch, &
       qs_ch, &
       qe_ch
    implicit none
 
    real(RP), intent(inout) :: SFLX_QTRC(LIA,LJA,QA)
 
    !---------------------------------------------------------------------------
 
    select case( atmos_phy_ch_type )
    case( 'RN222' )
       !$acc update host(SFLX_QTRC(:,:,QS_CH:QE_CH))
       call atmos_sfc_ch_rn222_land_flux( lia, lis, lie,             & ! [IN]
                                          lja, ljs, lje,             & ! [IN]
                                          qa_ch,                     & ! [IN]
                                          time_nowdate(:),           & ! [IN]
                                          sflx_qtrc(:,:,qs_ch:qe_ch) ) ! [INOUT]
       !$acc update device(SFLX_QTRC(:,:,QS_CH:QE_CH))
    end select
 
    return

References mod_atmos_admin::atmos_phy_ch_type, scale_atmos_sfc_ch_rn222::atmos_sfc_ch_rn222_land_flux(), scale_land_grid_cartesc_index::lia, scale_land_grid_cartesc_index::lie, scale_land_grid_cartesc_index::lis, scale_land_grid_cartesc_index::lja, scale_land_grid_cartesc_index::lje, scale_land_grid_cartesc_index::ljs, mod_atmos_phy_ch_vars::qa_ch, mod_atmos_phy_ch_vars::qe_ch, mod_atmos_phy_ch_vars::qs_ch, and scale_time::time_nowdate.

Referenced by mod_land_driver::land_driver_calc_tendency().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ch_driver_urban_flux()

subroutine, public mod_atmos_phy_ch_driver::atmos_phy_ch_driver_urban_flux

(

real(rp), dimension(uia,uja,qa), intent(inout)

SFLX_QTRC

)

Driver.

Definition at line 328 of file mod_atmos_phy_ch_driver.F90.

    use scale_urban_grid_cartesc_index
    use scale_time, only: &
       time_nowdate
    use scale_atmos_sfc_ch_rn222, only: &
       atmos_sfc_ch_rn222_land_flux
    use mod_atmos_admin, only: &
       atmos_phy_ch_type
    use mod_atmos_phy_ch_vars, only: &
       qa_ch, &
       qs_ch, &
       qe_ch
    implicit none
 
    real(RP), intent(inout) :: SFLX_QTRC(UIA,UJA,QA)
 
    !---------------------------------------------------------------------------
 
    select case( atmos_phy_ch_type )
    case( 'RN222' )
       !$acc update host(SFLX_QTRC(:,:,QS_CH:QE_CH))
       call atmos_sfc_ch_rn222_land_flux( uia, uis, uie,             & ! [IN]
                                          uja, ujs, uje,             & ! [IN]
                                          qa_ch,                     & ! [IN]
                                          time_nowdate(:),           & ! [IN]
                                          sflx_qtrc(:,:,qs_ch:qe_ch) ) ! [INOUT]
       !$acc update device(SFLX_QTRC(:,:,QS_CH:QE_CH))
    end select
 
    return

References mod_atmos_admin::atmos_phy_ch_type, scale_atmos_sfc_ch_rn222::atmos_sfc_ch_rn222_land_flux(), mod_atmos_phy_ch_vars::qa_ch, mod_atmos_phy_ch_vars::qe_ch, mod_atmos_phy_ch_vars::qs_ch, scale_time::time_nowdate, scale_urban_grid_cartesc_index::uia, scale_urban_grid_cartesc_index::uie, scale_urban_grid_cartesc_index::uis, scale_urban_grid_cartesc_index::uja, scale_urban_grid_cartesc_index::uje, and scale_urban_grid_cartesc_index::ujs.

Referenced by mod_urban_driver::urban_driver_calc_tendency().

Here is the call graph for this function:

Here is the caller graph for this function: