BoussinesqBubbleFieldInit.H

Go to the documentation of this file.

#ifndef BoussinesqBubbleFIELDINIT_H
#define BoussinesqBubbleFIELDINIT_H
 
#include "AMReX_Array.H"
#include "AMReX_Array4.H"
#include "AMReX_Box.H"
#include "AMReX_Geometry.H"
#include "AMReX_REAL.H"
#include "AMReX_Vector.H"
#include "AMReX_Gpu.H"
 
namespace amr_wind {
 
class BoussinesqBubbleFieldInit
{
    static_assert(
        AMREX_SPACEDIM == 3, "BoussinesqBubble requires 3 dimensional mesh");
 
public:
    BoussinesqBubbleFieldInit();
 
    void operator()(
        const amrex::Box& vbx,
        const amrex::Geometry& geom,
        const amrex::Array4<amrex::Real>& velocity,
        const amrex::Array4<amrex::Real>& density,
        const amrex::Array4<amrex::Real>& tracer) const;
 
private:
    amrex::Vector<amrex::Real> m_loc{0.25, 0.25, 0.5};
 
    amrex::Real m_rho{1.0};
 
    amrex::Real m_tracer_inner{0.0};
 
    amrex::Real m_tracer_outer{0.01};
 
    amrex::Real m_tracer_radius{0.1};
};
 
} // namespace amr_wind
 
#endif /* BoussinesqBubbleFIELDINIT_H */