ERF_InitCustomPertVels_MovingTerrain.H File Reference

Go to the source code of this file.

Functions
ParmParse	pp_prob ("prob")
pp_prob	query ("Ampl", Ampl)
pp_prob	queryAdd ("wavelength", wavelength)
	AMREX_ALWAYS_ASSERT_WITH_MESSAGE (wavelength > zero, "prob.wavelength must be > 0")
	ParallelFor (xbx, [=] 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 *dx[0];Real z=fourth *(z_nd(i, j, k)+z_nd(i, j+1, k)+z_nd(i, j, k+1)+z_nd(i, j+1, k+1));Real z_base=Ampl *std::sin(kp *x);z -=z_base;Real fac=std::cosh(kp *(z - H))/std::sinh(kp *H);x_vel_pert(i, j, k)=-Ampl *omega *fac *std::sin(kp *x);})
	ParallelFor (zbx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const auto prob_lo=geomdata.ProbLo();const auto dx=geomdata.CellSize();Real x=prob_lo[0]+(i+myhalf) *dx[0];Real z=fourth *(z_nd(i, j, k)+z_nd(i+1, j, k)+z_nd(i, j+1, k)+z_nd(i+1, j+1, k));Real z_base=Ampl *std::sin(kp *x);z -=z_base;Real fac=std::sinh(kp *(z - H))/std::sinh(kp *H);z_vel_pert(i, j, k)=Ampl *omega *fac *std::cos(kp *x);})

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() [1/2]

ParallelFor	(	xbx	,
		[=] 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 *dx[0];Real z=fourth *(z_nd(i, j, k)+z_nd(i, j+1, k)+z_nd(i, j, k+1)+z_nd(i, j+1, k+1));Real z_base=Ampl *std::sin(kp *x);z -=z_base;Real fac=std::cosh(kp *(z - H))/std::sinh(kp *H);x_vel_pert(i, j, k)=-Ampl *omega *fac *std::sin(kp *x);}
	)

ParallelFor() [2/2]

ParallelFor	(	zbx	,
		[=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept { const auto prob_lo=geomdata.ProbLo();const auto dx=geomdata.CellSize();Real x=prob_lo[0]+(i+myhalf) *dx[0];Real z=fourth *(z_nd(i, j, k)+z_nd(i+1, j, k)+z_nd(i, j+1, k)+z_nd(i+1, j+1, k));Real z_base=Ampl *std::sin(kp *x);z -=z_base;Real fac=std::sinh(kp *(z - H))/std::sinh(kp *H);z_vel_pert(i, j, k)=Ampl *omega *fac *std::cos(kp *x);}
	)

pp_prob()

ParmParse pp_prob

(

"prob"

)

query()

pp_prob query	(	"Ampl"	,
		Ampl
	)

queryAdd()

pp_prob queryAdd	(	"wavelength"	,
		wavelength
	)

Variable Documentation

Ampl

Real Ampl = zero

Referenced by ParallelFor().

g

Real g = CONST_GRAV

H

Real H = phys_zhi(geomdata, lev)

Referenced by ParallelFor().

kp

Real kp = two * PI / wavelength

Referenced by ParallelFor().

omega

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

Referenced by ParallelFor().

wavelength

Real wavelength = Real(100.)