CoriolisForcing.H

Go to the documentation of this file.

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

class CoriolisForcing : public MomentumSource::Register<CoriolisForcing>
{
public:
    static std::string identifier() { return "CoriolisForcing"; }
 
    explicit CoriolisForcing(const CFDSim& /*sim*/);
 
    ~CoriolisForcing() override;
 
    void operator()(
        int lev, FieldState fstate, amrex::MultiFab& src_term) const override;
 
private:
    const Field& m_velocity;


    amrex::Vector<amrex::Real> m_east{1.0_rt, 0.0_rt, 0.0_rt};
    amrex::Vector<amrex::Real> m_north{0.0_rt, 1.0_rt, 0.0_rt};
    amrex::Vector<amrex::Real> m_up{0.0_rt, 0.0_rt, 1.0_rt};

    amrex::Real m_latitude{0.0_rt};
 
    amrex::Real m_sinphi{0.0_rt};
    amrex::Real m_cosphi{0.0_rt};

    amrex::Real m_omega{7.3e-5_rt};

    amrex::Real m_coriolis_factor{0.0_rt};
 
    bool m_is_horizontal{false};
};
 
} // namespace kynema_sgf::pde::icns
 
#endif /* CORIOLISFORCING_H */