SpaceWeather.hpp

浏览该文件的文档.

 
 
#pragma once
 
#include "AstGlobal.h"
#include "AstUtil/StringView.hpp"
#include "SpaceWeatherProvider.hpp"
#include <vector>
 
AST_NAMESPACE_BEGIN
 
class AST_CORE_API SpaceWeather: public SpaceWeatherProvider
{
public:
    struct Entry{
        int mjd;
        int BSRN;
        int ND;
        int Kp[8];
        int KpSum;
        int Ap[8];
        int ApAvg;
        double Cp;
        int C9;
        int ISN;
        int fluxQualifier;
        double F10p7Adj;
        double F10p7AdjCtr81;
        double F10p7AdjLst81;
        double F10p7Obs;
        double F10p7ObsCtr81;
        double F10p7ObsLst81;
    };
    SpaceWeather() = default;
    ~SpaceWeather() = default;
 
    errc_t loadDefault();
 
    errc_t load(StringView filepath);
 
    size_t size() const { return data_.size(); }
 
    const Entry* getEntry(int mjd) const;
 
    errc_t setEntry(int mjd, const Entry& entry);
 
public:
    double getApDaily(const TimePoint& tp) const override;
 
    double getApDaily_UTCMJD(double mjdUTC) const;
 
    int getAp3HourlyList(const TimePoint& tp, double* apList, int maxLookback) const override;
 
    double getKpDaily(const TimePoint& tp) const override;
 
    double getKpDaily_UTCMJD(double mjdUTC) const;
 
 
    double getF10p7Daily(const TimePoint& tp) const override;
 
    double getF10p7Daily_UTCMJD(double mjdUTC) const;
 
 
    double getF10p7Average(const TimePoint& tp) const override;
 
    double getF10p7Average_UTCMJD(double mjdUTC) const;
 
    static errc_t load(StringView filepath, std::vector<Entry>& data);
 
protected:
    void findEntryIndex(double mjdUTC, int& index, double& frac) const;
 
    void findFluxIndex(double mjdUTC, int& index, double& frac) const;
 
protected:
    std::vector<Entry> data_;
    int startMJD_ = 0;
    int endMJD_ = 0;
};
 
AST_NAMESPACE_END