derived Namespace Reference

Functions
void	erf_derrhodivide (const Box &bx, FArrayBox &derfab, const FArrayBox &datfab, const int scalar_index)
void	erf_dernull (const Box &, FArrayBox &, int, int, const FArrayBox &, const Geometry &, Real, const int *, const int)
void	erf_dersoundspeed (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const Geometry &, Real, const int *, const int)
void	erf_dertemp (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const Geometry &, Real, const int *, const int)
void	erf_dermoisttemp (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const Geometry &, Real, const int *, const int)
void	erf_dertheta (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const Geometry &, Real, const int *, const int)
void	erf_derscalar (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const Geometry &, Real, const int *, const int)
void	erf_derKE (const Box &bx, FArrayBox &derfab, int, int, const FArrayBox &datfab, const Geometry &, Real, const int *, const int)
void	erf_dervortx (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real, const int *, const int)
void	erf_dervorty (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real, const int *, const int)
void	erf_dervortz (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real, const int *, const int)
void	erf_derenstrophysq (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real, const int *, const int)
void	erf_dermagvel (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &, amrex::Real, const int *, const int)
void	erf_dermagvelsq (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &, amrex::Real, const int *, const int)
void	erf_derrhodivide (const amrex::Box &bx, amrex::FArrayBox &derfab, const amrex::FArrayBox &datfab, const int scalar_index)
void	erf_dernull (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
void	erf_dersoundspeed (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
void	erf_dertemp (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
void	erf_dermoisttemp (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
void	erf_dertheta (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
void	erf_derscalar (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)
void	erf_derKE (const amrex::Box &bx, amrex::FArrayBox &derfab, int dcomp, int ncomp, const amrex::FArrayBox &datfab, const amrex::Geometry &geomdata, amrex::Real time, const int *bcrec, const int level)

Function Documentation

erf_derenstrophysq()

void derived::erf_derenstrophysq	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	,
		const int *	,
		const int
	)

{
    AMREX_ALWAYS_ASSERT(dcomp == 0);
    AMREX_ALWAYS_ASSERT(ncomp == 1);
 
    auto const dat = datfab.array(); // cell-centered velocity
    auto tfab      = derfab.array(); // cell-centered vorticity x-component
 
    const Real dx = geomdata.CellSize(0);
    const Real dy = geomdata.CellSize(1);
    const Real dz = geomdata.CellSize(2);
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        Real vortx = (dat(i,j+1,k,2) - dat(i,j-1,k,2)) / (2.0*dy)  // dw/dy
                    -(dat(i,j,k+1,1) - dat(i,j,k-1,1)) / (2.0*dz); // dv/dz
        Real vorty = (dat(i,j,k+1,0) - dat(i,j,k-1,0)) / (2.0*dz)  // du/dz
                    -(dat(i+1,j,k,2) - dat(i-1,j,k,2)) / (2.0*dx); // dw/dx
        Real vortz = (dat(i+1,j,k,1) - dat(i-1,j,k,1)) / (2.0*dx)  // dv/dx
                    -(dat(i,j+1,k,0) - dat(i,j-1,k,0)) / (2.0*dy); // du/dy
 
        tfab(i,j,k,dcomp) = vortx*vortx + vorty*vorty + vortz*vortz;
    });
}

Referenced by ERF::sum_derived_quantities().

Here is the caller graph for this function:

erf_derKE() [1/2]

void derived::erf_derKE	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	time,
		const int *	bcrec,
		const int	level
	)

erf_derKE() [2/2]

void derived::erf_derKE	(	const Box &	bx,
		FArrayBox &	derfab,
		int	,
		int	,
		const FArrayBox &	datfab,
		const Geometry &	,
		Real	,
		const int *	,
		const int
	)

Function to define the kinetic energy KE by dividing (rho KE) by rho

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds KE @params[in] datfab array of data used to construct derived quantity

{
    erf_derrhodivide(bx, derfab, datfab, RhoKE_comp);
}

Referenced by ERF::Write3DPlotFile().

Here is the call graph for this function:

Here is the caller graph for this function:

erf_dermagvel()

void derived::erf_dermagvel	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	,
		amrex::Real	,
		const int *	,
		const int
	)

{
    AMREX_ALWAYS_ASSERT(dcomp == 0);
    AMREX_ALWAYS_ASSERT(ncomp == 1);
 
    auto const dat = datfab.array(); // cell-centered velocity
    auto tfab      = derfab.array(); // cell-centered magvel
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        Real u = dat(i,j,k,0);
        Real v = dat(i,j,k,1);
        Real w = dat(i,j,k,2);
        tfab(i,j,k,dcomp) = std::sqrt(u*u + v*v + w*w);
    });
}

Referenced by ERF::Write3DPlotFile().

Here is the caller graph for this function:

erf_dermagvelsq()

void derived::erf_dermagvelsq	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	,
		amrex::Real	,
		const int *	,
		const int
	)

{
    AMREX_ALWAYS_ASSERT(dcomp == 0);
    AMREX_ALWAYS_ASSERT(ncomp == 1);
 
    auto const dat = datfab.array(); // cell-centered velocity
    auto tfab      = derfab.array(); // cell-centered magvel
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        Real u = dat(i,j,k,0);
        Real v = dat(i,j,k,1);
        Real w = dat(i,j,k,2);
        tfab(i,j,k,dcomp) = u*u + v*v + w*w;
    });
}

Referenced by ERF::sum_derived_quantities().

Here is the caller graph for this function:

erf_dermoisttemp() [1/2]

void derived::erf_dermoisttemp	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	time,
		const int *	bcrec,
		const int	level
	)

erf_dermoisttemp() [2/2]

void derived::erf_dermoisttemp	(	const Box &	bx,
		FArrayBox &	derfab,
		int	,
		int	,
		const FArrayBox &	datfab,
		const Geometry &	,
		Real	,
		const int *	,
		const int
	)

{
    auto const dat = datfab.array();
    auto tfab      = derfab.array();
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        const Real rho = dat(i, j, k, Rho_comp);
        const Real rhotheta = dat(i, j, k, RhoTheta_comp);
        AMREX_ALWAYS_ASSERT(rhotheta > 0.);
        const Real qv = dat(i, j, k, RhoQ1_comp) / rho;
        tfab(i,j,k) = getTgivenRandRTh(rho,rhotheta,qv);
    });
}

Referenced by ERF::Write3DPlotFile(), and WriteBndryPlanes::write_planes().

Here is the call graph for this function:

Here is the caller graph for this function:

erf_dernull() [1/2]

void derived::erf_dernull	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	time,
		const int *	bcrec,
		const int	level
	)

erf_dernull() [2/2]

void derived::erf_dernull	(	const Box &	,
		FArrayBox &	,
		int	,
		int	,
		const FArrayBox &	,
		const Geometry &	,
		Real	,
		const int *	,
		const int
	)

Placeholder function that does nothing

{ }

Referenced by ERF::Write3DPlotFile().

Here is the caller graph for this function:

erf_derrhodivide() [1/2]

void derived::erf_derrhodivide	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		const amrex::FArrayBox &	datfab,
		const int	scalar_index
	)

erf_derrhodivide() [2/2]

void derived::erf_derrhodivide	(	const Box &	bx,
		FArrayBox &	derfab,
		const FArrayBox &	datfab,
		const int	scalar_index
	)

Function to define a derived quantity by dividing by density (analogous to cons_to_prim)

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity @params[in] datfab array of data used to construct derived quantity @params[in] scalar_index index of quantity to be divided by density

{
    // This routine divides any cell-centered conserved quantity by density
    auto const dat = datfab.array();
    auto primitive  = derfab.array();
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        const Real rho       = dat(i, j, k, Rho_comp);
        const Real conserved = dat(i, j, k, scalar_index);
        primitive(i,j,k) = conserved / rho;
    });
}

Referenced by erf_derKE(), erf_derscalar(), erf_dertheta(), and WriteBndryPlanes::write_planes().

Here is the caller graph for this function:

erf_derscalar() [1/2]

void derived::erf_derscalar	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	time,
		const int *	bcrec,
		const int	level
	)

erf_derscalar() [2/2]

void derived::erf_derscalar	(	const Box &	bx,
		FArrayBox &	derfab,
		int	,
		int	,
		const FArrayBox &	datfab,
		const Geometry &	,
		Real	,
		const int *	,
		const int
	)

Function to define a scalar s by dividing (rho s) by rho

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds scalar s @params[in] datfab array of data used to construct derived quantity

{
    erf_derrhodivide(bx, derfab, datfab, RhoScalar_comp);
}

Referenced by ERF::ErrorEst(), and ERF::Write3DPlotFile().

Here is the call graph for this function:

Here is the caller graph for this function:

erf_dersoundspeed() [1/2]

void derived::erf_dersoundspeed	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	time,
		const int *	bcrec,
		const int	level
	)

erf_dersoundspeed() [2/2]

void derived::erf_dersoundspeed	(	const Box &	bx,
		FArrayBox &	derfab,
		int	,
		int	,
		const FArrayBox &	datfab,
		const Geometry &	,
		Real	,
		const int *	,
		const int
	)

Function to define the sound speed by calling an EOS routine

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds pressure @params[in] datfab array of data used to construct derived quantity

{
    auto const dat = datfab.array();
    auto cfab      = derfab.array();
 
    // NOTE: we compute the soundspeed of dry air -- we do not account for any moisture effects here
    Real qv = 0.;
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        const Real rhotheta = dat(i, j, k, RhoTheta_comp);
        const Real rho      = dat(i, j, k, Rho_comp);
        AMREX_ALWAYS_ASSERT(rhotheta > 0.);
        cfab(i,j,k) = std::sqrt(Gamma * getPgivenRTh(rhotheta,qv) / rho);
    });
}

Referenced by ERF::Write3DPlotFile().

Here is the call graph for this function:

Here is the caller graph for this function:

erf_dertemp() [1/2]

void derived::erf_dertemp	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	time,
		const int *	bcrec,
		const int	level
	)

erf_dertemp() [2/2]

void derived::erf_dertemp	(	const Box &	bx,
		FArrayBox &	derfab,
		int	,
		int	,
		const FArrayBox &	datfab,
		const Geometry &	,
		Real	,
		const int *	,
		const int
	)

Function to define the temperature by calling an EOS routine

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds pressure @params[in] datfab array of data used to construct derived quantity

{
    auto const dat = datfab.array();
    auto tfab      = derfab.array();
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        const Real rho = dat(i, j, k, Rho_comp);
        const Real rhotheta = dat(i, j, k, RhoTheta_comp);
        AMREX_ALWAYS_ASSERT(rhotheta > 0.);
        tfab(i,j,k) = getTgivenRandRTh(rho,rhotheta);
    });
}

Referenced by ERF::Write3DPlotFile(), and WriteBndryPlanes::write_planes().

Here is the call graph for this function:

Here is the caller graph for this function:

erf_dertheta() [1/2]

void derived::erf_dertheta	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	time,
		const int *	bcrec,
		const int	level
	)

erf_dertheta() [2/2]

void derived::erf_dertheta	(	const Box &	bx,
		FArrayBox &	derfab,
		int	,
		int	,
		const FArrayBox &	datfab,
		const Geometry &	,
		Real	,
		const int *	,
		const int
	)

Function to define the potential temperature by calling an EOS routine

@params[in] bx box on which to divide by density @params[out] derfab array of derived quantity – here it holds pressure @params[in] datfab array of data used to construct derived quantity

{
    erf_derrhodivide(bx, derfab, datfab, RhoTheta_comp);
}

Referenced by ERF::ErrorEst(), and ERF::Write3DPlotFile().

Here is the call graph for this function:

Here is the caller graph for this function:

erf_dervortx()

void derived::erf_dervortx	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	,
		const int *	,
		const int
	)

{
    AMREX_ALWAYS_ASSERT(dcomp == 0);
    AMREX_ALWAYS_ASSERT(ncomp == 1);
 
    auto const dat = datfab.array(); // cell-centered velocity
    auto tfab      = derfab.array(); // cell-centered vorticity x-component
 
    const Real dy = geomdata.CellSize(1);
    const Real dz = geomdata.CellSize(2);
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        tfab(i,j,k,dcomp) = (dat(i,j+1,k,2) - dat(i,j-1,k,2)) / (2.0*dy)  // dw/dy
                          - (dat(i,j,k+1,1) - dat(i,j,k-1,1)) / (2.0*dz); // dv/dz
    });
}

Referenced by ERF::Write3DPlotFile().

Here is the caller graph for this function:

erf_dervorty()

void derived::erf_dervorty	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	,
		const int *	,
		const int
	)

{
    AMREX_ALWAYS_ASSERT(dcomp == 0);
    AMREX_ALWAYS_ASSERT(ncomp == 1);
 
    auto const dat = datfab.array(); // cell-centered velocity
    auto tfab      = derfab.array(); // cell-centered vorticity y-component
 
    const Real dx = geomdata.CellSize(0);
    const Real dz = geomdata.CellSize(2);
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        tfab(i,j,k,dcomp) = (dat(i,j,k+1,0) - dat(i,j,k-1,0)) / (2.0*dz)  // du/dz
                          - (dat(i+1,j,k,2) - dat(i-1,j,k,2)) / (2.0*dx); // dw/dx
    });
}

Referenced by ERF::Write3DPlotFile().

Here is the caller graph for this function:

erf_dervortz()

void derived::erf_dervortz	(	const amrex::Box &	bx,
		amrex::FArrayBox &	derfab,
		int	dcomp,
		int	ncomp,
		const amrex::FArrayBox &	datfab,
		const amrex::Geometry &	geomdata,
		amrex::Real	,
		const int *	,
		const int
	)

{
    AMREX_ALWAYS_ASSERT(dcomp == 0);
    AMREX_ALWAYS_ASSERT(ncomp == 1);
 
    auto const dat = datfab.array(); // cell-centered velocity
    auto tfab      = derfab.array(); // cell-centered vorticity z-component
 
    const Real dx = geomdata.CellSize(0);
    const Real dy = geomdata.CellSize(2);
 
    ParallelFor(bx, [=] AMREX_GPU_DEVICE(int i, int j, int k) noexcept
    {
        tfab(i,j,k,dcomp) = (dat(i+1,j,k,1) - dat(i-1,j,k,1)) / (2.0*dx)  // dv/dx
                          - (dat(i,j+1,k,0) - dat(i,j-1,k,0)) / (2.0*dy); // du/dy
    });
}

Referenced by ERF::ErrorEst(), and ERF::Write3DPlotFile().

Here is the caller graph for this function: