mod_atmos_phy_ae_driver Module Reference

module ATMOSPHERE / Physics Aerosol Microphysics More...

Functions/Subroutines
subroutine, public	atmos_phy_ae_driver_tracer_setup
	Setup. More...
subroutine, public	atmos_phy_ae_driver_setup
	Setup. More...
subroutine, public	atmos_phy_ae_driver_adjustment
	adjustment More...
subroutine, public	atmos_phy_ae_driver_calc_tendency (update_flag)
	Driver. More...

Detailed Description

module ATMOSPHERE / Physics Aerosol Microphysics

Description

Aerosol Microphysics driver

Author

Team SCALE

NAMELIST

• No namelist group

History Output

name	description	unit	variable
CCN	cloud condensation nucrei	num/cc	CCN*1.E-6_RP
CN	condensation nucrei	num/cc	CN *1.E-6_RP

Function/Subroutine Documentation

atmos_phy_ae_driver_tracer_setup()

subroutine, public mod_atmos_phy_ae_driver::atmos_phy_ae_driver_tracer_setup

Setup.

Definition at line 50 of file mod_atmos_phy_ae_driver.F90.

    use mod_atmos_admin, only: &
       atmos_phy_ae_type, &
       atmos_sw_phy_ae
    use scale_tracer, only: &
       tracer_regist
    use scale_atmos_phy_ae_kajino13, only: &
       atmos_phy_ae_kajino13_tracer_setup, &
       atmos_phy_ae_kajino13_name, &
       atmos_phy_ae_kajino13_desc, &
       atmos_phy_ae_kajino13_unit
    use scale_prc, only: &
       prc_abort
    use mod_atmos_phy_ae_vars, only: &
       qa_ae, &
       qs_ae, &
       qe_ae
    implicit none
 
    log_newline
    log_info("ATMOS_PHY_AE_driver_tracer_setup",*) 'Setup'
 
    if ( atmos_sw_phy_ae ) then
       select case ( atmos_phy_ae_type )
       case ( 'OFF', 'NONE' )
          log_info("ATMOS_PHY_AE_driver_tracer_setup",*) 'this component is never called.'
       case ( 'KAJINO13' )
          call atmos_phy_ae_kajino13_tracer_setup( qa_ae ) ! [OUT]
 
          call tracer_regist( qs_ae,                         & ! [OUT]
                              qa_ae,                         & ! [IN]
                              atmos_phy_ae_kajino13_name(:), & ! [IN]
                              atmos_phy_ae_kajino13_desc(:), & ! [IN]
                              atmos_phy_ae_kajino13_unit(:)  ) ! [IN]
       case ( 'OFFLINE' )
          log_info("ATMOS_PHY_AE_driver_tracer_setup",*) 'offline aerosol module has no tracers'
       case default
          log_error("ATMOS_PHY_AE_driver_tracer_setup",*) 'invalid aerosol type(', atmos_phy_ae_type, '). CHECK!'
          call prc_abort
       end select
 
       qe_ae = qs_ae + qa_ae - 1
 
    else
       qa_ae = 0
       qs_ae = -1
       qe_ae = -2
    end if
 
    return

References scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_desc, scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_name, scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_tracer_setup(), scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_unit, mod_atmos_admin::atmos_phy_ae_type, mod_atmos_admin::atmos_sw_phy_ae, scale_prc::prc_abort(), mod_atmos_phy_ae_vars::qa_ae, mod_atmos_phy_ae_vars::qe_ae, mod_atmos_phy_ae_vars::qs_ae, 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_ae_driver_setup()

subroutine, public mod_atmos_phy_ae_driver::atmos_phy_ae_driver_setup

Setup.

Definition at line 104 of file mod_atmos_phy_ae_driver.F90.

    use mod_atmos_admin, only: &
       atmos_phy_ae_type, &
       atmos_sw_phy_ae
    use scale_atmos_phy_ae_kajino13, only: &
        atmos_phy_ae_kajino13_setup
    use scale_atmos_phy_ae_offline, only: &
        atmos_phy_ae_offline_setup
    use scale_prc, only: &
       prc_abort
    implicit none
    !---------------------------------------------------------------------------
 
    log_newline
    log_info("ATMOS_PHY_AE_driver_setup",*) 'Setup'
 
    if ( atmos_sw_phy_ae ) then
 
       select case ( atmos_phy_ae_type )
       case ( 'KAJINO13' )
          call atmos_phy_ae_kajino13_setup
       case ( 'OFFLINE' )
          call atmos_phy_ae_offline_setup
       case default
          log_error("ATMOS_PHY_AE_driver_setup",*) 'invalid aerosol type(', atmos_phy_ae_type, '). CHECK!'
          call prc_abort
       end select
 
    endif
 
    return

References scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_setup(), scale_atmos_phy_ae_offline::atmos_phy_ae_offline_setup(), mod_atmos_admin::atmos_phy_ae_type, mod_atmos_admin::atmos_sw_phy_ae, and scale_prc::prc_abort().

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_ae_driver_adjustment()

subroutine, public mod_atmos_phy_ae_driver::atmos_phy_ae_driver_adjustment

adjustment

Definition at line 139 of file mod_atmos_phy_ae_driver.F90.

    use mod_atmos_vars, only: &
       qtrc
    use mod_atmos_phy_ae_vars, only: &
       qa_ae, &
       qs_ae, &
       qe_ae
    use mod_atmos_admin, only: &
       atmos_phy_ae_type, &
       atmos_sw_phy_ae
    use scale_atmos_phy_ae_kajino13, only: &
        atmos_phy_ae_kajino13_negative_fixer
    implicit none
 
    if ( atmos_sw_phy_ae ) then
 
       select case ( atmos_phy_ae_type )
       case ( 'KAJINO13' )
          call atmos_phy_ae_kajino13_negative_fixer( ka, ks, ke, ia, is, ie, ja, js, je, qa_ae, &
                                                     qtrc(:,:,:,qs_ae:qe_ae) ) ! [INOUT]
       end select
 
    end if
 
    return

References scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_negative_fixer(), mod_atmos_admin::atmos_phy_ae_type, mod_atmos_admin::atmos_sw_phy_ae, 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_ae_vars::qa_ae, mod_atmos_phy_ae_vars::qe_ae, mod_atmos_phy_ae_vars::qs_ae, and mod_atmos_vars::qtrc.

Referenced by mod_atmos_driver::atmos_driver_update().

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_ae_driver_calc_tendency()

subroutine, public mod_atmos_phy_ae_driver::atmos_phy_ae_driver_calc_tendency

(

logical, intent(in)

update_flag

)

Driver.

Definition at line 168 of file mod_atmos_phy_ae_driver.F90.

    use scale_tracer, only: &
       tracer_name
    use scale_prc, only: &
       prc_abort
    use scale_time, only: &
       dt_ae => time_dtsec_atmos_phy_ae, &
       time_nowdaysec
    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, &
       qdry, &
       pres, &
       temp, &
       qv,   &
       rhoq_t => rhoq_tp
    use mod_atmos_phy_ae_vars, only: &
       qa_ae, &
       qs_ae, &
       qe_ae, &
       rhoq_t_ae => atmos_phy_ae_rhoq_t, &
       ccn       => atmos_phy_ae_ccn,   &
       ccn_t     => atmos_phy_ae_ccn_t, &
       ae_emit   => atmos_phy_ae_emit
    use mod_atmos_phy_mp_vars, only: &
       evaporate => atmos_phy_mp_evaporate
    use mod_atmos_admin, only: &
       atmos_phy_ae_type
    use scale_atmos_phy_ae_kajino13, only: &
       atmos_phy_ae_kajino13_tendency
    use scale_atmos_phy_ae_offline, only: &
       atmos_phy_ae_offline_tendency
    implicit none
 
    logical, intent(in) :: update_flag
 
    real(RP) :: CN(KA,IA,JA)
    real(RP) :: NREG(KA,IA,JA)
 
    integer  :: k, i, j, iq
    !---------------------------------------------------------------------------
 
    if ( update_flag ) then
 
