ERF_InitCustomPert_MovingTerrain.H File Reference

Go to the source code of this file.

Functions
ParmParse	pp ("prob")
pp	query ("Ampl", Ampl)
pp	queryAdd ("wavelength", wavelength)
	AMREX_ALWAYS_ASSERT_WITH_MESSAGE (wavelength > zero, "prob.wavelength must be > 0")
	ParallelFor (bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const auto prob_lo=geomdata.ProbLo();const auto dx=geomdata.CellSize();const Real x=prob_lo[0]+(i+myhalf) *dx[0];Real z=z_cc(i, j, k);Real z_base=Ampl *std::sin(kp *x);z -=z_base;Real fac=std::cosh(kp *(z - H))/std::sinh(kp *H);Real p_prime=-(Ampl *omega *omega/kp) *fac *std::sin(kp *x) *r_hse(i, j, k);Real p_total=p_hse(i, j, k)+p_prime;state_pert(i, j, k, RhoTheta_comp)=getRhoThetagivenP(p_total) - getRhoThetagivenP(p_hse(i, j, k));state_pert(i, j, k, RhoScalar_comp)=zero;})

Variables
Real	Ampl = zero
Real	wavelength = Real(100.)
Real	H = phys_zhi(geomdata, lev)
Real	kp = two * PI / wavelength
Real	g = CONST_GRAV
Real	omega = std::sqrt(g * kp)

Function Documentation

AMREX_ALWAYS_ASSERT_WITH_MESSAGE()

AMREX_ALWAYS_ASSERT_WITH_MESSAGE	(	wavelength	,
		zero	,
		"prob.wavelength must	be,
		0"
	)

ParallelFor()

ParallelFor	(	bx	,
		[=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const auto prob_lo=geomdata.ProbLo();const auto dx=geomdata.CellSize();const Real x=prob_lo[0]+(i+myhalf) *dx[0];Real z=z_cc(i, j, k);Real z_base=Ampl *std::sin(kp *x);z -=z_base;Real fac=std::cosh(kp *(z - H))/std::sinh(kp *H);Real p_prime=-(Ampl *omega *omega/kp) *fac *std::sin(kp *x) *r_hse(i, j, k);Real p_total=p_hse(i, j, k)+p_prime;state_pert(i, j, k, RhoTheta_comp)=getRhoThetagivenP(p_total) - getRhoThetagivenP(p_hse(i, j, k));state_pert(i, j, k, RhoScalar_comp)=zero;}
	)

Here is the call graph for this function:

pp()

ParmParse pp

(

"prob"

)

query()

pp query	(	"Ampl"	,
		Ampl
	)

queryAdd()

pp queryAdd	(	"wavelength"	,
		wavelength
	)

Variable Documentation

Ampl

Real Ampl = zero

Referenced by init_my_custom_terrain(), and ParallelFor().

g

Real g = CONST_GRAV

H

Real H = phys_zhi(geomdata, lev)

Referenced by ERF::cf_flux_history(), ERF::cf_flux_history_record(), ERF::cf_flux_history_targets(), HSEutils::compute_F_for_temp_in_zone(), IBFaceSet::dump_faces(), erf_sbm::SpectralGrid::endpoints_to_two_moment(), ERF::FillPlot3DVars(), erf_sbm::inverse_two_moment(), make_mom_sources(), ParallelFor(), erf_sbm::try_endpoints_to_two_moment(), and erf_sbm::try_two_moment_to_endpoints().

kp

Real kp = two * PI / wavelength

Referenced by init_my_custom_terrain(), and ParallelFor().

omega

Real omega = std::sqrt(g * kp)

Referenced by ParallelFor().

wavelength

Real wavelength = Real(100.)

Referenced by rrtmgp::get_wavelength_index(), rrtmgp::get_wavelength_index_lw(), rrtmgp::get_wavelength_index_sw(), and init_my_custom_terrain().