ERF_ScalarLimitPost.H File Reference

This graph shows which files directly or indirectly include this file:

Go to the source code of this file.

Functions
	if (l_limit_Q)
	if (l_limit_KE||l_limit_Q)

Variables
const TurbChoice &	l_tc = solverChoice.turbChoice[level]
const bool	l_limit_KE = l_tc.use_tke
const bool	l_limit_Q = (solverChoice.moisture_type != MoistureType::None)
const bool	l_sbm = (solverChoice.moisture_type == MoistureType::SBM)
const bool	l_dirichlet_k
const int	l_klo = fine_geom.Domain().smallEnd(2)
const Real	l_tke_floor = l_tc.tke_floor
const Real	l_eps = std::numeric_limits<Real>::epsilon()
const int	l_nq = (l_sbm) ? 1 : S_new[IntVars::cons].nComp() - RhoQ1_comp

Function Documentation

if() [1/2]

if

(

l_limit_KE||

l_limit_Q

)

                                         {
#ifdef _OPENMP
#pragma omp parallel if (Gpu::notInLaunchRegion())
#endif
                for ( MFIter mfi(S_new[IntVars::cons],TilingIfNotGPU()); mfi.isValid(); ++mfi)
                {
                    const Box& tbx = mfi.tilebox();
                    const Array4<Real>& cons = S_new[IntVars::cons].array(mfi);
 
                    if (l_limit_KE) {
                        const bool skip_klo = l_dirichlet_k && (tbx.smallEnd(2) <= l_klo);
                        ParallelFor(tbx, [=] AMREX_GPU_DEVICE (int i, int j, int k) noexcept
                        {
                            if (skip_klo && k == l_klo) { return; }
                            const Real ke_floor = (l_tke_floor > zero) ? cons(i,j,k,Rho_comp) * l_tke_floor : l_eps;
                            cons(i,j,k,RhoKE_comp) = amrex::max(cons(i,j,k,RhoKE_comp), ke_floor);
                        });
                    }
 
                    if (l_limit_Q) {
                        ParallelFor(tbx, l_nq, [=,zero_d=zero] AMREX_GPU_DEVICE (int i, int j, int k, int nn) noexcept
                        {
                            const int n = RhoQ1_comp + nn;
                            cons(i,j,k,n) = amrex::max(cons(i,j,k,n), zero_d);
                        });
                    }
                } // mfi
            } // limit KE or Q

Here is the call graph for this function:

if() [2/2]

if

(

l_limit_Q

)

                           {
                erf_sbm::require_host_write_range_allowed(
                    l_sbm, RhoQ1_comp, l_nq, erf_sbm::HostWritePath::Positivity);
            }

Here is the call graph for this function:

Variable Documentation

l_dirichlet_k

const bool l_dirichlet_k

Initial value:

= ( l_tc.rans_type == RANSType::kEqn && l_tc.dirichlet_k &&
                                         (m_SurfaceLayer[Orientation(Direction::z, Orientation::low)] != nullptr) )

Referenced by erf_slow_rhs_post(), and if().

l_eps

const Real l_eps = std::numeric_limits<Real>::epsilon()

Referenced by if().

l_klo

const int l_klo = fine_geom.Domain().smallEnd(2)

Referenced by if().

l_limit_KE

const bool l_limit_KE = l_tc.use_tke

Referenced by if().

l_limit_Q

const bool l_limit_Q = (solverChoice.moisture_type != MoistureType::None)

Referenced by if().

l_nq

const int l_nq = (l_sbm) ? 1 : S_new[IntVars::cons].nComp() - RhoQ1_comp

Referenced by if().

l_sbm

const bool l_sbm = (solverChoice.moisture_type == MoistureType::SBM)

Referenced by if().

l_tc

const TurbChoice& l_tc = solverChoice.turbChoice[level]

l_tke_floor

const Real l_tke_floor = l_tc.tke_floor

Referenced by if().