module FElib / Fluid dyn solver / Atmosphere / Physics turbulence / Common More...

Functions/Subroutines
subroutine, public	atm_phy_tb_dgm_common_setup_variables (tends, auxvars, diagvars, tends_manager, auxvar_manager, diagvar_manager, tends_num_tot, mesh3d)

subroutine, public	atm_phy_tb_dgm_common_calc_lambda (lambda, cs, filter_fac, lmesh, elem, lmesh2d, elem2d)

subroutine, public	atm_phy_tb_dgm_common_cal_grad_qtrc (dfq1, dfq2, dfq3, drdx, drdy, drdz, kh, qtrc, ddens, dens_hyd, dx, dy, dz, sx, sy, sz, lift, lmesh, elem, lmesh2d, elem2d, is_bound, cal_grad_dens)
	Calculate parameterized diffusive mass flux of tracer with turbulent model.

subroutine, public	atm_phy_tb_dgm_common_cal_tend (momx_t, momy_t, momz_t, rhot_t, t11, t12, t13, t21, t22, t23, t31, t32, t33, df1, df2, df3, nu, kh, ddens_, momx_, momy_, momz_, drhot_, dens_hyd, pres_hyd, pres_, pt_, dx, dy, dz, sx, sy, sz, lift, lmesh, elem, lmesh2d, elem2d, is_bound)
	Calculate tendecies with turbulent model.

subroutine, public	atm_phy_tb_dgm_common_cal_tend_qtrc (rhoq_t, dfq1, dfq2, dfq3, kh, ddens_, dens_hyd, dx, dy, dz, sx, sy, sz, lift, lmesh, elem, lmesh2d, elem2d, is_bound)
	Calculate tendecies of tracer density with turbulent model.

Variables
integer, parameter, public	atmos_phy_tb_aux_t13_id = 1

integer, parameter, public	atmos_phy_tb_aux_t23_id = 2

integer, parameter, public	atmos_phy_tb_aux_t33_id = 3

integer, parameter, public	atmos_phy_tb_aux_diffflx3_id = 4

integer, parameter, public	atmos_phy_tb_aux_t31_id = 5

integer, parameter, public	atmos_phy_tb_aux_t32_id = 6

integer, parameter, public	atmos_phy_tb_aux_diffflx1_id = 7

integer, parameter, public	atmos_phy_tb_aux_diffflx2_id = 8

integer, parameter, public	atmos_phy_tb_aux_t11_id = 9

integer, parameter, public	atmos_phy_tb_aux_t12_id = 10

integer, parameter, public	atmos_phy_tb_aux_t21_id = 11

integer, parameter, public	atmos_phy_tb_aux_t22_id = 12

integer, parameter, public	atmos_phy_tb_aux_scalar_num = 4

integer, parameter, public	atmos_phy_tb_aux_hvec_num = 2

integer, parameter, public	atmos_phy_tb_aux_htensor_num = 1

integer, parameter, public	atmos_phy_tb_aux_num = 12

integer, parameter, public	atmos_phy_tb_diag_tke_id = 1

integer, parameter, public	atmos_phy_tb_diag_nu_id = 2

integer, parameter, public	atmos_phy_tb_diag_kh_id = 3

integer, parameter, public	atmos_phy_tb_diag_num = 3

integer, parameter, public	atmos_phy_tb_momx_t_id = 1

integer, parameter, public	atmos_phy_tb_momy_t_id = 2

integer, parameter, public	atmos_phy_tb_momz_t_id = 3

integer, parameter, public	atmos_phy_tb_rhot_t_id = 4

integer, parameter, public	atmos_phy_tb_tends_num1 = 4

Detailed Description

module FElib / Fluid dyn solver / Atmosphere / Physics turbulence / Common

Description: A coomon modules for atmospheric turbulent parameterization

Author: Yuta Kawai, Team SCALE

NAMELIST

No namelist group

History Output

name	description	unit	variable
TB_T13	stress tensor (T13)	kg.m-3.m2.s-2	TB_T13
TB_T23	stress tensor (T23)	kg.m-3.m2.s-2	TB_T23
TB_T33	stress tensor (T33)	kg.m-3.m2.s-2	TB_T33
DIFF_FLX3	diffusive heat flux (z) / density	m2/s.K/m	DIFF_FLX3
TB_T31	stress tensor (T31)	kg.m-3.m2.s-2	TB_T31
TB_T32	stress tensor (T32)	kg.m-3.m2.s-2	TB_T32
DIFF_FLX1	diffusive heat flux (x) / density	m2/s.K/m	DIFF_FLX1
DIFF_FLX2	diffusive heat flux (y) / density	m2/s.K/m	DIFF_FLX2
TB_T11	stress tensor (T11)	kg.m-3.m2.s-2	TB_T11
TB_T12	stress tensor (T12)	kg.m-3.m2.s-2	TB_T12
TB_T21	stress tensor (T21)	kg.m-3.m2.s-2	TB_T21
TB_T22	stress tensor (T22)	kg.m-3.m2.s-2	TB_T22
TKE	SGS turbluence kinetic energy	m2/s2	TKE
NU	eddy viscosity	m2/s	NU
KH	eddy diffusion	m2/s	KH
TB_MOMX_t	tendency of x-momentum in TB process	kg/m2/s2	TB_MOMX_t
TB_MOMY_t	tendency of y-momentum in TB process	kg/m2/s2	TB_MOMY_t
TB_MOMZ_t	tendency of z-momentum in TB process	kg/m2/s2	TB_MOMZ_t
TB_RHOT_t	tendency of rho*PT in TB process	kg/m3.K/s	TB_RHOT_t
TB_{TRACER_NAME}_t	tendency of {TRACER_NAME} in TB process; {TRACER_NAME} depends on the physics schemes, e.g., QV, QC, QR.	TRACER_UNIT	TRACER_NAME

Function/Subroutine Documentation

◆ atm_phy_tb_dgm_common_setup_variables()

subroutine, public scale_atm_phy_tb_dgm_common::atm_phy_tb_dgm_common_setup_variables	(	type(meshfield3d), dimension(tends_num_tot), intent(inout)	tends,
		type(meshfield3d), dimension(atmos_phy_tb_aux_num), intent(inout)	auxvars,
		type(meshfield3d), dimension(atmos_phy_tb_diag_num), intent(inout)	diagvars,
		type(modelvarmanager), intent(inout)	tends_manager,
		type(modelvarmanager), intent(inout)	auxvar_manager,
		type(modelvarmanager), intent(inout)	diagvar_manager,
		integer, intent(in)	tends_num_tot,
		class(meshbase3d), intent(in)	mesh3d )

Definition at line 93 of file scale_atm_phy_tb_dgm_common.F90.

    use scale_tracer, only: &
      qa, tracer_name, tracer_desc, tracer_unit    
    use scale_mesh_base3d, only: meshbase3d
    use scale_meshfield_base, only: meshfield3d
    implicit none
    integer, intent(in) :: TENDS_NUM_TOT
    type(MeshField3D), intent(inout) :: tends(TENDS_NUM_TOT)
    type(MeshField3D), intent(inout) :: auxvars(ATMOS_PHY_TB_AUX_NUM)
    type(MeshField3D), intent(inout) :: diagvars(ATMOS_PHY_TB_DIAG_NUM)
    type(ModelVarManager), intent(inout) :: tends_manager
    type(ModelVarManager), intent(inout) :: auxvar_manager
    type(ModelVarManager), intent(inout) :: diagvar_manager
    class(MeshBase3D), intent(in) :: mesh3D
 
    ! type(VariableInfo), intent(out) :: auxvar_info(ATMOS_PHY_TB_AUX_NUM)
    ! type(VariableInfo), intent(out) :: diagvar_info(ATMOS_PHY_TB_DIAG_NUM)
    ! type(VariableInfo), intent(out) :: tend_info(ATMOS_PHY_TB_TENDS_NUM1)
 
    type(VariableInfo) :: ATMOS_PHY_TB_AUX_VINFO(ATMOS_PHY_TB_AUX_NUM)
    DATA atmos_phy_tb_aux_vinfo / &
      variableinfo( atmos_phy_tb_aux_t13_id, 'TB_T13', 'stress tensor (T13)',                         & ! rho x nu x S : [kg/m3] x [m2/s] x [s-1]
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                ), &
      variableinfo( atmos_phy_tb_aux_t23_id, 'TB_T23', 'stress tensor (T23)',                         &
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                ), &
      variableinfo( atmos_phy_tb_aux_t33_id, 'TB_T33', 'stress tensor (T33)',                         &
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                ), &
      variableinfo( atmos_phy_tb_aux_diffflx3_id, 'DIFF_FLX3', 'diffusive heat flux (z) / density',   &
                    'm2/s.K/m',  3, 'XYZ',  ''                                                     ), &                                      
      variableinfo( atmos_phy_tb_aux_t31_id, 'TB_T31', 'stress tensor (T31)',                         &
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                ), &
      variableinfo( atmos_phy_tb_aux_t32_id, 'TB_T32', 'stress tensor (T32)',                         &
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                ), &
      variableinfo( atmos_phy_tb_aux_diffflx1_id, 'DIFF_FLX1', 'diffusive heat flux (x) / density',   &
                    'm2/s.K/m',  3, 'XYZ',  ''                                                     ), &
      variableinfo( atmos_phy_tb_aux_diffflx2_id, 'DIFF_FLX2', 'diffusive heat flux (y) / density',   &
                    'm2/s.K/m',  3, 'XYZ',  ''                                                     ), &
      variableinfo( atmos_phy_tb_aux_t11_id, 'TB_T11', 'stress tensor (T11)',                         &
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                ), &
      variableinfo( atmos_phy_tb_aux_t12_id, 'TB_T12', 'stress tensor (T12)',                         &
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                ), &
      variableinfo( atmos_phy_tb_aux_t21_id, 'TB_T21', 'stress tensor (T21)',                         &
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                ), &
      variableinfo( atmos_phy_tb_aux_t22_id, 'TB_T22', 'stress tensor (T22)',                         &
                    'kg.m-3.m2.s-2',  3, 'XYZ',  ''                                                )  /
 
    type(VariableInfo) :: ATMOS_PHY_TB_DIAG_VINFO(ATMOS_PHY_TB_DIAG_NUM)
    DATA atmos_phy_tb_diag_vinfo / &
      variableinfo( atmos_phy_tb_diag_tke_id, 'TKE', 'SGS turbluence kinetic energy',   &
                    'm2/s2',  3, 'XYZ',  ''                                          ), &
      variableinfo( atmos_phy_tb_diag_nu_id, 'NU', 'eddy viscosity',                    &
                    'm2/s',  3, 'XYZ',  ''                                           ), &
      variableinfo( atmos_phy_tb_diag_kh_id, 'KH', 'eddy diffusion',                    &
                    'm2/s',  3, 'XYZ',  ''                                           )  /
                    
    type(VariableInfo) :: ATMOS_PHY_TB_TEND_VINFO(ATMOS_PHY_TB_TENDS_NUM1)
    DATA atmos_phy_tb_tend_vinfo / &
      variableinfo( atmos_phy_tb_momx_t_id, 'TB_MOMX_t', 'tendency of x-momentum in TB process',    &
                    'kg/m2/s2',  3, 'XYZ',  ''                                                   ), &
      variableinfo( atmos_phy_tb_momy_t_id, 'TB_MOMY_t', 'tendency of y-momentum in TB process',    &
                    'kg/m2/s2',  3, 'XYZ',  ''                                                   ), &
      variableinfo( atmos_phy_tb_momz_t_id, 'TB_MOMZ_t', 'tendency of z-momentum in TB process',    &
                    'kg/m2/s2',  3, 'XYZ',  ''                                                   ), &
      variableinfo( atmos_phy_tb_rhot_t_id, 'TB_RHOT_t', 'tendency of rho*PT in TB process',        &
                    'kg/m3.K/s', 3, 'XYZ',  ''                                                   )  / 
 
    type(VariableInfo) :: qtrc_vinfo_tmp
 
    integer :: iv
    integer :: iq
    logical :: reg_file_hist
    !----------------------------------------------------------
 
    ! auxvar_info(:) = ATMOS_PHY_TB_AUX_VINFO
    ! diagvar_info(:) = ATMOS_PHY_TB_DIAG_VINFO
    ! tend_info(:) = ATMOS_PHY_TB_TEND_VINFO
 
    !- Initialize an object to manage tendencies with turbulent model
 
    reg_file_hist = .true.    
    do iv = 1, atmos_phy_tb_tends_num1
      call tends_manager%Regist( &
        atmos_phy_tb_tend_vinfo(iv), mesh3d,     & ! (in) 
        tends(iv), reg_file_hist,                & ! (out)
        fill_zero=.true.                         )
    end do
 
    qtrc_vinfo_tmp%ndims    = 3
    qtrc_vinfo_tmp%dim_type = 'XYZ'
    qtrc_vinfo_tmp%STDNAME  = ''
 
    do iq = 1, qa
      iv = atmos_phy_tb_tends_num1 + iq 
      qtrc_vinfo_tmp%keyID = iv
      qtrc_vinfo_tmp%NAME  = 'TB_'//trim(tracer_name(iq))//'_t'
      qtrc_vinfo_tmp%DESC  = 'tendency of '//trim(tracer_desc(iq))//' in TB process'
      qtrc_vinfo_tmp%UNIT  = trim(tracer_unit(iq))//'/s'
 
      call tends_manager%Regist( &
        qtrc_vinfo_tmp, mesh3d,                  & ! (in) 
        tends(iv), reg_file_hist,                & ! (out)
        fill_zero=.true.                         )
    end do
 
    !- Initialize an object to manage auxiliary variables with turbulent model
 
    reg_file_hist = .true.    
    do iv = 1, atmos_phy_tb_aux_num
      call auxvar_manager%Regist( &
        atmos_phy_tb_aux_vinfo(iv), mesh3d, & ! (in) 
        auxvars(iv), reg_file_hist,         & ! (out)
        fill_zero=.true.                    )
    end do
 
    !- Initialize an object to manage diagnostic variables with turbulent model
    
    reg_file_hist = .true.    
    do iv = 1, atmos_phy_tb_diag_num
      call diagvar_manager%Regist( &
        atmos_phy_tb_diag_vinfo(iv), mesh3d, & ! (in) 
        diagvars(iv), reg_file_hist,         & ! (out)
        fill_zero=.true.                     )
    end do
 
    return

