module atmosphere / physics / pbl / common More...

Functions/Subroutines
subroutine, public	atmos_phy_bl_tendency_tracer (KA, KS, KE, IA, IS, IE, JA, JS, JE, DENS, QTRC, SFLX_Q, Kh, MASS, CZ, FZ, F2H, DDT, TRACER_NAME, RHOQ_t)
	ATMOS_PHY_BL_tendency_tracer calculate tendency of tracers by the eddy viscosity. More...

Detailed Description

module atmosphere / physics / pbl / common

Description: Common routines for boundary layer turbulence

Author: Team SCALE

NAMELIST

No namelist group

History Output

name	description	unit	variable
ZFLX_{TRACER_NAME}_BL	Z FLUX of DENS * {TRACER_NAME} (PBL); {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	kg/m2/s	flx

Function/Subroutine Documentation

◆ atmos_phy_bl_tendency_tracer()

subroutine, public scale_atmos_phy_bl_common::atmos_phy_bl_tendency_tracer	(	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,
		real(rp), dimension (ka,ia,ja), intent(in)	DENS,
		real(rp), dimension (ka,ia,ja), intent(in)	QTRC,
		real(rp), dimension( ia,ja), intent(in)	SFLX_Q,
		real(rp), dimension (ka,ia,ja), intent(in)	Kh,
		real(rp), intent(in)	MASS,
		real(rp), dimension ( ka,ia,ja), intent(in)	CZ,
		real(rp), dimension (0:ka,ia,ja), intent(in)	FZ,
		real(rp), dimension(ka,2,ia,ja), intent(in)	F2H,
		real(dp), intent(in)	DDT,
		character(len=*), intent(in)	TRACER_NAME,
		real(rp), dimension(ka,ia,ja), intent(out)	RHOQ_t
	)

ATMOS_PHY_BL_tendency_tracer calculate tendency of tracers by the eddy viscosity.

Definition at line 62 of file scale_atmos_phy_bl_common.F90.

     use scale_matrix, only: &
        matrix_solver_tridiagonal
     use scale_file_history, only: &
        file_history_in
     integer, intent(in) :: KA, KS, KE
     integer, intent(in) :: IA, IS, IE
     integer, intent(in) :: JA, JS, JE
  
     real(RP),         intent(in) :: DENS  (KA,IA,JA)
     real(RP),         intent(in) :: QTRC  (KA,IA,JA)
     real(RP),         intent(in) :: SFLX_Q(   IA,JA)
     real(RP),         intent(in) :: Kh    (KA,IA,JA)
     real(RP),         intent(in) :: MASS
     real(RP),         intent(in) :: CZ (  KA,IA,JA)
     real(RP),         intent(in) :: FZ (0:KA,IA,JA)
     real(RP),         intent(in) :: F2H(KA,2,IA,JA)
     real(DP),         intent(in) :: DDT
     character(len=*), intent(in) :: TRACER_NAME
  
     real(RP), intent(out) :: RHOQ_t(KA,IA,JA)
  
     real(RP) :: QTRC_n(KA)
     real(RP) :: RHO   (KA)
     real(RP) :: RHOKh (KA)
     real(RP) :: a(KA)
     real(RP) :: b(KA)
     real(RP) :: c(KA)
     real(RP) :: d(KA)
     real(RP) :: rho_h
     real(RP) :: ap
     real(RP) :: sf_t
  
     real(RP) :: flx(KA,IA,JA)
  
     real(RP) :: CDZ(KA)
     real(RP) :: FDZ(KA)
  
     real(RP) :: dt
  
     integer :: KE_PBL
     integer :: k, i, j
  
 #ifdef _OPENACC
     real(RP) :: work(KS:KE,4)
 #endif
     dt = real( ddt, kind=rp )
  
     !$acc data create(flx)
  
 !OCL INDEPENDENT
     !$omp parallel do default(none) OMP_SCHEDULE_ collapse(2) &
     !$omp shared(KA,KS,KE,IS,IE,JS,JE) &
     !$omp shared(RHOQ_t,DENS,QTRC,SFLX_Q,Kh,MASS,CZ,FZ,F2H,DT,flx) &
     !$omp private(QTRC_n,RHO,RHOKh,rho_h,a,b,c,d,ap,sf_t,CDZ,FDZ) &
     !$omp private(KE_PBL,k,i,j)
     !$acc kernels
     !$acc loop independent collapse(2) &
     !$acc private(QTRC_n,RHO,RHOKh,rho_h,a,b,c,d,ap,sf_t,CDZ,FDZ, &
     !$acc         KE_PBL, &
     !$acc         work)
     do j = js, je
     do i = is, ie
  
        ke_pbl = ke-1
        !$acc loop seq
        do k = ke-2, ks+1, -1
           if ( kh(k,i,j) > 0.0_rp ) then
              ke_pbl = k + 1
           else
              exit
           end if
        end do
  
        do k = ks, ke_pbl
           cdz(k) = fz(k  ,i,j) - fz(k-1,i,j)
           fdz(k) = cz(k+1,i,j) - cz(k  ,i,j)
        end do
  
        sf_t = sflx_q(i,j) / cdz(ks)
  
        rho(ks) = dens(ks,i,j) + dt * sf_t * mass
        do k = ks+1, ke_pbl
           rho(k) = dens(k,i,j)
        end do
  
        ! dens * coefficient at the half level
        do k = ks, ke_pbl-1
           rho_h = f2h(k,1,i,j) * dens(k+1,i,j) + f2h(k,2,i,j) * dens(k,i,j)
           rhokh(k) = rho_h * kh(k,i,j)
        end do
  
        d(ks) = ( qtrc(ks,i,j) * dens(ks,i,j) + dt * sf_t ) / rho(ks)
        do k = ks+1, ke_pbl
           d(k) = qtrc(k,i,j)
        end do
  
        c(ks) = 0.0_rp
        do k = ks, ke_pbl-1
           ap = - dt * rhokh(k) / fdz(k)
           a(k) = ap / ( rho(k) * cdz(k) )
 #ifdef _OPENACC
           if ( k==ks ) then
              b(k) = - a(k) + 1.0_rp
           else
              b(k) = - a(k) + dt * rhokh(k-1) / ( fdz(k-1) * rho(k) * cdz(k) ) + 1.0_rp
           end if
 #else
           b(k) = - a(k) - c(k) + 1.0_rp
 #endif
           c(k+1) = ap / ( rho(k+1) * cdz(k+1) )
        end do
        a(ke_pbl) = 0.0_rp
        b(ke_pbl) = - c(ke_pbl) + 1.0_rp
  
        call matrix_solver_tridiagonal( &
                ka, ks, ke_pbl, &
 #ifdef _OPENACC
                work(:,:), & ! (wrok)
 #endif
                a(:), b(:), c(:), d(:), & ! (in)
                qtrc_n(:)               ) ! (out)
  
        rhoq_t(ks,i,j) = ( qtrc_n(ks) * rho(ks) - qtrc(ks,i,j) * dens(ks,i,j) ) / dt - sf_t
        do k = ks+1, ke_pbl
           rhoq_t(k,i,j) = ( qtrc_n(k) - qtrc(k,i,j) ) * rho(k) / dt
        end do
        do k = ke_pbl+1, ke
           rhoq_t(k,i,j) = 0.0_rp
        end do
  
        do k = ks, ke_pbl-1
           flx(k,i,j) = - rhokh(k) * ( qtrc_n(k+1) - qtrc_n(k) ) / fdz(k)
        end do
  
     end do
     end do
     !$acc end kernels
  
     call file_history_in(flx(:,:,:), 'ZFLX_'//trim(tracer_name)//'_BL', 'Z FLUX of DENS * '//trim(tracer_name)//' (PBL)', 'kg/m2/s', fill_halo=.true.)
  
     !$acc end data
  
     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::ka, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_tracer::mass, scale_precision::rp, and scale_tracer::tracer_name.

Referenced by mod_atmos_phy_bl_driver::atmos_phy_bl_driver_calc_tendency().

Here is the caller graph for this function:

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ atmos_phy_bl_tendency_tracer()