!OCL XFILL
       ccn(:,:,:)   = 0.0_rp ! reset
!OCL XFILL
       rhoq_t_ae(:,:,:,:) = 0.0_rp ! reset
 
       do j  = js, je
       do i  = is, ie
       do k  = ks, ke
          nreg(k,i,j) = evaporate(k,i,j) * dt_ae
       enddo
       enddo
       enddo
 
       select case ( atmos_phy_ae_type )
       case ( 'KAJINO13' )
          call atmos_phy_ae_kajino13_tendency( ka, ks, ke, ia, is, ie, ja, js, je, qa_ae, &
                                               temp(:,:,:),             & ! [IN]
                                               pres(:,:,:),             & ! [IN]
                                               qdry(:,:,:),             & ! [IN]
                                               nreg(:,:,:),             & ! [IN]
                                               dens(:,:,:),             & ! [IN]
                                               qv(:,:,:),             & ! [IN]
                                               qtrc(:,:,:,qs_ae:qe_ae), & ! [IN]
                                               ae_emit(:,:,:,qs_ae:qe_ae), & ! [IN]
                                               dt_ae,                        & ! [IN]
                                               rhoq_t_ae(:,:,:,qs_ae:qe_ae), & ! [OUT]
                                               cn(:,:,:),             & ! [OUT]
                                               ccn(:,:,:)              ) ! [OUT]
       case ( 'OFFLINE' )
          call atmos_phy_ae_offline_tendency ( ka, ks, ke, ia, is, ie, ja, js, je, &
                                               time_nowdaysec,               & ! [IN]
                                               ccn(:,:,:)              ) ! [OUT]
          cn(:,:,:) = 0.0_rp ! not supported
       end select
 
       ccn_t(:,:,:) = ccn(:,:,:) / dt_ae
 
       call file_history_in( cn(:,:,:)*1.e-6_rp, 'CN',  'condensation nucrei',       'num/cc' )
       call file_history_in( ccn(:,:,:)*1.e-6_rp, 'CCN', 'cloud condensation nucrei', 'num/cc' )
 
    endif
 
    do iq = qs_ae, qe_ae
       !$omp parallel do private(i,j,k) OMP_SCHEDULE_
       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_ae(k,i,j,iq)
       enddo
       enddo
       enddo
    enddo
 
    if ( statistics_checktotal ) then
       do iq = qs_ae, qe_ae
          call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                 rhoq_t_ae(:,:,:,iq), trim(tracer_name(iq))//'_t_AE', &
                                 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_ae_vars::atmos_phy_ae_ccn, mod_atmos_phy_ae_vars::atmos_phy_ae_ccn_t, mod_atmos_phy_ae_vars::atmos_phy_ae_emit, scale_atmos_phy_ae_kajino13::atmos_phy_ae_kajino13_tendency(), scale_atmos_phy_ae_offline::atmos_phy_ae_offline_tendency(), mod_atmos_phy_ae_vars::atmos_phy_ae_rhoq_t, mod_atmos_admin::atmos_phy_ae_type, mod_atmos_phy_mp_vars::atmos_phy_mp_evaporate, 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, scale_prc::prc_abort(), mod_atmos_vars::pres, mod_atmos_phy_ae_vars::qa_ae, mod_atmos_vars::qdry, mod_atmos_phy_ae_vars::qe_ae, mod_atmos_phy_ae_vars::qs_ae, mod_atmos_vars::qtrc, mod_atmos_vars::qtrc_av, mod_atmos_vars::qv, mod_atmos_vars::rhoq_tp, scale_statistics::statistics_checktotal, mod_atmos_vars::temp, scale_time::time_dtsec_atmos_phy_ae, scale_time::time_nowdaysec, 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: