HurricaneForcing.H

Go to the documentation of this file.

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

 
class HurricaneForcing : public MomentumSource::Register<HurricaneForcing>
{
public:
    static std::string identifier() { return "HurricaneForcing"; }
 
    explicit HurricaneForcing(const CFDSim& sim);
 
    ~HurricaneForcing() override;
 
    void operator()(
        int lev, FieldState fstate, amrex::MultiFab& src_term) const override;
 
    void mean_velocity_init(const VelPlaneAveraging& /*vavg*/);
 
    void mean_velocity_update(const VelPlaneAveraging& /*vavg*/);
 
private:
    const amrex::AmrCore& m_mesh;
 
    amrex::Gpu::DeviceVector<amrex::Real> m_vel_ht;
    amrex::Gpu::DeviceVector<amrex::Real> m_vel_vals;

    int m_axis{2};

    amrex::Real m_V{40.0_rt}; // V=40 m/s is the default value

    amrex::Real m_dVdR{-0.0008_rt};

    amrex::Real m_Vzh{18000.};

    amrex::Real m_R{40000.0_rt}; // 40km is taken as the default value

    amrex::Real m_coriolis_factor{0.0005_rt};
};
 
} // namespace kynema_sgf::pde::icns
 
#endif /* HURRICANE_FORCING_H */