EnergyBalanceModel.h

Go to the documentation of this file.

 
#ifndef ENERGYBALANCEMODEL
#define ENERGYBALANCEMODEL
 
#include "Context.h"
 
// -- constants -- //
 
constexpr float cp_air_mol = 29.25f;
constexpr float lambda_mol = 44000.f;
constexpr float von_Karman_constant = 0.41f;
constexpr float R = 8.314462618f;
 
inline float esat_Pa(float T_K) {
    float Tc = T_K - 273.15f;
    // This is Clausius-Clapeyron equation (See Campbell and Norman pp. 41 Eq. 3.8).  Note that temperature must be in Kelvin, and result is in Pascals
    return 611.0f * expf(17.502f * Tc / (Tc + 240.97f));
}
 
 
class EnergyBalanceModel {
public:
 
    EnergyBalanceModel(helios::Context *context);
 
 
    static int selfTest(int argc = 0, char **argv = nullptr);
 
    void enableMessages();
 
    void disableMessages();
 
    void run();
 
 
    void run(float dt);
 
 
    void run(const std::vector<uint> &UUIDs);
 
 
    void run(const std::vector<uint> &UUIDs, float dt);
 
 
    void addRadiationBand(const char *band);
 
 
    void addRadiationBand(const std::vector<std::string> &bands);
 
 
    void enableAirEnergyBalance();
 
 
    void enableAirEnergyBalance(float canopy_height_m, float reference_height_m);
 
 
    void evaluateAirEnergyBalance(float dt_sec, float time_advance_sec);
 
 
    void evaluateAirEnergyBalance(const std::vector<uint> &UUIDs, float dt_sec, float time_advance_sec);
 
 
    void optionalOutputPrimitiveData(const char *label);
 
#ifdef HELIOS_CUDA_AVAILABLE
 
    void enableGPUAcceleration();
 
 
    void disableGPUAcceleration();
 
 
    bool isGPUAccelerationEnabled() const;
#endif
 
    void printDefaultValueReport() const;
 
 
    void printDefaultValueReport(const std::vector<uint> &UUIDs) const;
 
private:
    void evaluateSurfaceEnergyBalance(const std::vector<uint> &UUIDs, float dt);
 
#ifdef HELIOS_CUDA_AVAILABLE
    void evaluateSurfaceEnergyBalance_GPU(const std::vector<uint> &UUIDs, float dt);
 
    void initializeGPUAcceleration();
#endif
 
    void evaluateSurfaceEnergyBalance_CPU(const std::vector<uint> &UUIDs, float dt);
 
    helios::Context *context;
 
#ifdef HELIOS_CUDA_AVAILABLE
    bool gpu_acceleration_enabled = false;
#endif
 
    float temperature_default;
 
    float wind_speed_default;
 
    float air_temperature_default;
 
    float air_humidity_default;
 
    float pressure_default;
 
    float gS_default;
 
    float heatcapacity_default;
 
    float Qother_default;
 
    float surface_humidity_default;
 
    bool air_energy_balance_enabled;
 
    helios::vec3 canopy_dimensions;
 
    float canopy_height_m;
    float reference_height_m;
 
    bool message_flag;
 
    std::vector<std::string> radiation_bands;
 
    std::vector<std::string> output_prim_data;
};
 
#endif