mod_atmos_phy_tb_driver Module Reference

module ATMOSPHERE / Physics Turbulence More...

Functions/Subroutines
subroutine, public	atmos_phy_tb_driver_setup
	Setup. More...
subroutine, public	atmos_phy_tb_driver_resume
	Resume. More...
subroutine, public	atmos_phy_tb_driver (update_flag)
	Driver. More...

Detailed Description

module ATMOSPHERE / Physics Turbulence

Description

Sub-grid scale turbulence process

Author

Team SCALE

History

• 2013-12-05 (S.Nishizawa) [new]

NAMELIST

• No namelist group

History Output

name	description	unit	variable
MOMX_t_TB	MOMX tendency (TB)	kg/m2/s2	MOMX_t_TB
MOMY_t_TB	MOMY tendency (TB)	kg/m2/s2	MOMY_t_TB
MOMZ_t_TB	MOMZ tendency (TB)	kg/m2/s2	MOMZ_t_TB
N2	squared Brunt-Vaisala frequency	1/s2	N2
NU	eddy viscosity	m2/s	NU
Pr	Prantle number	NIL	Pr
RHOT_t_TB	RHOT tendency (TB)	K.kg/m3/s	RHOT_t_TB
Ri	Richardson number	NIL	Ri
SGS_XFLX_MOMX	SGS X FLUX of MOMX	kg/m/s2	QFLX_MOMX
SGS_XFLX_MOMY	SGS X FLUX of MOMY	kg/m/s2	QFLX_MOMY
SGS_XFLX_MOMZ	SGS X FLUX of MOMZ	kg/m/s2	QFLX_MOMZ
SGS_XFLX_QC	SGS X FLUX of QC	kg/m2/s	QFLX_RHOQ
SGS_XFLX_QG	SGS X FLUX of QG	kg/m2/s	QFLX_RHOQ
SGS_XFLX_QI	SGS X FLUX of QI	kg/m2/s	QFLX_RHOQ
SGS_XFLX_QR	SGS X FLUX of QR	kg/m2/s	QFLX_RHOQ
SGS_XFLX_QS	SGS X FLUX of QS	kg/m2/s	QFLX_RHOQ
SGS_XFLX_QV	SGS X FLUX of QV	kg/m2/s	QFLX_RHOQ
SGS_XFLX_RHOT	SGS X FLUX of RHOT	K*kg/m2/s	QFLX_RHOT
SGS_YFLX_MOMX	SGS Y FLUX of MOMX	kg/m/s2	QFLX_MOMX
SGS_YFLX_MOMY	SGS Y FLUX of MOMY	kg/m/s2	QFLX_MOMY
SGS_YFLX_MOMZ	SGS Y FLUX of MOMZ	kg/m/s2	QFLX_MOMZ
SGS_YFLX_QC	SGS Y FLUX of QC	kg/m2/s	QFLX_RHOQ
SGS_YFLX_QG	SGS Y FLUX of QG	kg/m2/s	QFLX_RHOQ
SGS_YFLX_QI	SGS Y FLUX of QI	kg/m2/s	QFLX_RHOQ
SGS_YFLX_QR	SGS Y FLUX of QR	kg/m2/s	QFLX_RHOQ
SGS_YFLX_QS	SGS Y FLUX of QS	kg/m2/s	QFLX_RHOQ
SGS_YFLX_QV	SGS Y FLUX of QV	kg/m2/s	QFLX_RHOQ
SGS_YFLX_RHOT	SGS Y FLUX of RHOT	K*kg/m2/s	QFLX_RHOT
SGS_ZFLX_MOMX	SGS Z FLUX of MOMX	kg/m/s2	QFLX_MOMX
SGS_ZFLX_MOMY	SGS Z FLUX of MOMY	kg/m/s2	QFLX_MOMY
SGS_ZFLX_MOMZ	SGS Z FLUX of MOMZ	kg/m/s2	QFLX_MOMZ
SGS_ZFLX_QC	SGS Z FLUX of QC	kg/m2/s	QFLX_RHOQ
SGS_ZFLX_QG	SGS Z FLUX of QG	kg/m2/s	QFLX_RHOQ
SGS_ZFLX_QI	SGS Z FLUX of QI	kg/m2/s	QFLX_RHOQ
SGS_ZFLX_QR	SGS Z FLUX of QR	kg/m2/s	QFLX_RHOQ
SGS_ZFLX_QS	SGS Z FLUX of QS	kg/m2/s	QFLX_RHOQ
SGS_ZFLX_QV	SGS Z FLUX of QV	kg/m2/s	QFLX_RHOQ
SGS_ZFLX_RHOT	SGS Z FLUX of RHOT	K*kg/m2/s	QFLX_RHOT
TKE	turburent kinetic energy	m2/s2	TKE
trim(AQ_NAME(iq))//'_t_TB'	RHO*'//trim(AQ_NAME(iq))//' tendency (TB)	kg/m3/s	RHOQ_t_TB

Function/Subroutine Documentation

atmos_phy_tb_driver_setup()

subroutine, public mod_atmos_phy_tb_driver::atmos_phy_tb_driver_setup

(

)

Setup.

Definition at line 52 of file mod_atmos_phy_tb_driver.f90.

References mod_atmos_phy_tb_vars::atmos_phy_tb_momz_t, scale_atmos_phy_tb::atmos_phy_tb_setup(), mod_atmos_phy_tb_vars::atmos_phy_tb_tke_t, mod_atmos_admin::atmos_phy_tb_type, mod_atmos_admin::atmos_sw_phy_tb, scale_grid::grid_cdx, scale_grid::grid_cdy, scale_grid::grid_cdz, scale_grid_index::ie, scale_stdio::io_fid_log, scale_stdio::io_l, scale_grid_index::is, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ke, scale_grid_index::ks, scale_process::prc_mpistop(), and scale_grid_real::real_cz.

Referenced by mod_atmos_driver::atmos_driver_setup().

    use scale_grid, only: &
       cdz => grid_cdz, &
       cdx => grid_cdx, &
       cdy => grid_cdy
    use scale_grid_real, only: &
       cz  => real_cz
    use scale_process, only: &
       prc_mpistop
    use mod_atmos_admin, only: &
       atmos_phy_tb_type, &
       atmos_sw_phy_tb
    use scale_atmos_phy_tb, only: &
       atmos_phy_tb_setup
    use mod_atmos_phy_tb_vars, only: &
       tke_t_tb  => atmos_phy_tb_tke_t,  &
       momz_t_tb => atmos_phy_tb_momz_t
    implicit none

!    NAMELIST / PARAM_ATMOS_PHY_TB / &

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

    if( io_l ) write(io_fid_log,*)
    if( io_l ) write(io_fid_log,*) '++++++ Module[DRIVER] / Categ[ATMOS PHY_TB] / Origin[SCALE-RM]'

!    !--- read namelist
!    rewind(IO_FID_CONF)
!    read(IO_FID_CONF,nml=PARAM_ATMOS_PHY_TB,iostat=ierr)
!    if( ierr < 0 ) then !--- missing
!       if( IO_L ) write(IO_FID_LOG,*) '*** Not found namelist. Default used.'
!    elseif( ierr > 0 ) then !--- fatal error
!       write(*,*) 'xxx Not appropriate names in namelist PARAM_ATMOS_PHY_TB. Check!'
!       call PRC_MPIstop
!    endif
!    if( IO_LNML ) write(IO_FID_LOG,nml=PARAM_ATMOS_PHY_TB)


    ! initialize
    do j = js, je
    do i = is, ie
       momz_t_tb(ks-1,i,j) = 0.0_rp
       momz_t_tb(ke  ,i,j) = 0.0_rp
    enddo
    enddo

    do j = js, je
    do i = is, ie
    do k = ks, ke
       tke_t_tb(k,i,j) = 0.0_rp
    enddo
    enddo
    enddo

    if ( atmos_sw_phy_tb ) then

       ! setup library component
       call atmos_phy_tb_setup( atmos_phy_tb_type, & ! [IN]
                                cdz, cdx, cdy, cz  ) ! [IN]
    else
       if( io_l ) write(io_fid_log,*) '*** this component is never called.'
    endif

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_tb_driver_resume()

subroutine, public mod_atmos_phy_tb_driver::atmos_phy_tb_driver_resume

(

)

Resume.

Definition at line 122 of file mod_atmos_phy_tb_driver.f90.

References atmos_phy_tb_driver(), mod_atmos_admin::atmos_sw_phy_tb, scale_prof::prof_rapend(), and scale_prof::prof_rapstart().

Referenced by mod_atmos_driver::atmos_driver_resume2().

    use mod_atmos_admin, only: &
       atmos_sw_phy_tb
    implicit none

    if ( atmos_sw_phy_tb ) then

       ! run once (only for the diagnostic value)
       call prof_rapstart('ATM_Turbulence', 1)
       call atmos_phy_tb_driver( update_flag = .true. )
       call prof_rapend  ('ATM_Turbulence', 1)

    end if

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_tb_driver()

subroutine, public mod_atmos_phy_tb_driver::atmos_phy_tb_driver

(

logical, intent(in)

update_flag

)

Driver.

Definition at line 141 of file mod_atmos_phy_tb_driver.f90.

References scale_tracer::aq_name, 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_mw, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_qtrc, mod_atmos_phy_sf_vars::atmos_phy_sf_sflx_sh, scale_atmos_phy_tb::atmos_phy_tb, scale_atmos_phy_tb_common::atmos_phy_tb_calc_tend_momx(), scale_atmos_phy_tb_common::atmos_phy_tb_calc_tend_momy(), scale_atmos_phy_tb_common::atmos_phy_tb_calc_tend_momz(), scale_atmos_phy_tb_common::atmos_phy_tb_calc_tend_phi(), mod_atmos_phy_tb_vars::atmos_phy_tb_momx_t, mod_atmos_phy_tb_vars::atmos_phy_tb_momy_t, mod_atmos_phy_tb_vars::atmos_phy_tb_momz_t, mod_atmos_phy_tb_vars::atmos_phy_tb_nu, mod_atmos_phy_tb_vars::atmos_phy_tb_rhoq_t, mod_atmos_phy_tb_vars::atmos_phy_tb_rhot_t, mod_atmos_phy_tb_vars::atmos_phy_tb_tke, mod_atmos_phy_tb_vars::atmos_phy_tb_tke_t, mod_atmos_vars::dens, mod_atmos_vars::dens_av, scale_gridtrans::gtrans_gsqrt, scale_gridtrans::gtrans_j13g, scale_gridtrans::gtrans_j23g, scale_gridtrans::gtrans_j33g, scale_gridtrans::gtrans_mapf, scale_tracer::i_qc, scale_tracer::i_qg, scale_tracer::i_qi, scale_tracer::i_qr, scale_tracer::i_qs, scale_tracer::i_qv, scale_grid_index::iblock, scale_grid_index::ie, scale_grid_index::is, scale_grid_index::jblock, scale_grid_index::je, scale_grid_index::js, scale_grid_index::ka, scale_grid_index::ke, scale_grid_index::ks, mod_atmos_vars::momx, mod_atmos_vars::momx_av, mod_atmos_vars::momx_tp, mod_atmos_vars::momy, mod_atmos_vars::momy_av, mod_atmos_vars::momy_tp, mod_atmos_vars::momz, mod_atmos_vars::momz_av, mod_atmos_vars::momz_tp, scale_tracer::qa, mod_atmos_vars::qtrc, mod_atmos_vars::qtrc_av, mod_atmos_vars::rhoq_tp, mod_atmos_vars::rhot, mod_atmos_vars::rhot_av, mod_atmos_vars::rhot_tp, scale_rm_statistics::statistics_checktotal, scale_time::time_dtsec_atmos_phy_tb, scale_grid_index::xdir, scale_grid_index::ydir, and scale_grid_index::zdir.

Referenced by mod_atmos_driver::atmos_driver(), and atmos_phy_tb_driver_resume().

    use scale_gridtrans, only: &
       gsqrt => gtrans_gsqrt, &
       j13g  => gtrans_j13g,  &
       j23g  => gtrans_j23g,  &
       j33g  => gtrans_j33g,  &
       mapf  => gtrans_mapf
    use scale_rm_statistics, only: &
       statistics_checktotal, &
       stat_total
    use scale_comm, only: &
       comm_vars8, &
       comm_wait
    use scale_history, only: &
       hist_in
    use scale_time, only: &
       dt_tb => time_dtsec_atmos_phy_tb
    use scale_atmos_phy_tb, only: &
       atmos_phy_tb
    use scale_atmos_phy_tb_common, only: &
       calc_tend_momz => atmos_phy_tb_calc_tend_momz, &
       calc_tend_momx => atmos_phy_tb_calc_tend_momx, &
       calc_tend_momy => atmos_phy_tb_calc_tend_momy, &
       calc_tend_phi  => atmos_phy_tb_calc_tend_phi
    use mod_atmos_vars, only: &
       dens => dens_av,   &
       momz => momz_av,   &
       momx => momx_av,   &
       momy => momy_av,   &
       rhot => rhot_av,   &
       qtrc => qtrc_av,   &
       momz_t => momz_tp, &
       momx_t => momx_tp, &
       momy_t => momy_tp, &
       rhot_t => rhot_tp, &
       rhoq_t => rhoq_tp
    use mod_atmos_phy_tb_vars, only: &
       momz_t_tb => atmos_phy_tb_momz_t, &
       momx_t_tb => atmos_phy_tb_momx_t, &
       momy_t_tb => atmos_phy_tb_momy_t, &
       rhot_t_tb => atmos_phy_tb_rhot_t, &
       rhoq_t_tb => atmos_phy_tb_rhoq_t, &
       tke_t_tb  => atmos_phy_tb_tke_t,  &
       tke       => atmos_phy_tb_tke,    &
       nu        => atmos_phy_tb_nu
    use mod_atmos_phy_sf_vars, only: &
       sflx_mw => atmos_phy_sf_sflx_mw, &
       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
    implicit none

    logical, intent(in) :: update_flag

    ! eddy viscosity/diffusion flux
    real(RP) :: qflx_momz(ka,ia,ja,3)
    real(RP) :: qflx_momx(ka,ia,ja,3)
    real(RP) :: qflx_momy(ka,ia,ja,3)
    real(RP) :: qflx_rhot(ka,ia,ja,3)
    real(RP) :: qflx_rhoq(ka,ia,ja,3,qa)

    real(RP) :: ri(ka,ia,ja) ! Richardson number
    real(RP) :: pr(ka,ia,ja) ! Prandtl number
    real(RP) :: n2(ka,ia,ja) ! squared Brunt-Vaisala frequency

    integer :: jjs, jje
    integer :: iis, iie

    real(RP) :: total ! dummy

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

    if ( update_flag ) then

       do j = js, je
       do i = is, ie
       do k = ks, ke
          tke(k,i,j) = max( tke(k,i,j) + tke_t_tb(k,i,j) * dt_tb, 0.0_rp )
       end do
       end do
       end do

       do j = js, je
       do i = is, ie
          tke(1:ks-1,i,j) = tke(ks,i,j)
          tke(ke+1:ka,i,j) = tke(ke,i,j)
       end do
       end do

       call comm_vars8( tke(:,:,:), 1 )
       call comm_wait ( tke(:,:,:), 1 )

       call atmos_phy_tb( qflx_momz, & ! [OUT]
                          qflx_momx, & ! [OUT]
                          qflx_momy, & ! [OUT]
                          qflx_rhot, & ! [OUT]
                          qflx_rhoq, & ! [OUT]
                          tke,       & ! [INOUT]
                          tke_t_tb,  & ! [OUT]
                          nu,        & ! [OUT]
                          ri,        & ! [OUT]
                          pr,        & ! [OUT]
                          n2,        & ! [OUT]
                          momz,      & ! [IN]
                          momx,      & ! [IN]
                          momy,      & ! [IN]
                          rhot,      & ! [IN]
                          dens,      & ! [IN]
                          qtrc,      & ! [IN]
                          sflx_mw,   & ! [IN]
                          sflx_mu,   & ! [IN]
                          sflx_mv,   & ! [IN]
                          sflx_sh,   & ! [IN]
                          sflx_q(:,:,i_qv), & ! [IN]
                          gsqrt,     & ! [IN]
                          j13g,      & ! [IN]
                          j23g,      & ! [IN]
                          j33g,      & ! [IN]
                          mapf,      & ! [IN]
                          dt_tb      ) ! [IN]

       do jjs = js, je, jblock
       jje = jjs+jblock-1
       do iis = is, ie, iblock
       iie = iis+iblock-1
          call calc_tend_momz( momz_t_tb,   & ! (out)
                               qflx_momz,   & ! (in)
                               gsqrt, j13g, j23g, j33g, mapf, & ! (in)
                               iis, iie, jjs, jje ) ! (in)
       end do
       end do

       do jjs = js, je, jblock
       jje = jjs+jblock-1
       do iis = is, ie, iblock
       iie = iis+iblock-1
          call calc_tend_momx( momx_t_tb,   & ! (out)
                               qflx_momx,   & ! (in)
                               gsqrt, j13g, j23g, j33g, mapf, & ! (in)
                               iis, iie, jjs, jje ) ! (in)
       end do
       end do

       do jjs = js, je, jblock
       jje = jjs+jblock-1
       do iis = is, ie, iblock
       iie = iis+iblock-1
          call calc_tend_momy( momy_t_tb,   & ! (out)
                               qflx_momy,   & ! (in)
                               gsqrt, j13g, j23g, j33g, mapf, & ! (in)
                               iis, iie, jjs, jje ) ! (in)
       end do
       end do

       do jjs = js, je, jblock
       jje = jjs+jblock-1
       do iis = is, ie, iblock
       iie = iis+iblock-1
          call calc_tend_phi ( rhot_t_tb,  & ! (out)
                               qflx_rhot,  & ! (in)
                               gsqrt, j13g, j23g, j33g, mapf, & ! (in)
                               iis, iie, jjs, jje ) ! (in)
       end do
       end do

       do iq = 1, qa
          do jjs = js, je, jblock
          jje = jjs+jblock-1
          do iis = is, ie, iblock
          iie = iis+iblock-1
             call calc_tend_phi ( rhoq_t_tb(:,:,:,iq),    & ! (out)
                                  qflx_rhoq(:,:,:,:,iq),  & ! (in)
                                  gsqrt, j13g, j23g, j33g, mapf, & ! (in)
                                  iis, iie, jjs, jje ) ! (in)
          end do
          end do
       end do

       call hist_in( tke(:,:,:), 'TKE', 'turburent kinetic energy', 'm2/s2', nohalo=.true. )
       call hist_in( nu(:,:,:), 'NU',  'eddy viscosity',           'm2/s' , nohalo=.true. )
       call hist_in( ri(:,:,:), 'Ri',  'Richardson number',        'NIL'  , nohalo=.true. )
       call hist_in( pr(:,:,:), 'Pr',  'Prantle number',           'NIL'  , nohalo=.true. )
       call hist_in( n2(:,:,:), 'N2',  'squared Brunt-Vaisala frequency', '1/s2', nohalo=.true. )

       call hist_in( momz_t_tb(:,:,:), 'MOMZ_t_TB', 'MOMZ tendency (TB)', 'kg/m2/s2',  nohalo=.true. )
       call hist_in( momx_t_tb(:,:,:), 'MOMX_t_TB', 'MOMX tendency (TB)', 'kg/m2/s2',  nohalo=.true. )
       call hist_in( momy_t_tb(:,:,:), 'MOMY_t_TB', 'MOMY tendency (TB)', 'kg/m2/s2',  nohalo=.true. )
       call hist_in( rhot_t_tb(:,:,:), 'RHOT_t_TB', 'RHOT tendency (TB)', 'K.kg/m3/s', nohalo=.true. )

       do iq = 1, qa
          call hist_in( rhoq_t_tb(:,:,:,iq), trim(aq_name(iq))//'_t_TB',                      &
                        'RHO*'//trim(aq_name(iq))//' tendency (TB)', 'kg/m3/s', nohalo=.true. )
       enddo

       call hist_in( qflx_momz(:,:,:,zdir), 'SGS_ZFLX_MOMZ', 'SGS Z FLUX of MOMZ', 'kg/m/s2', &
                     nohalo=.true.)
       call hist_in( qflx_momz(:,:,:,xdir), 'SGS_XFLX_MOMZ', 'SGS X FLUX of MOMZ', 'kg/m/s2', &
                     xdim='half', zdim='half', nohalo=.true.)
       call hist_in( qflx_momz(:,:,:,ydir), 'SGS_YFLX_MOMZ', 'SGS Y FLUX of MOMZ', 'kg/m/s2', &
                     ydim='half', zdim='half', nohalo=.true.)

       call hist_in( qflx_momx(:,:,:,zdir), 'SGS_ZFLX_MOMX', 'SGS Z FLUX of MOMX', 'kg/m/s2', &
                     xdim='half', zdim='half', nohalo=.true.)
       call hist_in( qflx_momx(:,:,:,xdir), 'SGS_XFLX_MOMX', 'SGS X FLUX of MOMX', 'kg/m/s2', &
                     nohalo=.true.)
       call hist_in( qflx_momx(:,:,:,ydir), 'SGS_YFLX_MOMX', 'SGS Y FLUX of MOMX', 'kg/m/s2', &
                     xdim='half', ydim='half', nohalo=.true.)

       call hist_in( qflx_momy(:,:,:,zdir), 'SGS_ZFLX_MOMY', 'SGS Z FLUX of MOMY', 'kg/m/s2', &
                     ydim='half', zdim='half', nohalo=.true.)
       call hist_in( qflx_momy(:,:,:,xdir), 'SGS_XFLX_MOMY', 'SGS X FLUX of MOMY', 'kg/m/s2', &
                     xdim='half', ydim='half', nohalo=.true.)
       call hist_in( qflx_momy(:,:,:,ydir), 'SGS_YFLX_MOMY', 'SGS Y FLUX of MOMY', 'kg/m/s2', &
                     nohalo=.true.)

       call hist_in( qflx_rhot(:,:,:,zdir), 'SGS_ZFLX_RHOT', 'SGS Z FLUX of RHOT', 'K*kg/m2/s', &
                     zdim='half', nohalo=.true.)
       call hist_in( qflx_rhot(:,:,:,xdir), 'SGS_XFLX_RHOT', 'SGS X FLUX of RHOT', 'K*kg/m2/s', &
                     xdim='half', nohalo=.true.)
       call hist_in( qflx_rhot(:,:,:,ydir), 'SGS_YFLX_RHOT', 'SGS Y FLUX of RHOT', 'K*kg/m2/s', &
                     ydim='half', nohalo=.true.)

       if ( i_qv > 0 ) then
          call hist_in( qflx_rhoq(:,:,:,zdir,i_qv), 'SGS_ZFLX_QV', 'SGS Z FLUX of QV', 'kg/m2/s', &
                        zdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,xdir,i_qv), 'SGS_XFLX_QV', 'SGS X FLUX of QV', 'kg/m2/s', &
                        xdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,ydir,i_qv), 'SGS_YFLX_QV', 'SGS Y FLUX of QV', 'kg/m2/s', &
                        ydim='half', nohalo=.true.)
       endif

#ifndef DRY
       if ( i_qc > 0 ) then
          call hist_in( qflx_rhoq(:,:,:,zdir,i_qc), 'SGS_ZFLX_QC', 'SGS Z FLUX of QC', 'kg/m2/s', &
                        zdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,xdir,i_qc), 'SGS_XFLX_QC', 'SGS X FLUX of QC', 'kg/m2/s', &
                        xdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,ydir,i_qc), 'SGS_YFLX_QC', 'SGS Y FLUX of QC', 'kg/m2/s', &
                        ydim='half', nohalo=.true.)
       endif

       if ( i_qr > 0 ) then
          call hist_in( qflx_rhoq(:,:,:,zdir,i_qr), 'SGS_ZFLX_QR', 'SGS Z FLUX of QR', 'kg/m2/s', &
                        zdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,xdir,i_qr), 'SGS_XFLX_QR', 'SGS X FLUX of QR', 'kg/m2/s', &
                        xdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,ydir,i_qr), 'SGS_YFLX_QR', 'SGS Y FLUX of QR', 'kg/m2/s', &
                        ydim='half', nohalo=.true.)
       endif

       if ( i_qi > 0 ) then
          call hist_in( qflx_rhoq(:,:,:,zdir,i_qi), 'SGS_ZFLX_QI', 'SGS Z FLUX of QI', 'kg/m2/s', &
                        zdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,xdir,i_qi), 'SGS_XFLX_QI', 'SGS X FLUX of QI', 'kg/m2/s', &
                        xdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,ydir,i_qi), 'SGS_YFLX_QI', 'SGS Y FLUX of QI', 'kg/m2/s', &
                        ydim='half', nohalo=.true.)
       endif

       if ( i_qs > 0 ) then
          call hist_in( qflx_rhoq(:,:,:,zdir,i_qs), 'SGS_ZFLX_QS', 'SGS Z FLUX of QS', 'kg/m2/s', &
                        zdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,xdir,i_qs), 'SGS_XFLX_QS', 'SGS X FLUX of QS', 'kg/m2/s', &
                        xdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,ydir,i_qs), 'SGS_YFLX_QS', 'SGS Y FLUX of QS', 'kg/m2/s', &
                        ydim='half', nohalo=.true.)
       endif

       if ( i_qg > 0 ) then
          call hist_in( qflx_rhoq(:,:,:,zdir,i_qg), 'SGS_ZFLX_QG', 'SGS Z FLUX of QG', 'kg/m2/s', &
                        zdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,xdir,i_qg), 'SGS_XFLX_QG', 'SGS X FLUX of QG', 'kg/m2/s', &
                        xdim='half', nohalo=.true.)
          call hist_in( qflx_rhoq(:,:,:,ydir,i_qg), 'SGS_YFLX_QG', 'SGS Y FLUX of QG', 'kg/m2/s', &
                        ydim='half', nohalo=.true.)
       endif
#endif

       if ( statistics_checktotal ) then
          call stat_total( total, momz_t_tb(:,:,:), 'MOMZ_t_TB' )
          call stat_total( total, momx_t_tb(:,:,:), 'MOMX_t_TB' )
          call stat_total( total, momy_t_tb(:,:,:), 'MOMY_t_TB' )
          call stat_total( total, rhot_t_tb(:,:,:), 'RHOT_t_TB' )
          call stat_total( total, tke(:,:,:), 'TKE' )
          call stat_total( total, nu(:,:,:), 'Nu' )
          call stat_total( total, ri(:,:,:), 'Ri' )
          call stat_total( total, pr(:,:,:), 'Pr' )

          do iq = 1, qa
             call stat_total( total, rhoq_t_tb(:,:,:,iq), trim(aq_name(iq))//'_t_TB' )
          enddo
       endif

    endif

    !$omp parallel do private(i,j,k) OMP_SCHEDULE_ collapse(2)
    do j = js, je
    do i = is, ie
    do k = ks, ke
       momz_t(k,i,j) = momz_t(k,i,j) + momz_t_tb(k,i,j)
       momx_t(k,i,j) = momx_t(k,i,j) + momx_t_tb(k,i,j)
       momy_t(k,i,j) = momy_t(k,i,j) + momy_t_tb(k,i,j)
       rhot_t(k,i,j) = rhot_t(k,i,j) + rhot_t_tb(k,i,j)
    enddo
    enddo
    enddo

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

    return

Here is the call graph for this function:

Here is the caller graph for this function: