GeostrophicForcing.H

Go to the documentation of this file.

#ifndef GEOSTROPHICFORCING_H
#define GEOSTROPHICFORCING_H
 
#include "src/equation_systems/icns/MomentumSource.H"
#include "src/core/SimTime.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf::pde::icns {

class GeostrophicForcing : public MomentumSource::Register<GeostrophicForcing>
{
public:
    static std::string identifier() { return "GeostrophicForcing"; }
 
    explicit GeostrophicForcing(const CFDSim& /*unused*/);
 
    ~GeostrophicForcing() override;
 
    void operator()(
        int lev,
        const amrex::MFIter& mfi,
        const amrex::Box& bx,
        FieldState fstate,
        const amrex::Array4<amrex::Real>& src_term) const override;
 
private:
    const SimTime& m_time;
    const amrex::AmrCore& m_mesh;

    std::string m_vel_timetable;

    amrex::Vector<amrex::Real> m_time_table;

    amrex::Vector<amrex::Real> m_speed_table;

    amrex::Vector<amrex::Real> m_direction_table;

    amrex::Real m_coriolis_factor;

    bool m_use_phase_ramp{false};

    int m_n_band{2};

    amrex::Vector<amrex::Real> m_target_vel{0.0_rt, 0.0_rt, 0.0_rt};

    amrex::Vector<amrex::Real> m_g_forcing{0.0_rt, 0.0_rt, 0.0_rt};
 
    bool m_is_horizontal{false};

    amrex::Real m_forcing_mphase0;
    amrex::Real m_forcing_mphase1;

    amrex::Real m_water_level;

    const Field* m_vof;
};
 
} // namespace kynema_sgf::pde::icns
 
#endif /* GEOSTROPHICFORCING_H */