HurricaneTempForcing.H

Go to the documentation of this file.

#ifndef HURRICANE_TEMP_FORCING_H
#define HURRICANE_TEMP_FORCING_H
 
#include "src/core/FieldRepo.H"
#include "src/equation_systems/temperature/TemperatureSource.H"
#include "src/utilities/FieldPlaneAveraging.H"
#include "AMReX_REAL.H"
 
using namespace amrex::literals;
 
namespace kynema_sgf::pde::temperature {
 
class HurricaneTempForcing
    : public TemperatureSource::Register<HurricaneTempForcing>
{
 
public:
    static std::string identifier() { return "HurricaneTempForcing"; }
 
    explicit HurricaneTempForcing(const CFDSim& /*sim*/);
 
    ~HurricaneTempForcing() 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_dTdR{0.001_rt};

    amrex::Real m_dTzh{18000.};
};
 
} // namespace kynema_sgf::pde::temperature
 
#endif