Referenced by mod_atmos_phy_tb_vars::atmosphytbvars_init().

◆ atm_phy_tb_dgm_common_calc_lambda()

subroutine, public scale_atm_phy_tb_dgm_common::atm_phy_tb_dgm_common_calc_lambda	(	real(rp), dimension(elem%np,lmesh%ne), intent(out)	lambda,
		real(rp), intent(in)	cs,
		real(rp), intent(in)	filter_fac,
		class(localmesh3d), intent(in)	lmesh,
		class(elementbase3d), intent(in)	elem,
		class(localmesh2d), intent(in)	lmesh2d,
		class(elementbase2d), intent(in)	elem2d )

Definition at line 224 of file scale_atm_phy_tb_dgm_common.F90.

 
    use scale_const, only: &
      karman  => const_karman
    implicit none
 
    class(LocalMesh3D), intent(in) :: lmesh
    class(ElementBase3D), intent(in) :: elem
    class(LocalMesh2D), intent(in) :: lmesh2D
    class(ElementBase2D), intent(in) :: elem2D
    real(RP), intent(out) :: lambda(elem%Np,lmesh%Ne)
    real(RP), intent(in) :: Cs
    real(RP), intent(in) :: filter_fac
 
    integer :: ke
 
    real(RP) :: vol
    real(RP) :: he
    real(RP) :: lambda0
    real(RP) :: Zs(elem2D%Np)
    real(RP) :: Z1(elem2D%Np)
    real(RP) :: dz(elem  %Np)
    real(RP) :: FZ
    real(RP) :: elem_aspect_eff
 
    real(RP), parameter :: OneOverThree  = 1.0_rp / 3.0_rp
    !--------------------------------------------------------------------
 
    !$omp parallel do private( &
    !$omp vol, lambda0, Zs, Z1, dz, he, FZ, elem_aspect_eff )
    do ke=lmesh%NeS, lmesh%NeE
      vol = sum( elem%IntWeight_lgl(:) * lmesh%J(:,ke) * lmesh%Gsqrt(:,ke) )
      he = ( vol / ( dble(elem%PolyOrder_h+1)**2 * dble(elem%PolyOrder_v+1) ) )**oneoverthree
 
      fz = lmesh%zlev(elem%Colmask(elem%Nnode_v,1),ke) - lmesh%zlev(elem%Colmask(1,1),ke)      
      elem_aspect_eff = fact( fz/dble(elem%PolyOrder_v+1), sqrt(vol/fz)/dble(elem%PolyOrder_h+1), sqrt(vol/fz)/dble(elem%PolyOrder_h+1) )
      lambda0 = elem_aspect_eff * cs * filter_fac * he
 
      zs(:) = lmesh%zlev(elem%Hslice(:,1),lmesh%EMap3Dto2D(ke))
      z1(:) = lmesh%zlev(elem%Hslice(:,2),lmesh%EMap3Dto2D(ke))
      dz(:) = max( lmesh%zlev(:,ke) - zs(elem%IndexH2Dto3D(:)), he )
 
      !lambda(:,ke) = sqrt( 1.0_RP / (1.0_RP / lambda0**2 + 1.0_RP / ( KARMAN * max( dz(:), EPS ) )**2 ) )
      lambda(:,ke) = sqrt( 1.0_rp / (1.0_rp / lambda0**2 + 1.0_rp / ( karman * ( dz(:) + 1.0e-4_rp ) )**2 ) )    
    end do
 
    return

Referenced by scale_atm_phy_tb_dgm_globalsmg::atm_phy_tb_dgm_globalsmg_cal_grad(), and scale_atm_phy_tb_dgm_smg::atm_phy_tb_dgm_smg_cal_grad().

◆ atm_phy_tb_dgm_common_cal_grad_qtrc()

subroutine, public scale_atm_phy_tb_dgm_common::atm_phy_tb_dgm_common_cal_grad_qtrc	(	real(rp), dimension(elem%np,lmesh%nea), intent(out)	dfq1,
		real(rp), dimension(elem%np,lmesh%nea), intent(out)	dfq2,
		real(rp), dimension(elem%np,lmesh%nea), intent(out)	dfq3,
		real(rp), dimension(elem%np,lmesh%nea), intent(inout)	drdx,
		real(rp), dimension(elem%np,lmesh%nea), intent(inout)	drdy,
		real(rp), dimension(elem%np,lmesh%nea), intent(inout)	drdz,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	kh,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	qtrc,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	ddens,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	dens_hyd,
		type(sparsemat), intent(in)	dx,
		type(sparsemat), intent(in)	dy,
		type(sparsemat), intent(in)	dz,
		type(sparsemat), intent(in)	sx,
		type(sparsemat), intent(in)	sy,
		type(sparsemat), intent(in)	sz,
		type(sparsemat), intent(in)	lift,
		class(localmesh3d), intent(in)	lmesh,
		class(elementbase3d), intent(in)	elem,
		class(localmesh2d), intent(in)	lmesh2d,
		class(elementbase2d), intent(in)	elem2d,
		logical, dimension(elem%nfptot,lmesh%ne), intent(in)	is_bound,
		logical, intent(in)	cal_grad_dens )

Calculate parameterized diffusive mass flux of tracer with turbulent model.

Parameters

[out]	dfq1	Diffusive mass flux in x1 direction / density (Kh dq/dx1)
[out]	dfq2	Diffusive mass flux in x2 direction / density (Kh dq/dx2)
[out]	dfq3	Diffusive mass flux in x3 direction / density (Kh dq/dx3)
[in,out]	drdx	Spatial gradient of density in x1 direction
[in,out]	drdy	Spatial gradient of density in x2 direction
[in,out]	drdz	Spatial gradient of density in x3 direction
[in]	kh	Eddy diffusivity
[in]	qtrc	Mass faction of tracer
[in]	ddens	Density perturbation
[in]	dens_hyd	Reference desity in hydrostatic state
[in]	dz	Differential matrix managed by sparse matrix type
[in]	sz	Stiffness matrix managed by sparse matrix type
[in]	lift	Lifting matrix managed by sparse matrix type
[in]	is_bound	Flag whether nodes are located at domain boundaries
[in]	cal_grad_dens	Flag whether spatial gradients of density are calcuated

Definition at line 278 of file scale_atm_phy_tb_dgm_common.F90.

 
    implicit none
 
    class(LocalMesh3D), intent(in) :: lmesh
    class(ElementBase3D), intent(in) :: elem
    class(LocalMesh2D), intent(in) :: lmesh2D
    class(ElementBase2D), intent(in) :: elem2D
    real(RP), intent(out) :: DFQ1(elem%Np,lmesh%NeA)
    real(RP), intent(out) :: DFQ2(elem%Np,lmesh%NeA)
    real(RP), intent(out) :: DFQ3(elem%Np,lmesh%NeA)
    real(RP), intent(inout) :: dRdx(elem%Np,lmesh%NeA)
    real(RP), intent(inout) :: dRdy(elem%Np,lmesh%NeA)
    real(RP), intent(inout) :: dRdz(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: Kh(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: QTRC(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DDENS(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DENS_hyd(elem%Np,lmesh%NeA)
    type(SparseMat), intent(in) :: Dx, Dy, Dz
    type(SparseMat), intent(in) :: Sx, Sy, Sz
    type(SparseMat), intent(in) :: Lift
    logical, intent(in) :: is_bound(elem%NfpTot,lmesh%Ne)
    logical, intent(in) :: cal_grad_dens
 
    real(RP) :: Fx(elem%Np), Fy(elem%Np), Fz(elem%Np), LiftDelFlx(elem%Np)
    real(RP) :: bnd_flux(elem%NfpTot,lmesh%Ne,3)
    real(RP) :: bnd_flux_rho (elem%NfpTot,lmesh%Ne,3)    
    
    real(RP) :: DENS(elem%Np), RDENS(elem%Np)
    real(RP) :: RHOxQTRC(elem%Np)
 
    integer :: ke
    integer :: p
    !--------------------------------------------------------------------
 
    call cal_bnd_flux_grad_qtrc( bnd_flux, bnd_flux_rho,                      & ! (out)
      qtrc, ddens, dens_hyd,                                                  & ! (in)
      lmesh%normal_fn(:,:,1), lmesh%normal_fn(:,:,2), lmesh%normal_fn(:,:,3), & ! (in)
      lmesh%vmapM, lmesh%vmapP,                                               & ! (in)
      lmesh, elem, is_bound, cal_grad_dens )                                    ! (in)
 
    !$omp parallel private( ke, Fx, Fy, Fz, LiftDelFlx, RHOxQTRC, DENS, RDENS )
 
    ! Calculate gradient of density
    if ( cal_grad_dens ) then
      !$omp do
      do ke=lmesh%NeS, lmesh%NeE
        dens(:) = dens_hyd(:,ke) + ddens(:,ke)
 
        call sparsemat_matmul( dx, dens, fx )
        call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux_rho(:,ke,1), liftdelflx )
        drdx(:,ke) = lmesh%Escale(:,ke,1,1) * fx(:) + liftdelflx(:)
  
        call sparsemat_matmul( dy, dens, fy )
        call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux_rho(:,ke,2), liftdelflx )
        drdy(:,ke) = lmesh%Escale(:,ke,2,2) * fy(:) + liftdelflx(:)
  
        call sparsemat_matmul( dz, dens, fz )
        call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux_rho(:,ke,3), liftdelflx )
        drdz(:,ke) = lmesh%Escale(:,ke,3,3) * fz(:) + liftdelflx(:)        
      end do
    end if
 
    ! Calculate gradient of tracer
    !$omp do
    do ke=lmesh%NeS, lmesh%NeE
      dens(:) = dens_hyd(:,ke) + ddens(:,ke)
      rdens(:) = 1.0_rp / dens(:)
      rhoxqtrc(:) = dens(:) * qtrc(:,ke)
 
      !---
      call sparsemat_matmul( dx, rhoxqtrc(:), fx )
      call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux(:,ke,1), liftdelflx )
      dfq1(:,ke) = kh(:,ke) * ( lmesh%Escale(:,ke,1,1) * fx(:) + liftdelflx(:) - qtrc(:,ke) * drdx(:,ke) ) * rdens(:)
 
      call sparsemat_matmul( dy, rhoxqtrc(:), fy )
      call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux(:,ke,2), liftdelflx )
      dfq2(:,ke) = kh(:,ke) * ( lmesh%Escale(:,ke,2,2) * fy(:) + liftdelflx(:) - qtrc(:,ke) * drdy(:,ke) ) * rdens(:)
 
      call sparsemat_matmul( dz, rhoxqtrc(:), fz )
      call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux(:,ke,3), liftdelflx )
      dfq3(:,ke) = kh(:,ke) * ( lmesh%Escale(:,ke,3,3) * fz(:) + liftdelflx(:) - qtrc(:,ke) * drdz(:,ke) ) * rdens(:)
    end do
 
    !$omp end parallel
 
    return

◆ atm_phy_tb_dgm_common_cal_tend()

subroutine, public scale_atm_phy_tb_dgm_common::atm_phy_tb_dgm_common_cal_tend	(	real(rp), dimension(elem%np,lmesh%nea), intent(out)	momx_t,
		real(rp), dimension(elem%np,lmesh%nea), intent(out)	momy_t,
		real(rp), dimension(elem%np,lmesh%nea), intent(out)	momz_t,
		real(rp), dimension(elem%np,lmesh%nea), intent(out)	rhot_t,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t11,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t12,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t13,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t21,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t22,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t23,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t31,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t32,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	t33,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	df1,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	df2,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	df3,
		real(rp), dimension (elem%np,lmesh%nea), intent(in)	nu,
		real(rp), dimension (elem%np,lmesh%nea), intent(in)	kh,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	ddens_,
		real(rp), dimension (elem%np,lmesh%nea), intent(in)	momx_,
		real(rp), dimension (elem%np,lmesh%nea), intent(in)	momy_,
		real(rp), dimension (elem%np,lmesh%nea), intent(in)	momz_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	drhot_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	dens_hyd,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	pres_hyd,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	pres_,
		real(rp), dimension (elem%np,lmesh%nea), intent(in)	pt_,
		type(sparsemat), intent(in)	dx,
		type(sparsemat), intent(in)	dy,
		type(sparsemat), intent(in)	dz,
		type(sparsemat), intent(in)	sx,
		type(sparsemat), intent(in)	sy,
		type(sparsemat), intent(in)	sz,
		type(sparsemat), intent(in)	lift,
		class(localmesh3d), intent(in)	lmesh,
		class(elementbase3d), intent(in)	elem,
		class(localmesh2d), intent(in)	lmesh2d,
		class(elementbase2d), intent(in)	elem2d,
		logical, dimension(elem%nfptot,lmesh%ne), intent(in)	is_bound )

Calculate tendecies with turbulent model.

Parameters

[out]	momx_t	Tendency of momentum in x1 direction with turbulent model
[out]	momy_t	Tendency of momentum in x2 direction with turbulent model
[out]	momz_t	Tendency of momentum in x3 direction with turbulent model
[out]	rhot_t	Tendency of density x potential temperature with turbulent model
[in]	t11	(1,1) component of stress tensor
[in]	t12	(1,2) component of stress tensor
[in]	t13	(1,3) component of stress tensor
[in]	t21	(2,1) component of stress tensor
[in]	t22	(2,2) component of stress tensor
[in]	t23	(2,3) component of stress tensor
[in]	t31	(3,1) component of stress tensor
[in]	t32	(3,2) component of stress tensor
[in]	t33	(3,3) component of stress tensor
[in]	df1	Diffusive heat flux in x1 direction / density
[in]	df2	Diffusive heat flux in x2 direction / density
[in]	df3	Diffusive heat flux in x3 direction / density
[in]	nu	Eddy viscosity
[in]	kh	Eddy diffusivity
[in]	ddens_	Density perturbation
[in]	momx_	Momentum in x1 direction
[in]	momy_	Momentum in x2 direction
[in]	momz_	Momentum in x3 direction
[in]	drhot_	Density x potential temperature perturbation
[in]	dens_hyd	Reference pressure in hydrostatic balance
[in]	pres_hyd	Reference density in hydrostatic balance
[in]	pres_	Pressure
[in]	pt_	Potential temperature
[in]	dz	Differential matrix managed by sparse matrix type
[in]	sz	Stiffness matrix managed by sparse matrix type
[in]	lift	Lifting matrix managed by sparse matrix type
[in]	is_bound	Flag whether nodes are located at domain boundaries

Definition at line 375 of file scale_atm_phy_tb_dgm_common.F90.

 
    implicit none
 
    class(LocalMesh3D), intent(in) :: lmesh
    class(ElementBase3D), intent(in) :: elem
    class(LocalMesh2D), intent(in) :: lmesh2D
    class(ElementBase2D), intent(in) :: elem2D
    real(RP), intent(out) :: MOMX_t(elem%Np,lmesh%NeA)
    real(RP), intent(out) :: MOMY_t(elem%Np,lmesh%NeA)
    real(RP), intent(out) :: MOMZ_t(elem%Np,lmesh%NeA)
    real(RP), intent(out) :: RHOT_t(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T11(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T12(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T13(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T21(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T22(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T23(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T31(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T32(elem%Np,lmesh%NeA)
    real(RP), intent(in) :: T33(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DF1(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DF2(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DF3(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: Nu   (elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: Kh   (elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DDENS_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: MOMX_ (elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: MOMY_ (elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: MOMZ_ (elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DRHOT_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DENS_hyd(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: PRES_hyd(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: PRES_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: PT_  (elem%Np,lmesh%NeA)
    type(SparseMat), intent(in) :: Dx, Dy, Dz
    type(SparseMat), intent(in) :: Sx, Sy, Sz
    type(SparseMat), intent(in) :: Lift
    logical, intent(in) :: is_bound(elem%NfpTot,lmesh%Ne)
 
    integer :: ke
 
    real(RP) :: Fx(elem%Np), Fy(elem%Np), Fz(elem%Np), LiftDelFlx(elem%Np)
    real(RP) :: DENS(elem%Np)
    real(RP) :: bnd_flux_mom(elem%NfpTot,lmesh%Ne,3)
    real(RP) :: bnd_flux_rhot(elem%NfpTot,lmesh%Ne)    
    !--------------------------------------------------------------------
 
    call cal_bnd_flux( bnd_flux_mom, bnd_flux_rhot,                           & ! (out)
      t11, t12, t13, t21, t22, t23, t31, t32, t33,                            & ! (in)
      df1, df2, df3,                                                          & ! (in)
      nu, kh,                                                                 & ! (in)
      ddens_, momx_, momy_, momz_, drhot_, dens_hyd, pres_hyd,                & ! (in)
      lmesh%normal_fn(:,:,1), lmesh%normal_fn(:,:,2), lmesh%normal_fn(:,:,3), & ! (in)
      lmesh%vmapM, lmesh%vmapP, lmesh, elem, is_bound   )                       ! (in)
 
    !$omp parallel do private( &
    !$omp Fx, Fy, Fz, LiftDelFlx,                    &
    !$omp DENS                                       )
    do ke=lmesh%NeS, lmesh%NeE
      dens(:) = dens_hyd(:,ke) + ddens_(:,ke)
 
      ! MOMX
      call sparsemat_matmul( dx, t11(:,ke), fx )
      call sparsemat_matmul( dy, t12(:,ke), fy )
      call sparsemat_matmul( dz, t13(:,ke), fz )
      call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux_mom(:,ke,1), liftdelflx )
 
      momx_t(:,ke) = lmesh%Escale(:,ke,1,1) * fx(:) + lmesh%Escale(:,ke,2,2) * fy(:) &
                   + lmesh%Escale(:,ke,3,3) * fz(:) + liftdelflx(:)
 
      ! MOMY
      call sparsemat_matmul( dx, t21(:,ke), fx )
      call sparsemat_matmul( dy, t22(:,ke), fy )
      call sparsemat_matmul( dz, t23(:,ke), fz )
      call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux_mom(:,ke,2), liftdelflx )
 
      momy_t(:,ke) = lmesh%Escale(:,ke,1,1) * fx(:) + lmesh%Escale(:,ke,2,2) * fy(:) &
                   + lmesh%Escale(:,ke,3,3) * fz(:) + liftdelflx(:)
 
      ! MOMZ
      call sparsemat_matmul( dx, t31(:,ke), fx )
      call sparsemat_matmul( dy, t32(:,ke), fy )
      call sparsemat_matmul( dz, t33(:,ke), fz )
      call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux_mom(:,ke,3), liftdelflx )
 
      momz_t(:,ke) = lmesh%Escale(:,ke,1,1) * fx(:) + lmesh%Escale(:,ke,2,2) * fy(:) &
                   + lmesh%Escale(:,ke,3,3) * fz(:) + liftdelflx(:)
 
      ! RHOT
      call sparsemat_matmul( dx, dens(:) * df1(:,ke), fx )
      call sparsemat_matmul( dy, dens(:) * df2(:,ke), fy )
      call sparsemat_matmul( dz, dens(:) * df3(:,ke), fz )
      call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * bnd_flux_rhot(:,ke), liftdelflx )
 
      rhot_t(:,ke)  = ( lmesh%Escale(:,ke,1,1) * fx(:) &
                      + lmesh%Escale(:,ke,2,2) * fy(:) &
                      + lmesh%Escale(:,ke,3,3) * fz(:) &
                      + liftdelflx(:) ) / lmesh%Gsqrt(:,ke)
 
    end do
 
    return

◆ atm_phy_tb_dgm_common_cal_tend_qtrc()

subroutine, public scale_atm_phy_tb_dgm_common::atm_phy_tb_dgm_common_cal_tend_qtrc	(	real(rp), dimension(elem%np,lmesh%nea), intent(out)	rhoq_t,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	dfq1,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	dfq2,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	dfq3,
		real(rp), dimension (elem%np,lmesh%nea), intent(in)	kh,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	ddens_,
		real(rp), dimension(elem%np,lmesh%nea), intent(in)	dens_hyd,
		type(sparsemat), intent(in)	dx,
		type(sparsemat), intent(in)	dy,
		type(sparsemat), intent(in)	dz,
		type(sparsemat), intent(in)	sx,
		type(sparsemat), intent(in)	sy,
		type(sparsemat), intent(in)	sz,
		type(sparsemat), intent(in)	lift,
		class(localmesh3d), intent(in)	lmesh,
		class(elementbase3d), intent(in)	elem,
		class(localmesh2d), intent(in)	lmesh2d,
		class(elementbase2d), intent(in)	elem2d,
		logical, dimension(elem%nfptot,lmesh%ne), intent(in)	is_bound )

Calculate tendecies of tracer density with turbulent model.

Parameters

[out]	rhoq_t	Tendency of tracer mass fraction
[in]	dfq1	Diffusive mass flux in x1 direction / density (Kh dq/dx1)
[in]	dfq2	Diffusive mass flux in x2 direction / density (Kh dq/dx2)
[in]	dfq3	Diffusive mass flux in x3 direction / density (Kh dq/dx3)
[in]	kh	Eddy diffusivity
[in]	ddens_	Density perturbation
[in]	dens_hyd	Reference pressure in hydrostatic state
[in]	dz	Differential matrix managed by sparse matrix type
[in]	sz	Stiffness matrix managed by sparse matrix type
[in]	lift	Lifting matrix managed by sparse matrix type
[in]	is_bound	Flag whether nodes are located at domain boundaries

Definition at line 491 of file scale_atm_phy_tb_dgm_common.F90.

 
    implicit none
 
    class(LocalMesh3D), intent(in) :: lmesh
    class(ElementBase3D), intent(in) :: elem
    class(LocalMesh2D), intent(in) :: lmesh2D
    class(ElementBase2D), intent(in) :: elem2D
    real(RP), intent(out) :: RHOQ_t(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DFQ1(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DFQ2(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DFQ3(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: Kh   (elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DDENS_(elem%Np,lmesh%NeA)
    real(RP), intent(in)  :: DENS_hyd(elem%Np,lmesh%NeA)
    type(SparseMat), intent(in) :: Dx, Dy, Dz
    type(SparseMat), intent(in) :: Sx, Sy, Sz
    type(SparseMat), intent(in) :: Lift
    logical, intent(in) :: is_bound(elem%NfpTot,lmesh%Ne)
 
    integer :: ke
 
    real(RP) :: Fx(elem%Np), Fy(elem%Np), Fz(elem%Np), LiftDelFlx(elem%Np)
    real(RP) :: DENS(elem%Np), RHOT(elem%Np)
    real(RP) :: del_flux(elem%NfpTot,lmesh%Ne)
    !--------------------------------------------------------------------
 
    call cal_bnd_flux_qtrc( del_flux,                                         & ! (out)
      dfq1, dfq2, dfq3,                                                       & ! (in)
      kh, ddens_, dens_hyd,                                                   & ! (in)
      lmesh%normal_fn(:,:,1), lmesh%normal_fn(:,:,2), lmesh%normal_fn(:,:,3), & ! (in)
      lmesh%vmapM, lmesh%vmapP,                                               & ! (in)
      lmesh, elem, is_bound )                                                   ! (in)
 
    !$omp parallel do private( &
    !$omp Fx, Fy, Fz, LiftDelFlx,            &
    !$omp DENS                               )
    do ke=lmesh%NeS, lmesh%NeE
      dens(:) = dens_hyd(:,ke) + ddens_(:,ke)
 
      ! RHOQ
      call sparsemat_matmul( dx, dens(:) * dfq1(:,ke), fx )
      call sparsemat_matmul( dy, dens(:) * dfq2(:,ke), fy )
      call sparsemat_matmul( dz, dens(:) * dfq3(:,ke), fz )
      call sparsemat_matmul( lift, lmesh%Fscale(:,ke) * del_flux(:,ke), liftdelflx )
 
      rhoq_t(:,ke) = lmesh%Escale(:,ke,1,1) * fx(:) + lmesh%Escale(:,ke,2,2) * fy(:) &
                   + lmesh%Escale(:,ke,3,3) * fz(:) + liftdelflx(:)
    end do
 
    return

Variable Documentation

◆ atmos_phy_tb_aux_t13_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t13_id = 1

Definition at line 50 of file scale_atm_phy_tb_dgm_common.F90.

50 integer, public, parameter :: ATMOS_PHY_TB_AUX_T13_ID = 1

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_t23_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t23_id = 2

Definition at line 51 of file scale_atm_phy_tb_dgm_common.F90.

51 integer, public, parameter :: ATMOS_PHY_TB_AUX_T23_ID = 2

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_t33_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t33_id = 3

Definition at line 52 of file scale_atm_phy_tb_dgm_common.F90.

52 integer, public, parameter :: ATMOS_PHY_TB_AUX_T33_ID = 3

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_diffflx3_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_diffflx3_id = 4

Definition at line 53 of file scale_atm_phy_tb_dgm_common.F90.

53 integer, public, parameter :: ATMOS_PHY_TB_AUX_DIFFFLX3_ID = 4

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_t31_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t31_id = 5

Definition at line 54 of file scale_atm_phy_tb_dgm_common.F90.

54 integer, public, parameter :: ATMOS_PHY_TB_AUX_T31_ID = 5

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_t32_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t32_id = 6

Definition at line 55 of file scale_atm_phy_tb_dgm_common.F90.

55 integer, public, parameter :: ATMOS_PHY_TB_AUX_T32_ID = 6

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_diffflx1_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_diffflx1_id = 7

Definition at line 56 of file scale_atm_phy_tb_dgm_common.F90.

56 integer, public, parameter :: ATMOS_PHY_TB_AUX_DIFFFLX1_ID = 7

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_diffflx2_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_diffflx2_id = 8

Definition at line 57 of file scale_atm_phy_tb_dgm_common.F90.

57 integer, public, parameter :: ATMOS_PHY_TB_AUX_DIFFFLX2_ID = 8

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_t11_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t11_id = 9

Definition at line 58 of file scale_atm_phy_tb_dgm_common.F90.

58 integer, public, parameter :: ATMOS_PHY_TB_AUX_T11_ID = 9

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_t12_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t12_id = 10

Definition at line 59 of file scale_atm_phy_tb_dgm_common.F90.

59 integer, public, parameter :: ATMOS_PHY_TB_AUX_T12_ID = 10

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_t21_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t21_id = 11

Definition at line 60 of file scale_atm_phy_tb_dgm_common.F90.

60 integer, public, parameter :: ATMOS_PHY_TB_AUX_T21_ID = 11

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_t22_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_t22_id = 12

Definition at line 61 of file scale_atm_phy_tb_dgm_common.F90.

61 integer, public, parameter :: ATMOS_PHY_TB_AUX_T22_ID = 12

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_aux_scalar_num

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_scalar_num = 4

Definition at line 62 of file scale_atm_phy_tb_dgm_common.F90.

62 integer, public, parameter :: ATMOS_PHY_TB_AUX_SCALAR_NUM = 4

Referenced by mod_atmos_phy_tb_vars::atmosphytbvars_init().

◆ atmos_phy_tb_aux_hvec_num

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_hvec_num = 2

Definition at line 63 of file scale_atm_phy_tb_dgm_common.F90.

63 integer, public, parameter :: ATMOS_PHY_TB_AUX_HVEC_NUM = 2

Referenced by mod_atmos_phy_tb_vars::atmosphytbvars_init().

◆ atmos_phy_tb_aux_htensor_num

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_htensor_num = 1

Definition at line 64 of file scale_atm_phy_tb_dgm_common.F90.

64 integer, public, parameter :: ATMOS_PHY_TB_AUX_HTENSOR_NUM = 1

Referenced by mod_atmos_phy_tb_vars::atmosphytbvars_init().

◆ atmos_phy_tb_aux_num

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_aux_num = 12

Definition at line 65 of file scale_atm_phy_tb_dgm_common.F90.

65 integer, public, parameter :: ATMOS_PHY_TB_AUX_NUM = 12

Referenced by atm_phy_tb_dgm_common_setup_variables(), and mod_atmos_phy_tb_vars::atmosphytbvars_init().

◆ atmos_phy_tb_diag_tke_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_diag_tke_id = 1

Definition at line 67 of file scale_atm_phy_tb_dgm_common.F90.

67 integer, public, parameter :: ATMOS_PHY_TB_DIAG_TKE_ID = 1

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_diag_nu_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_diag_nu_id = 2

Definition at line 68 of file scale_atm_phy_tb_dgm_common.F90.

68 integer, public, parameter :: ATMOS_PHY_TB_DIAG_NU_ID = 2

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_diag_kh_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_diag_kh_id = 3

Definition at line 69 of file scale_atm_phy_tb_dgm_common.F90.

69 integer, public, parameter :: ATMOS_PHY_TB_DIAG_KH_ID = 3

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_diag_num

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_diag_num = 3

Definition at line 70 of file scale_atm_phy_tb_dgm_common.F90.

70 integer, public, parameter :: ATMOS_PHY_TB_DIAG_NUM = 3

Referenced by atm_phy_tb_dgm_common_setup_variables(), and mod_atmos_phy_tb_vars::atmosphytbvars_init().

◆ atmos_phy_tb_momx_t_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_momx_t_id = 1

Definition at line 72 of file scale_atm_phy_tb_dgm_common.F90.

72 integer, public, parameter :: ATMOS_PHY_TB_MOMX_t_ID = 1

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_momy_t_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_momy_t_id = 2

Definition at line 73 of file scale_atm_phy_tb_dgm_common.F90.

73 integer, public, parameter :: ATMOS_PHY_TB_MOMY_t_ID = 2

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_momz_t_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_momz_t_id = 3

Definition at line 74 of file scale_atm_phy_tb_dgm_common.F90.

74 integer, public, parameter :: ATMOS_PHY_TB_MOMZ_t_ID = 3

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_rhot_t_id

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_rhot_t_id = 4

Definition at line 75 of file scale_atm_phy_tb_dgm_common.F90.

75 integer, public, parameter :: ATMOS_PHY_TB_RHOT_t_ID = 4

Referenced by atm_phy_tb_dgm_common_setup_variables().

◆ atmos_phy_tb_tends_num1

integer, parameter, public scale_atm_phy_tb_dgm_common::atmos_phy_tb_tends_num1 = 4

Definition at line 76 of file scale_atm_phy_tb_dgm_common.F90.

76 integer, public, parameter :: ATMOS_PHY_TB_TENDS_NUM1 = 4

Referenced by atm_phy_tb_dgm_common_setup_variables(), and mod_atmos_phy_tb_vars::atmosphytbvars_init().

Functions/Subroutines

Variables

Detailed Description

Function/Subroutine Documentation

◆ atm_phy_tb_dgm_common_setup_variables()

◆ atm_phy_tb_dgm_common_calc_lambda()

◆ atm_phy_tb_dgm_common_cal_grad_qtrc()

◆ atm_phy_tb_dgm_common_cal_tend()

◆ atm_phy_tb_dgm_common_cal_tend_qtrc()

Variable Documentation

◆ atmos_phy_tb_aux_t13_id

◆ atmos_phy_tb_aux_t23_id

◆ atmos_phy_tb_aux_t33_id

◆ atmos_phy_tb_aux_diffflx3_id

◆ atmos_phy_tb_aux_t31_id

◆ atmos_phy_tb_aux_t32_id

◆ atmos_phy_tb_aux_diffflx1_id

◆ atmos_phy_tb_aux_diffflx2_id

◆ atmos_phy_tb_aux_t11_id

◆ atmos_phy_tb_aux_t12_id

◆ atmos_phy_tb_aux_t21_id

◆ atmos_phy_tb_aux_t22_id

◆ atmos_phy_tb_aux_scalar_num

◆ atmos_phy_tb_aux_hvec_num

◆ atmos_phy_tb_aux_htensor_num

◆ atmos_phy_tb_aux_num

◆ atmos_phy_tb_diag_tke_id

◆ atmos_phy_tb_diag_nu_id

◆ atmos_phy_tb_diag_kh_id

◆ atmos_phy_tb_diag_num

◆ atmos_phy_tb_momx_t_id

◆ atmos_phy_tb_momy_t_id

◆ atmos_phy_tb_momz_t_id

◆ atmos_phy_tb_rhot_t_id

◆ atmos_phy_tb_tends_num1