module ATMOSPHERE / Physics Turbulence More...

Functions/Subroutines
subroutine, public	atmos_phy_tb_dummy_setup (TYPE_TB, CDZ, CDX, CDY, CZ)

subroutine, public	atmos_phy_tb_dummy (qflx_sgs_momz, qflx_sgs_momx, qflx_sgs_momy, qflx_sgs_rhot, qflx_sgs_rhoq, tke, tke_t, nu_C, Ri, Pr, N2, MOMZ, MOMX, MOMY, RHOT, DENS, QTRC, sflx_mw, sflx_mu, sflx_mv, sflx_sh, sflx_qv, GSQRT, J13G, J23G, J33G, MAPF, dt)

Detailed Description

module ATMOSPHERE / Physics Turbulence

Description: dummy code

Author: Team SCALE

History

2013-01-22 (S.Nishizawa) [new]

Function/Subroutine Documentation

◆ atmos_phy_tb_dummy_setup()

subroutine, public scale_atmos_phy_tb_dummy::atmos_phy_tb_dummy_setup	(	character(len=*), intent(in)	TYPE_TB,
		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 53 of file scale_atmos_phy_tb_dummy.F90.

References scale_stdio::io_fid_log, scale_stdio::io_l, and scale_process::prc_mpistop().

Referenced by scale_atmos_phy_tb::atmos_phy_tb_setup().

     use scale_process, only: &
        prc_mpistop
     implicit none
 
     character(len=*), intent(in) :: type_tb
 
     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)
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*) '*** Turbulence Dummy'
 
     if ( type_tb /= 'DUMMY' ) then
        write(*,*) 'xxx ATMOS_PHY_TB_TYPE is not Dummy. Check!'
        call prc_mpistop
     endif
 
     return

Here is the call graph for this function:

Here is the caller graph for this function:

◆ atmos_phy_tb_dummy()

subroutine, public scale_atmos_phy_tb_dummy::atmos_phy_tb_dummy	(	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,qa,3), intent(out)	qflx_sgs_rhoq,
		real(rp), dimension(ka,ia,ja), intent(inout)	tke,
		real(rp), dimension(ka,ia,ja), intent(out)	tke_t,
		real(rp), dimension(ka,ia,ja), intent(out)	nu_C,
		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(out)	N2,
		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(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), intent(in)	sflx_qv,
		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 factorp

Definition at line 84 of file scale_atmos_phy_tb_dummy.F90.

References scale_stdio::io_fid_log, and scale_stdio::io_l.

Referenced by scale_atmos_phy_tb::atmos_phy_tb_setup().

     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,qa,3)
 
     real(RP), intent(inout) :: tke(ka,ia,ja)
     real(RP), intent(out)   :: tke_t(ka,ia,ja) ! tendency TKE
 
     real(RP), intent(out)   :: nu_c(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(out)   :: n2  (ka,ia,ja) ! squared Brunt-Vaisala frequency
 
     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)    :: 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_qv(ia,ja)
 
     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
     !---------------------------------------------------------------------------
 
     if( io_l ) write(io_fid_log,*) '*** Physics step: Turbulence(dummy)'
 
     ! do nothing
     qflx_sgs_momz(:,:,:,:)   = 0.0_rp
     qflx_sgs_momx(:,:,:,:)   = 0.0_rp
     qflx_sgs_momy(:,:,:,:)   = 0.0_rp
     qflx_sgs_rhot(:,:,:,:)   = 0.0_rp
     qflx_sgs_rhoq(:,:,:,:,:) = 0.0_rp
 
     tke_t(:,:,:) = 0.0_rp
 
     tke(:,:,:) = 0.0_rp
     nu_c(:,:,:) = 0.0_rp
     ri(:,:,:) = 0.0_rp
     pr(:,:,:) = 1.0_rp
     n2(:,:,:) = 0.0_rp
 
     return

Here is the caller graph for this function:

Functions/Subroutines

Detailed Description

Function/Subroutine Documentation

◆ atmos_phy_tb_dummy_setup()

◆ atmos_phy_tb_dummy()