scale_atmos_phy_tb_dns Module Reference

module ATMOSPHERE / Physics Turbulence More...

Functions/Subroutines
subroutine, public	atmos_phy_tb_dns_config (TYPE_TB, I_TKE_out)
	Config. More...
subroutine, public	atmos_phy_tb_dns_setup (CDZ, CDX, CDY, CZ)
subroutine, public	atmos_phy_tb_dns (qflx_sgs_MOMZ, qflx_sgs_MOMX, qflx_sgs_MOMY, qflx_sgs_rhot, qflx_sgs_rhoq, RHOQ_t, nu, Ri, Pr, MOMZ, MOMX, MOMY, RHOT, DENS, QTRC, N2, SFLX_MW, SFLX_MU, SFLX_MV, SFLX_SH, SFLX_Q, GSQRT, J13G, J23G, J33G, MAPF, dt)

Detailed Description

module ATMOSPHERE / Physics Turbulence

Description

Sub-grid scale turbulent process for DNS

Author

Team SCALE

NAMELIST

• PARAM_ATMOS_PHY_TB_DNS

  name	type	default value	comment
  ATMOS_PHY_TB_DNS_NU	real(RP)	1.512E-5_RP	[m2/s] kinematic viscosity coefficient for air at 20degC
  ATMOS_PHY_TB_DNS_MU	real(RP)	1.512E-5_RP	same as NU (needed based on hyposes. see Mellado 2010)

  
  

History Output

No history output

Function/Subroutine Documentation

atmos_phy_tb_dns_config()

subroutine, public scale_atmos_phy_tb_dns::atmos_phy_tb_dns_config	(	character(len=*), intent(in)	TYPE_TB,
		integer, intent(out)	I_TKE_out
	)

Config.

Definition at line 64 of file scale_atmos_phy_tb_dns.F90.

References scale_prc::prc_abort().

Referenced by mod_atmos_phy_tb_driver::atmos_phy_tb_driver_tracer_setup().

    use scale_prc, only: &
       prc_abort
    implicit none

    character(len=*), intent(in)  :: type_tb
    integer,          intent(out) :: i_tke_out
    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_PHY_TB_dns_config",*) 'Setup'
    log_info("ATMOS_PHY_TB_dns_config",*) 'Tracers for Deardorff (1980) 1.5th TKE Model'

    if ( type_tb /= 'DNS' ) then
       log_error("ATMOS_PHY_TB_dns_config",*) 'ATMOS_PHY_TB_TYPE is not DNS. Check!'
       call prc_abort
    endif

    i_tke_out = -1

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_tb_dns_setup()

subroutine, public scale_atmos_phy_tb_dns::atmos_phy_tb_dns_setup	(	real(rp), dimension(ka), intent(in)	CDZ,
		real(rp), dimension(ia), intent(in)	CDX,
		real(rp), dimension(ja), intent(in)	CDY,
		real(rp), dimension (ka,ia,ja), intent(in)	CZ
	)

Definition at line 88 of file scale_atmos_phy_tb_dns.F90.

References scale_io::io_fid_conf, and scale_prc::prc_abort().

Referenced by mod_atmos_phy_tb_driver::atmos_phy_tb_driver_setup().

    use scale_prc, only: &
       prc_abort
    implicit none

    real(RP), intent(in) :: cdz(ka)
    real(RP), intent(in) :: cdx(ia)
    real(RP), intent(in) :: cdy(ja)
    real(RP), intent(in) :: cz (ka,ia,ja)

    namelist / param_atmos_phy_tb_dns / &
       atmos_phy_tb_dns_nu, &
       atmos_phy_tb_dns_mu

    integer :: ierr
    !---------------------------------------------------------------------------

    log_newline
    log_info("ATMOS_PHY_TB_dns_setup",*) 'Setup'
    log_info("ATMOS_PHY_TB_dns_setup",*) 'Eddy Viscocity Model for DNS'

    !--- read namelist
    rewind(io_fid_conf)
    read(io_fid_conf,nml=param_atmos_phy_tb_dns,iostat=ierr)
    if( ierr < 0 ) then !--- missing
       log_info("ATMOS_PHY_TB_dns_setup",*) 'Not found namelist. Default used.'
    elseif( ierr > 0 ) then !--- fatal error
       log_error("ATMOS_PHY_TB_dns_setup",*) 'Not appropriate names in namelist PARAM_ATMOS_PHY_TB_DNS. Check!'
       call prc_abort
    endif
    log_nml(param_atmos_phy_tb_dns)

    return

Here is the call graph for this function:

Here is the caller graph for this function:

atmos_phy_tb_dns()

subroutine, public scale_atmos_phy_tb_dns::atmos_phy_tb_dns	(	real(rp), dimension(ka,ia,ja,3), intent(out)	qflx_sgs_MOMZ,
		real(rp), dimension(ka,ia,ja,3), intent(out)	qflx_sgs_MOMX,
		real(rp), dimension(ka,ia,ja,3), intent(out)	qflx_sgs_MOMY,
		real(rp), dimension(ka,ia,ja,3), intent(out)	qflx_sgs_rhot,
		real(rp), dimension(ka,ia,ja,3,qa), intent(out)	qflx_sgs_rhoq,
		real(rp), dimension(ka,ia,ja,qa), intent(inout)	RHOQ_t,
		real(rp), dimension(ka,ia,ja), intent(out)	nu,
		real(rp), dimension(ka,ia,ja), intent(out)	Ri,
		real(rp), dimension(ka,ia,ja), intent(out)	Pr,
		real(rp), dimension(ka,ia,ja), intent(in)	MOMZ,
		real(rp), dimension(ka,ia,ja), intent(in)	MOMX,
		real(rp), dimension(ka,ia,ja), intent(in)	MOMY,
		real(rp), dimension(ka,ia,ja), intent(in)	RHOT,
		real(rp), dimension(ka,ia,ja), intent(in)	DENS,
		real(rp), dimension(ka,ia,ja,qa), intent(in)	QTRC,
		real(rp), dimension(ka,ia,ja), intent(in)	N2,
		real(rp), dimension(ia,ja), intent(in)	SFLX_MW,
		real(rp), dimension(ia,ja), intent(in)	SFLX_MU,
		real(rp), dimension(ia,ja), intent(in)	SFLX_MV,
		real(rp), dimension(ia,ja), intent(in)	SFLX_SH,
		real(rp), dimension (ia,ja,qa), intent(in)	SFLX_Q,
		real(rp), dimension (ka,ia,ja,7), intent(in)	GSQRT,
		real(rp), dimension (ka,ia,ja,7), intent(in)	J13G,
		real(rp), dimension (ka,ia,ja,7), intent(in)	J23G,
		real(rp), intent(in)	J33G,
		real(rp), dimension (ia,ja,2,4), intent(in)	MAPF,
		real(dp), intent(in)	dt
	)

Parameters

[in]	gsqrt	vertical metrics {G}^1/2
[in]	j13g	(1,3) element of Jacobian matrix
[in]	j23g	(1,3) element of Jacobian matrix
[in]	j33g	(3,3) element of Jacobian matrix
[in]	mapf	map factor

Definition at line 130 of file scale_atmos_phy_tb_dns.F90.

References scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdx, scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdy, scale_atmos_grid_cartesc::atmos_grid_cartesc_rcdz, scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdx, scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdy, scale_atmos_grid_cartesc::atmos_grid_cartesc_rfdz, scale_const::const_grav, scale_atmos_grid_cartesc_index::i_uy, scale_atmos_grid_cartesc_index::i_xv, scale_atmos_grid_cartesc_index::i_xy, scale_atmos_grid_cartesc_index::iblock, scale_atmos_grid_cartesc_index::ie, scale_atmos_grid_cartesc_index::is, scale_atmos_grid_cartesc_index::jblock, scale_atmos_grid_cartesc_index::je, scale_atmos_grid_cartesc_index::js, scale_atmos_grid_cartesc_index::ke, scale_atmos_grid_cartesc_index::ks, scale_tracer::qa, scale_tracer::tracer_advc, scale_atmos_grid_cartesc_index::xdir, scale_atmos_grid_cartesc_index::ydir, and scale_atmos_grid_cartesc_index::zdir.

Referenced by mod_atmos_phy_tb_driver::atmos_phy_tb_driver_calc_tendency().

    use scale_atmos_grid_cartesc_index
    use scale_tracer
    use scale_const, only: &
       grav => const_grav
    use scale_atmos_grid_cartesc, only: &
       rcdz => atmos_grid_cartesc_rcdz, &
       rcdx => atmos_grid_cartesc_rcdx, &
       rcdy => atmos_grid_cartesc_rcdy, &
       rfdz => atmos_grid_cartesc_rfdz, &
       rfdx => atmos_grid_cartesc_rfdx, &
       rfdy => atmos_grid_cartesc_rfdy
    implicit none

    ! SGS flux
    real(RP), intent(out) :: qflx_sgs_momz(ka,ia,ja,3)
    real(RP), intent(out) :: qflx_sgs_momx(ka,ia,ja,3)
    real(RP), intent(out) :: qflx_sgs_momy(ka,ia,ja,3)
    real(RP), intent(out) :: qflx_sgs_rhot(ka,ia,ja,3)
    real(RP), intent(out) :: qflx_sgs_rhoq(ka,ia,ja,3,qa)

    real(RP), intent(inout) :: rhoq_t(ka,ia,ja,qa) ! tendency of rho * QTRC

    real(RP), intent(out) :: nu(ka,ia,ja) ! eddy viscosity (center)
    real(RP), intent(out) :: ri(ka,ia,ja) ! Richardson number
    real(RP), intent(out) :: pr(ka,ia,ja) ! Prantle number

    real(RP), intent(in)  :: momz(ka,ia,ja)
    real(RP), intent(in)  :: momx(ka,ia,ja)
    real(RP), intent(in)  :: momy(ka,ia,ja)
    real(RP), intent(in)  :: rhot(ka,ia,ja)
    real(RP), intent(in)  :: dens(ka,ia,ja)
    real(RP), intent(in)  :: qtrc(ka,ia,ja,qa)
    real(RP), intent(in)  :: n2(ka,ia,ja)

    real(RP), intent(in)  :: sflx_mw(ia,ja)
    real(RP), intent(in)  :: sflx_mu(ia,ja)
    real(RP), intent(in)  :: sflx_mv(ia,ja)
    real(RP), intent(in)  :: sflx_sh(ia,ja)
    real(RP), intent(in)  :: sflx_q (ia,ja,qa)

    real(RP), intent(in)  :: gsqrt   (ka,ia,ja,7)
    real(RP), intent(in)  :: j13g    (ka,ia,ja,7)
    real(RP), intent(in)  :: j23g    (ka,ia,ja,7)
    real(RP), intent(in)  :: j33g
    real(RP), intent(in)  :: mapf    (ia,ja,2,4)
    real(DP), intent(in)  :: dt

    real(RP) :: pott(ka,ia,ja)

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

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

    log_progress(*) 'atmosphere / physics / turbulence / DNS'

#ifdef DEBUG
    qflx_sgs_momz(:,:,:,:)   = undef
    qflx_sgs_momx(:,:,:,:)   = undef
    qflx_sgs_momy(:,:,:,:)   = undef
    qflx_sgs_rhot(:,:,:,:)   = undef
    qflx_sgs_rhoq(:,:,:,:,:) = undef

    pott(:,:,:) = undef
#endif

    nu(:,:,:) = 0.0_rp
    ri(:,:,:) = 0.0_rp
    pr(:,:,:) = 1.0_rp

    ! potential temperature
    do j = js-1, je+1
    do i = is-1, ie+1
    do k = ks, ke
       pott(k,i,j) = rhot(k,i,j) / dens(k,i,j)
    enddo
    enddo
    enddo

    !##### Start Upadate #####
    do jjs = js, je, jblock
    jje = jjs+jblock-1
    do iis = is, ie, iblock
    iie = iis+iblock-1

       !##### momentum equation (z) #####
       ! (cell center)
       do j = jjs, jje
       do i = iis, iie
       do k = ks+1, ke-1
          qflx_sgs_momz(k,i,j,zdir) = -atmos_phy_tb_dns_nu * ( momz(k,i,j)-momz(k-1,i,j) ) * rcdz(k)
       enddo
       enddo
       enddo

       do j = jjs, jje
       do i = iis, iie
          qflx_sgs_momz(ks,i,j,zdir) = 0.0_rp ! bottom boundary
          qflx_sgs_momz(ke,i,j,zdir) = 0.0_rp ! top boundary
       enddo
       enddo

       ! (y edge)
       do j = jjs,   jje
       do i = iis-1, iie
       do k = ks, ke-1
          qflx_sgs_momz(k,i,j,xdir) = -atmos_phy_tb_dns_nu * ( momz(k,i+1,j)-momz(k,i,j) ) * rfdx(i) * mapf(i,j,1,i_xy)
       enddo
       enddo
       enddo

       ! (x edge)
       do j = jjs-1, jje
       do i = iis,   iie
       do k = ks, ke-1
          qflx_sgs_momz(k,i,j,ydir) = -atmos_phy_tb_dns_nu * ( momz(k,i,j+1)-momz(k,i,j) ) * rfdy(j) * mapf(i,j,2,i_xy)
       enddo
       enddo
       enddo

       !##### momentum equation (x) #####
       ! (y edge)
       do j = jjs, jje
       do i = iis, iie
       do k = ks, ke-1
          qflx_sgs_momx(k,i,j,zdir) = -atmos_phy_tb_dns_nu * ( momx(k+1,i,j)-momx(k,i,j) ) * rfdz(k)
       enddo
       enddo
       enddo

       do j = jjs, jje
       do i = iis, iie
          qflx_sgs_momx(ks-1,i,j,zdir) = 0.0_rp ! bottom boundary
          qflx_sgs_momx(ke  ,i,j,zdir) = 0.0_rp ! top boundary
       enddo
       enddo

       ! (cell center)
       do j = jjs, jje
       do i = iis, iie+1
       do k = ks, ke
          qflx_sgs_momx(k,i,j,xdir) = -atmos_phy_tb_dns_nu * ( momx(k,i,j)-momx(k,i-1,j) ) * rcdx(i) * mapf(i,j,1,i_uy)
       enddo
       enddo
       enddo

       ! (z edge)
       do j = jjs-1, jje
       do i = iis,   iie
       do k = ks, ke
          qflx_sgs_momx(k,i,j,ydir) = -atmos_phy_tb_dns_nu * ( momx(k,i,j+1)-momx(k,i,j) ) * rfdy(j) * mapf(i,j,2,i_uy)
       enddo
       enddo
       enddo

       !##### momentum equation (y) #####

       ! (x edge)
       do j = jjs, jje
       do i = iis, iie
       do k = ks, ke-1
          qflx_sgs_momy(k,i,j,zdir) = -atmos_phy_tb_dns_nu * ( momy(k+1,i,j)-momy(k,i,j) ) * rfdz(k)
       enddo
       enddo
       enddo

       do j = jjs, jje
       do i = iis, iie
          qflx_sgs_momy(ks-1,i,j,zdir) = 0.0_rp ! bottom boundary
          qflx_sgs_momy(ke  ,i,j,zdir) = 0.0_rp ! top boundary
       enddo
       enddo

       ! (z edge)
       do j = jjs,   jje
       do i = iis-1, iie
       do k = ks, ke
          qflx_sgs_momy(k,i,j,xdir) = -atmos_phy_tb_dns_nu * ( momy(k,i+1,j)-momy(k,i,j) ) * rfdx(i) * mapf(i,j,1,i_xv)
       enddo
       enddo
       enddo

       ! (z-x plane)
       do j = jjs, jje+1
       do i = iis, iie
       do k = ks, ke
          qflx_sgs_momy(k,i,j,ydir) = -atmos_phy_tb_dns_nu * ( momy(k,i,j)-momy(k,i,j-1) ) * rcdy(j) * mapf(i,j,2,i_xv)
       enddo
       enddo
       enddo

       !##### Thermodynamic Equation #####

       ! at x, y ,w
       do j = jjs, jje
       do i = iis, iie
       do k = ks, ke-1
          qflx_sgs_rhot(k,i,j,zdir) = -0.5_rp * ( dens(k+1,i,j)+dens(k,i,j) ) &
                                    * atmos_phy_tb_dns_mu * ( pott(k+1,i,j)-pott(k,i,j) ) * rfdz(k)
       enddo
       enddo
       enddo

       do j = jjs, jje
       do i = iis, iie
          qflx_sgs_rhot(ks-1,i,j,zdir) = 0.0_rp
          qflx_sgs_rhot(ke  ,i,j,zdir) = 0.0_rp
       enddo
       enddo

       ! at u, y, z
       do j = jjs,   jje
       do i = iis-1, iie
       do k = ks, ke
          qflx_sgs_rhot(k,i,j,xdir) = -0.5_rp * ( dens(k,i+1,j)+dens(k,i,j) ) &
                                    * atmos_phy_tb_dns_mu * ( pott(k,i+1,j)-pott(k,i,j) ) * rfdx(i) * mapf(i,j,1,i_xy)
       enddo
       enddo
       enddo

       ! at x, v, z
       do j = jjs-1, jje
       do i = iis,   iie
       do k = ks, ke
          qflx_sgs_rhot(k,i,j,ydir) = -0.5_rp * ( dens(k,i,j+1)+dens(k,i,j) ) &
                                    * atmos_phy_tb_dns_mu * ( pott(k,i,j+1)-pott(k,i,j) ) * rfdy(j) * mapf(i,j,2,i_xy)
       enddo
       enddo
       enddo

    enddo
    enddo

    !##### Tracers #####
    do iq = 1, qa

    if ( .not. tracer_advc(iq) ) cycle

    do jjs = js, je, jblock
    jje = jjs+jblock-1
    do iis = is, ie, iblock
    iie = iis+iblock-1

       ! at x, y ,w
       do j = jjs, jje
       do i = iis, iie
       do k = ks, ke-1
          qflx_sgs_rhoq(k,i,j,zdir,iq) = -0.5_rp * ( dens(k+1,i,j)+dens(k,i,j) ) &
                                       * atmos_phy_tb_dns_mu * ( qtrc(k+1,i,j,iq)-qtrc(k,i,j,iq) ) * rfdz(k)
       enddo
       enddo
       enddo
       do j = jjs, jje
       do i = iis, iie
          qflx_sgs_rhoq(ks-1,i,j,zdir,iq) = 0.0_rp
          qflx_sgs_rhoq(ke  ,i,j,zdir,iq) = 0.0_rp
       enddo
       enddo

       ! at u, y, z
       do j = jjs,   jje
       do i = iis-1, iie
       do k = ks,   ke
          qflx_sgs_rhoq(k,i,j,xdir,iq) = -0.5_rp * ( dens(k,i+1,j)+dens(k,i,j) ) &
                                       * atmos_phy_tb_dns_mu * ( qtrc(k,i+1,j,iq)-qtrc(k,i,j,iq) ) * rfdx(i) * mapf(i,j,1,i_xy)
       enddo
       enddo
       enddo

       ! at x, v, z
       do j = jjs-1, jje
       do i = iis,   iie
       do k = ks,   ke
          qflx_sgs_rhoq(k,i,j,ydir,iq) = -0.5_rp * ( dens(k,i,j+1)+dens(k,i,j) ) &
                                       * atmos_phy_tb_dns_mu * ( qtrc(k,i,j+1,iq)-qtrc(k,i,j,iq) ) * rfdy(j) * mapf(i,j,2,i_xy)
       enddo
       enddo
       enddo

    enddo
    enddo

    enddo ! scalar quantities loop

    return

Here is the caller graph for this function: