mod_atmos_phy_bl_driver Module Reference

module atmosphere / physics / PBL More...

Functions/Subroutines
subroutine, public	atmos_phy_bl_driver_tracer_setup
	Config. More...
subroutine, public	atmos_phy_bl_driver_setup
	Setup. More...
subroutine, public	atmos_phy_bl_driver_calc_tendency (update_flag)
	calculate tendency More...

Detailed Description

module atmosphere / physics / PBL

Description

Planetary boundary layer turbulence

Author

Team SCALE

NAMELIST

• No namelist group

History Output

name	description	unit	variable
Kh_BL	eddy diffusion	m2/s	Kh
Nu_BL	eddy viscosity	m2/s	Nu
RHOT_t_BL	RHOT tendency (BL)	K.kg/m3/s	RHOT_t_BL
RHOU_t_BL	MOMX tendency (BL)	kg/m2/s2	RHOU_t_BL
RHOV_t_BL	MOMY tendency (BL)	kg/m2/s2	RHOV_t_BL
{TRACER_NAME}_t_BL	RHO*{TRACER_NAME} tendency (BL); {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	kg/m3/s	RHOQ_t_BL

Function/Subroutine Documentation

atmos_phy_bl_driver_tracer_setup()

subroutine, public mod_atmos_phy_bl_driver::atmos_phy_bl_driver_tracer_setup

(

)

Config.

Definition at line 50 of file mod_atmos_phy_bl_driver.F90.

References mod_atmos_admin::atmos_phy_bl_type, mod_atmos_admin::atmos_sw_phy_bl, scale_prc::prc_abort(), mod_atmos_phy_bl_vars::qe, mod_atmos_phy_bl_vars::qs, and scale_tracer::tracer_regist().

Referenced by mod_atmos_driver::atmos_driver_tracer_setup().

    use scale_prc, only: &
       prc_abort
    use scale_tracer, only: &
       tracer_regist
    use scale_atmos_phy_bl_mynn, only: &
       atmos_phy_bl_mynn_tracer_setup, &
       atmos_phy_bl_mynn_ntracer, &
       atmos_phy_bl_mynn_name, &
       atmos_phy_bl_mynn_desc, &
       atmos_phy_bl_mynn_units
    use mod_atmos_admin, only: &
       atmos_phy_bl_type, &
       atmos_sw_phy_bl
    use mod_atmos_phy_bl_vars, only: &
       qs, qe
    implicit none
    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_PHY_BL_driver_tracer_setup",*) 'Setup'

    if ( atmos_sw_phy_bl ) then
       select case ( atmos_phy_bl_type )
       case ( 'MYNN' )
          call atmos_phy_bl_mynn_tracer_setup
          call tracer_regist( &
               qs, &
               atmos_phy_bl_mynn_ntracer, &
               atmos_phy_bl_mynn_name,    &
               atmos_phy_bl_mynn_desc,    &
               atmos_phy_bl_mynn_units    )
          qe = qs + atmos_phy_bl_mynn_ntracer - 1
       case default
          log_error("ATMOS_PHY_BL_driver_tracer_setup",*) 'ATMOS_PHY_BL_TYPE is invalid: ', trim(atmos_phy_bl_type)
          call prc_abort
       end select
    end if

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_bl_driver_setup()

subroutine, public mod_atmos_phy_bl_driver::atmos_phy_bl_driver_setup

(

)

Setup.

Definition at line 94 of file mod_atmos_phy_bl_driver.F90.

References scale_atmos_grid_cartesc::atmos_grid_cartesc_cdx, scale_atmos_grid_cartesc::atmos_grid_cartesc_cdy, scale_atmos_grid_cartesc::atmos_grid_cartesc_cdz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, mod_atmos_admin::atmos_phy_bl_type, mod_atmos_admin::atmos_sw_phy_bl, 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, and scale_atmos_grid_cartesc_index::ks.

Referenced by mod_atmos_driver::atmos_driver_setup().

    use scale_atmos_grid_cartesc, only: &
       cdz => atmos_grid_cartesc_cdz, &
       cdx => atmos_grid_cartesc_cdx, &
       cdy => atmos_grid_cartesc_cdy
    use scale_atmos_grid_cartesc_real, only: &
       real_cz  => atmos_grid_cartesc_real_cz
    use scale_atmos_phy_bl_mynn, only: &
       atmos_phy_bl_mynn_setup
    use mod_atmos_admin, only: &
       atmos_phy_bl_type, &
       atmos_sw_phy_bl
    implicit none
    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_PHY_BL_driver_setup",*) 'Setup'

    if ( atmos_sw_phy_bl ) then
       select case ( atmos_phy_bl_type )
       case ( 'MYNN' )
          call atmos_phy_bl_mynn_setup( &
               ka, ks, ke, ia, is, ie, ja, js, je, &
               real_cz ) ! (in)
       end select
    else
       log_info("ATMOS_PHY_BL_driver_setup",*) 'this component is never called.'
    endif

    return

Here is the caller graph for this function:

atmos_phy_bl_driver_calc_tendency()

subroutine, public mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency

(

logical, intent(in)

update_flag

)

calculate tendency

Definition at line 128 of file mod_atmos_phy_bl_driver.F90.

References scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_cz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_fz, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_totvol, scale_atmos_grid_cartesc_real::atmos_grid_cartesc_real_vol, mod_atmos_phy_bl_vars::atmos_phy_bl_rhoq_t, mod_atmos_phy_bl_vars::atmos_phy_bl_rhot_t, mod_atmos_phy_bl_vars::atmos_phy_bl_rhou_t, mod_atmos_phy_bl_vars::atmos_phy_bl_rhov_t, mod_atmos_admin::atmos_phy_bl_type, mod_atmos_phy_sf_vars::atmos_phy_sf_l_mo, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mu, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_mv, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qv, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_sh, mod_atmos_admin::atmos_sw_phy_bl, mod_atmos_vars::dens, mod_atmos_vars::dens_av, mod_atmos_vars::exner, scale_atmos_grid_cartesc_index::ia, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::ieb, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::isb, scale_atmos_grid_cartesc_index::ja, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::jeb, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::jsb, scale_atmos_grid_cartesc_index::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, mod_atmos_vars::pott, mod_atmos_vars::pres, scale_tracer::qa, mod_atmos_vars::qc, mod_atmos_vars::qdry, mod_atmos_phy_bl_vars::qe, mod_atmos_vars::qi, mod_atmos_phy_bl_vars::qs, mod_atmos_vars::qtrc, mod_atmos_vars::qtrc_av, mod_atmos_vars::qv, mod_atmos_vars::rhoq_tp, mod_atmos_vars::rhot_tp, mod_atmos_vars::rhou_tp, mod_atmos_vars::rhov_tp, scale_statistics::statistics_checktotal, scale_time::time_dtsec_atmos_phy_bl, scale_tracer::tracer_advc, scale_tracer::tracer_name, mod_atmos_vars::u, and mod_atmos_vars::v.

Referenced by mod_atmos_driver::atmos_driver_calc_tendency(), and mod_atmos_driver::atmos_driver_calc_tendency_from_sflux().

    use scale_statistics, only: &
       statistics_checktotal, &
       statistics_total
    use scale_file_history, only: &
       file_history_in
    use scale_time, only: &
       dt_bl => time_dtsec_atmos_phy_bl
    use scale_atmos_phy_bl_mynn, only: &
       atmos_phy_bl_mynn_tendency, &
       atmos_phy_bl_mynn_tendency_tracer
    use scale_atmos_grid_cartesc_real, only: &
       cz => atmos_grid_cartesc_real_cz, &
       fz => atmos_grid_cartesc_real_fz, &
       atmos_grid_cartesc_real_vol, &
       atmos_grid_cartesc_real_totvol
    use mod_atmos_admin, only: &
       atmos_phy_bl_type, &
       atmos_sw_phy_bl
    use mod_atmos_vars, only: &
       dens => dens_av, &
       qtrc => qtrc_av, &
       u,     &
       v,     &
       pott,  &
       pres,  &
       exner, &
       qdry,  &
       qv,    &
       qc,    &
       qi,    &
       rhou_t => rhou_tp, &
       rhov_t => rhov_tp, &
       rhot_t => rhot_tp, &
       rhoq_t => rhoq_tp, &
       atmos_vars_get_diagnostic
    use mod_atmos_phy_bl_vars, only: &
       qs, qe, &
       rhou_t_bl => atmos_phy_bl_rhou_t, &
       rhov_t_bl => atmos_phy_bl_rhov_t, &
       rhot_t_bl => atmos_phy_bl_rhot_t, &
       rhoq_t_bl => atmos_phy_bl_rhoq_t
    use mod_atmos_phy_sf_vars, only: &
       sflx_mu => atmos_phy_sf_sflx_mu, &
       sflx_mv => atmos_phy_sf_sflx_mv, &
       sflx_sh => atmos_phy_sf_sflx_sh, &
       sflx_q  => atmos_phy_sf_sflx_qtrc, &
       sflx_qv => atmos_phy_sf_sflx_qv, &
       l_mo    => atmos_phy_sf_l_mo
    implicit none

    logical, intent(in) :: update_flag

    real(RP) :: nu(ka,ia,ja)
    real(RP) :: kh(ka,ia,ja)
    real(RP) :: qw(ka,ia,ja)

    real(RP) :: n2  (ka,ia,ja)
    real(RP) :: potl(ka,ia,ja)
    real(RP) :: potv(ka,ia,ja)

    integer  :: k, i, j, iq
    !---------------------------------------------------------------------------

    if ( update_flag ) then

       rhoq_t_bl(:,:,:,:) = 0.0_rp

       select case ( atmos_phy_bl_type )
       case ( 'MYNN' )
          call atmos_vars_get_diagnostic( "N2",   n2   )
          call atmos_vars_get_diagnostic( "POTL", potl )
          call atmos_vars_get_diagnostic( "POTV", potv )
          do j = jsb, jeb
          do i = isb, ieb
          do k = ks, ke
             qw(k,i,j) = qv(k,i,j) + qc(k,i,j) + qi(k,i,j)
          end do
          end do
          end do
          call atmos_phy_bl_mynn_tendency( &
               ka, ks, ke, ia, is, ie, ja, js, je, &
               dens(:,:,:), u(:,:,:), v(:,:,:),                        & ! (in)
               pott(:,:,:), qtrc(:,:,:,qs:qe),                         & ! (in)
               pres(:,:,:), exner(:,:,:), n2(:,:,:),                   & ! (in)
               qdry(:,:,:), qv(:,:,:), qw(:,:,:),                      & ! (in)
               potl(:,:,:), potv(:,:,:),                               & ! (in)
               sflx_mu(:,:), sflx_mv(:,:), sflx_sh(:,:), sflx_qv(:,:), & ! (in)
               l_mo(:,:),                                              & ! (in)
               cz(:,:,:), fz(:,:,:), dt_bl,                            & ! (in)
               rhou_t_bl(:,:,:), rhov_t_bl(:,:,:),                     & ! (out)
               rhot_t_bl(:,:,:), rhoq_t_bl(:,:,:,qs:qe),               & ! (out)
               nu(:,:,:), kh(:,:,:)                                    ) ! (out)
          do iq = 1, qa
             if ( ( .not. tracer_advc(iq) ) .or. (iq>=qs .and. iq<=qe) ) cycle
             call atmos_phy_bl_mynn_tendency_tracer( &
                  ka, ks, ke, ia, is, ie, ja, js, je, &
                  dens(:,:,:), qtrc(:,:,:,iq), & ! (in)
                  sflx_q(:,:,iq), kh(:,:,:),   & ! (in)
                  cz(:,:,:), fz(:,:,:),        & ! (in)
                  dt_bl, tracer_name(iq),      & ! (in)
                  rhoq_t_bl(:,:,:,iq)          ) ! (out)
          end do
       end select

       call file_history_in( nu(:,:,:),        'Nu_BL',     'eddy viscosity',     'm2/s',      fill_halo=.true. )
       call file_history_in( kh(:,:,:),        'Kh_BL',     'eddy diffusion',     'm2/s',      fill_halo=.true. )

       call file_history_in( rhou_t_bl(:,:,:), 'RHOU_t_BL', 'MOMX tendency (BL)', 'kg/m2/s2',  fill_halo=.true. )
       call file_history_in( rhov_t_bl(:,:,:), 'RHOV_t_BL', 'MOMY tendency (BL)', 'kg/m2/s2',  fill_halo=.true. )
       call file_history_in( rhot_t_bl(:,:,:), 'RHOT_t_BL', 'RHOT tendency (BL)', 'K.kg/m3/s', fill_halo=.true. )

       do iq = 1, qa
          if ( .not. tracer_advc(iq) ) cycle
          call file_history_in( rhoq_t_bl(:,:,:,iq), trim(tracer_name(iq))//'_t_BL',                      &
                        'RHO*'//trim(tracer_name(iq))//' tendency (BL)', 'kg/m3/s', fill_halo=.true. )
       enddo

       if ( statistics_checktotal ) then
          call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                 rhou_t_bl(:,:,:), 'RHOU_t_BL',      &
                                 atmos_grid_cartesc_real_vol(:,:,:), &
                                 atmos_grid_cartesc_real_totvol      )
          call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                 rhov_t_bl(:,:,:), 'RHOV_t_BL',      &
                                 atmos_grid_cartesc_real_vol(:,:,:), &
                                 atmos_grid_cartesc_real_totvol      )
          call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                 rhot_t_bl(:,:,:), 'RHOT_t_BL',      &
                                 atmos_grid_cartesc_real_vol(:,:,:), &
                                 atmos_grid_cartesc_real_totvol      )
          call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                  nu(:,:,:),       'Nu_BL',          &
                                 atmos_grid_cartesc_real_vol(:,:,:), &
                                 atmos_grid_cartesc_real_totvol      )
          call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                 kh(:,:,:),         'Kh_BL',         &
                                 atmos_grid_cartesc_real_vol(:,:,:), &
                                 atmos_grid_cartesc_real_totvol      )

          do iq = 1, qa
             if ( .not. tracer_advc(iq) ) cycle
             call statistics_total( ka, ks, ke, ia, is, ie, ja, js, je, &
                                    rhoq_t_bl(:,:,:,iq), trim(tracer_name(iq))//'_t_BL', &
                                    atmos_grid_cartesc_real_vol(:,:,:),                  &
                                    atmos_grid_cartesc_real_totvol                       )
          enddo
       endif

    endif

    !$omp parallel do default(none) private(i,j,k) OMP_SCHEDULE_ collapse(2) &
    !$omp shared(JS,JE,IS,IE,KS,KE,RHOU_t,RHOU_t_BL,RHOV_t,RHOV_t_BL,RHOT_t,RHOT_t_BL)
    do j = js, je
    do i = is, ie
    do k = ks, ke
       rhou_t(k,i,j) = rhou_t(k,i,j) + rhou_t_bl(k,i,j)
       rhov_t(k,i,j) = rhov_t(k,i,j) + rhov_t_bl(k,i,j)
       rhot_t(k,i,j) = rhot_t(k,i,j) + rhot_t_bl(k,i,j)
    enddo
    enddo
    enddo

    do iq = 1,  qa
       if ( .not. tracer_advc(iq) ) cycle
       !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(3)
       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_bl(k,i,j,iq)
       enddo
       enddo
       enddo
    enddo

    return

Here is the caller graph for this function: