SolarPosition.h

Go to the documentation of this file.

 
#ifndef SOLARPOSITION
#define SOLARPOSITION
 
#include "Context.h"
 
class SolarPosition {
public:
 
    explicit SolarPosition(helios::Context *context_ptr);
 
 
    SolarPosition(float UTC_hrs, float latitude_deg, float longitude_deg, helios::Context *context_ptr);
 
    static int selfTest(int argc = 0, char **argv = nullptr);
 
    [[nodiscard]] helios::Time getSunriseTime() const;
 
    [[nodiscard]] helios::Time getSunsetTime() const;
 
    [[nodiscard]] float getSunElevation() const;
 
    [[nodiscard]] float getSunZenith() const;
 
    [[nodiscard]] float getSunAzimuth() const;
 
    [[nodiscard]] helios::vec3 getSunDirectionVector() const;
 
    [[nodiscard]] helios::SphericalCoord getSunDirectionSpherical() const;
 
 
    void setSunDirection(const helios::SphericalCoord &sundirection);
 
 
    [[nodiscard]] float getSolarFlux(float pressure_Pa, float temperature_K, float humidity_rel, float turbidity) const;
 
 
    [[nodiscard]] float getSolarFluxPAR(float pressure_Pa, float temperature_K, float humidity_rel, float turbidity) const;
 
 
    [[nodiscard]] float getSolarFluxNIR(float pressure_Pa, float temperature_K, float humidity_rel, float turbidity) const;
 
 
    [[nodiscard]] float getDiffuseFraction(float pressure_Pa, float temperature_K, float humidity_rel, float turbidity) const;
 
 
    [[nodiscard]] float getAmbientLongwaveFlux(float temperature_K, float humidity_rel) const;
 
 
    [[nodiscard]] float calibrateTurbidityFromTimeseries(const std::string &timeseries_shortwave_flux_label_Wm2) const;
 
 
    void enableCloudCalibration(const std::string &timeseries_shortwave_flux_label_Wm2);
 
    void disableCloudCalibration();
 
private:
    helios::Context *context;
 
    float UTC;
    float latitude;
    float longitude;
 
    bool issolarpositionoverridden = false;
    helios::SphericalCoord sun_direction;
 
    std::string cloudcalibrationlabel;
 
    [[nodiscard]] helios::SphericalCoord calculateSunDirection(const helios::Time &time, const helios::Date &date) const;
 
    void GueymardSolarModel(float pressure, float temperature, float humidity, float turbidity, float &Eb_PAR, float &Eb_NIR, float &fdiff) const;
 
    void applyCloudCalibration(float &R_calc_Wm2, float &fdiff_calc) const;
 
    static float turbidityResidualFunction(float turbidity, std::vector<float> &parameters, const void *a_solarpositionmodel);
};
 
#endif