ERF_InitCustomPert_KE.H File Reference

Go to the source code of this file.

Functions
ParmParse	pp ("prob")
pp	query ("KE_0", KE_0)
pp	query ("KE_decay_height", KE_decay_height)
pp	query ("KE_decay_order", KE_decay_order)
	if (KE_decay_height > 0 &&state_pert.nComp() > RhoKE_comp)
	if (state_pert.nComp() > RhoKE_comp)

Variables
Real	KE_0 = zero
Real	KE_decay_height = -1
Real	KE_decay_order = 1
const bool	use_eb = (sc.terrain_type == TerrainType::EB)
const int	klo = geomdata.Domain().smallEnd(2)

Function Documentation

if() [1/2]

if	(	KE_decay_height	,
		0 &&state_pert.	nComp(),
		RhoKE_comp
	)

                                         {
        amrex::Print() << "Initial KE profile (order " << KE_decay_order
                       << ") will extend up to " << KE_decay_height
                       << std::endl;
        AMREX_ALWAYS_ASSERT_WITH_MESSAGE(klo == bx.smallEnd(2),
        "Vertical decomposition with KE initialization is not allowed.");
    }

Here is the call graph for this function:

if() [2/2]

if	(	state_pert.	nComp(),
		RhoKE_comp
	)

                                         {
        AMREX_ALWAYS_ASSERT_WITH_MESSAGE(KE_0 >= zero,
            "KE_0 must be non-negative.");
 
        ParallelFor(bx, [=,one_d=one] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
        {
            // Set an initial value for SGS KE
            state_pert(i, j, k, RhoKE_comp) = r_hse(i,j,k) * KE_0;
            // scale initial SGS kinetic energy with height
            if (KE_decay_height > 0) {
                const Real z_sfc = (use_eb) ? zero :
                        fourth * ( z_nd(i  ,j  ,klo) + z_nd(i+1,j  ,klo)
                                 + z_nd(i  ,j+1,klo) + z_nd(i+1,j+1,klo) );
                const Real z_agl = z_cc(i,j,k) - z_sfc;
                state_pert(i, j, k, RhoKE_comp) *= max(
                    std::pow(1 - min(z_agl/KE_decay_height,one_d), KE_decay_order),
                    amrex::Real(1e-12));
            }
        });
    }

Here is the call graph for this function:

pp()

ParmParse pp

(

"prob"

)

query() [1/3]

pp query	(	"KE_0"	,
		KE_0
	)

query() [2/3]

pp query	(	"KE_decay_height"	,
		KE_decay_height
	)

query() [3/3]

pp query	(	"KE_decay_order"	,
		KE_decay_order
	)

Variable Documentation

KE_0

Real KE_0 = zero

Referenced by if().

KE_decay_height

Real KE_decay_height = -1

Referenced by if().

KE_decay_order

Real KE_decay_order = 1

Referenced by if().

klo

const int klo = geomdata.Domain().smallEnd(2)

Referenced by if().

use_eb

const bool use_eb = (sc.terrain_type == TerrainType::EB)

Referenced by if().