metkitlib 1.16.0.9__cp312-cp312-manylinux_2_28_aarch64.whl

metkitlib/include/metkit/metkit_config.h

@@ -0,0 +1,26 @@
+/*
+ * (C) Copyright 1996- ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ * In applying this licence, ECMWF does not waive the privileges and immunities
+ * granted to it by virtue of its status as an intergovernmental organisation
+ * nor does it submit to any jurisdiction.
+ */
+
+#ifndef metkit_config_h
+#define metkit_config_h
+
+#include "metkit_ecbuild_config.h" /* generated by ecbuild */
+
+/* variables */
+
+#define metkit_HAVE_GRIB
+/* #undef metkit_HAVE_NETCDF */
+#define metkit_HAVE_BUFR
+/* #undef metkit_HAVE_ODB */
+/* #undef metkit_HAVE_FAIL_ON_CCSDS */
+
+/* packages */
+
+#endif /* metkit_config_h */

metkitlib/include/metkit/metkit_ecbuild_config.h

@@ -0,0 +1,72 @@
+/*
+ * (C) Copyright 2011- ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ * In applying this licence, ECMWF does not waive the privileges and immunities
+ * granted to it by virtue of its status as an intergovernmental organisation nor
+ * does it submit to any jurisdiction.
+ */
+
+#ifndef METKIT_ecbuild_config_h
+#define METKIT_ecbuild_config_h
+
+/* ecbuild info */
+
+#ifndef ECBUILD_VERSION_STR
+#define ECBUILD_VERSION_STR "3.13.0"
+#endif
+#ifndef ECBUILD_VERSION
+#define ECBUILD_VERSION "3.13.0"
+#endif
+#ifndef ECBUILD_MACROS_DIR
+#define ECBUILD_MACROS_DIR  "/src/ecbuild/cmake"
+#endif
+
+/* config info */
+
+#define METKIT_OS_NAME          "Linux-6.14.0-1014-azure"
+#define METKIT_OS_BITS          64
+#define METKIT_OS_BITS_STR      "64"
+#define METKIT_OS_STR           "linux.64"
+#define METKIT_OS_VERSION       "6.14.0-1014-azure"
+#define METKIT_SYS_PROCESSOR    "aarch64"
+
+#define METKIT_BUILD_TIMESTAMP  "20260126164436"
+#define METKIT_BUILD_TYPE       "RelWithDebInfo"
+
+#define METKIT_C_COMPILER_ID      "GNU"
+#define METKIT_C_COMPILER_VERSION "14.2.1"
+
+#define METKIT_CXX_COMPILER_ID      "GNU"
+#define METKIT_CXX_COMPILER_VERSION "14.2.1"
+
+#define METKIT_C_COMPILER       "/opt/rh/gcc-toolset-14/root/usr/bin/cc"
+#define METKIT_C_FLAGS          " -pipe -O2 -g -DNDEBUG"
+
+#define METKIT_CXX_COMPILER     "/opt/rh/gcc-toolset-14/root/usr/bin/c++"
+#define METKIT_CXX_FLAGS        " -pipe -Wall -Wextra -Wno-unused-parameter -Wno-unused-variable -Wno-sign-compare -O2 -g -DNDEBUG"
+
+/* Needed for finding per package config files */
+
+#define METKIT_INSTALL_DIR       "/tmp/metkit/target/metkit"
+#define METKIT_INSTALL_BIN_DIR   "/tmp/metkit/target/metkit/bin"
+#define METKIT_INSTALL_LIB_DIR   "/tmp/metkit/target/metkit/lib64"
+#define METKIT_INSTALL_DATA_DIR  "/tmp/metkit/target/metkit/share/metkit"
+
+#define METKIT_DEVELOPER_SRC_DIR "/src/metkit"
+#define METKIT_DEVELOPER_BIN_DIR "/tmp/metkit/build"
+
+/* Fortran support */
+
+#if 0
+
+#define METKIT_Fortran_COMPILER_ID      ""
+#define METKIT_Fortran_COMPILER_VERSION ""
+
+#define METKIT_Fortran_COMPILER ""
+#define METKIT_Fortran_FLAGS    ""
+
+#endif
+
+#endif /* METKIT_ecbuild_config_h */

metkitlib/include/metkit/metkit_version.h

@@ -0,0 +1,30 @@
+#ifndef metkit_version_h
+#define metkit_version_h
+
+#define metkit_VERSION_STR "1.16.0"
+#define metkit_VERSION     "1.16.0"
+
+#define metkit_VERSION_MAJOR 1
+#define metkit_VERSION_MINOR 16
+#define metkit_VERSION_PATCH 0
+
+#define metkit_GIT_SHA1 "f104efea40d66180a4c57f092ab5836d2f9dd084"
+
+#ifdef __cplusplus
+extern "C" {
+#endif
+
+const char * metkit_version();
+
+unsigned int metkit_version_int();
+
+const char * metkit_version_str();
+
+const char * metkit_git_sha1();
+
+#ifdef __cplusplus
+}
+#endif
+
+
+#endif // metkit_version_h

metkitlib/include/metkit/pointdb/DataSource.h

@@ -0,0 +1,91 @@
+/*
+ * (C) Copyright 1996-2013 ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ * In applying this licence, ECMWF does not waive the privileges and immunities
+ * granted to it by virtue of its status as an intergovernmental organisation nor
+ * does it submit to any jurisdiction.
+ */
+
+#ifndef metkit_DataSource_H
+#define metkit_DataSource_H
+
+#include <iosfwd>
+#include <map>
+
+#include "eckit/memory/NonCopyable.h"
+
+namespace eckit {
+class JSON;
+class Value;
+}  // namespace eckit
+
+namespace metkit {
+
+class MarsRequest;
+
+namespace pointdb {
+
+class DataSource;
+
+struct PointResult {
+
+    double value_;
+    double lat_;
+    double lon_;
+    const DataSource* source_;  // Warning, the source must have a longer life time as the result
+
+
+    void print(std::ostream& s) const;
+
+
+    friend std::ostream& operator<<(std::ostream& s, const PointResult& f) {
+        f.print(s);
+        return s;
+    }
+};
+
+
+class DataSourceHandler {
+public:
+
+    virtual void handle(DataSource*) = 0;
+};
+
+class DataSource : public eckit::NonCopyable {
+public:
+
+    virtual ~DataSource();
+
+
+    virtual PointResult extract(double lat, double lon) const = 0;
+
+    // Encode a MARS-like request representing the field
+    virtual const std::map<std::string, eckit::Value>& request() const = 0;
+
+    // A key to group source togther, e.g. sources poiting to the same file
+    virtual std::string groupKey() const = 0;
+
+    // A key to sort sources of the same group, e.g. offset in the file
+    virtual std::string sortKey() const = 0;
+
+    // Used to throw away requests in case of restarted transactions
+    virtual size_t batch() const;
+
+private:
+
+    virtual void print(std::ostream& s) const = 0;
+
+    friend std::ostream& operator<<(std::ostream& s, const DataSource& f) {
+        f.print(s);
+        return s;
+    }
+};
+
+
+}  // namespace pointdb
+}  // namespace metkit
+
+
+#endif

metkitlib/include/metkit/pointdb/FieldIndexer.h

@@ -0,0 +1,128 @@
+/*
+ * (C) Copyright 1996- ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ * In applying this licence, ECMWF does not waive the privileges and immunities
+ * granted to it by virtue of its status as an intergovernmental organisation nor
+ * does it submit to any jurisdiction.
+ */
+
+// File FieldIndex.h
+// Baudouin Raoult - ECMWF Oct 19
+
+#ifndef metkit_FieldIndex_H
+#define metkit_FieldIndex_H
+
+#include <iosfwd>
+#include <vector>
+
+
+namespace eckit {
+class Value;
+class JSON;
+}  // namespace eckit
+
+namespace metkit {
+
+class MarsRequest;
+
+namespace pointdb {
+
+class DataSource;
+class DataSourceHandler;
+//----------------------------------------------------------------------------------------------------------------------
+
+struct FieldIndexerStatus {
+    size_t count_;
+    size_t batch_;
+    FieldIndexerStatus(size_t count, size_t batch) : count_(count), batch_(batch) {}
+};
+
+
+class FieldIndexer {
+public:
+
+    // -- Exceptions
+    // None
+
+    // -- Contructors
+
+
+    // -- Destructor
+
+    virtual ~FieldIndexer();
+
+    // -- Convertors
+    // None
+
+    // -- Operators
+    // None
+
+    // -- Methods
+
+    virtual FieldIndexerStatus lookup(const eckit::Value&, DataSourceHandler&) const = 0;
+
+    // -- Overridden methods
+    // None
+
+    // -- Class members
+    // None
+
+    // -- Class methods
+    // None
+
+    // Uncomment for persistent, remove otherwise
+
+protected:
+
+    // -- Members
+    // None
+
+    // -- Methods
+
+    virtual void print(std::ostream& s) const = 0;
+
+    // -- Overridden methods
+    // None
+
+    // -- Class members
+    // None
+
+    // -- Class methods
+    // None
+
+private:
+
+    // No copy allowed
+
+    // -- Members
+
+
+    // -- Methods
+    // None
+
+    // -- Overridden methods
+    // None
+
+    // -- Class members
+    // None
+
+    // -- Class methods
+    // None
+
+    // -- Friends
+
+    friend std::ostream& operator<<(std::ostream& s, const FieldIndexer& p) {
+        p.print(s);
+        return s;
+    }
+};
+
+
+//----------------------------------------------------------------------------------------------------------------------
+}  // namespace pointdb
+
+}  // namespace metkit
+
+#endif

metkitlib/include/metkit/pointdb/GribDataSource.h

@@ -0,0 +1,47 @@
+/*
+ * (C) Copyright 1996-2013 ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ * In applying this licence, ECMWF does not waive the privileges and immunities
+ * granted to it by virtue of its status as an intergovernmental organisation nor
+ * does it submit to any jurisdiction.
+ */
+
+#ifndef metkit_GribDataSource_H
+#define metkit_GribDataSource_H
+
+#include "metkit/pointdb/DataSource.h"
+#include "metkit/pointdb/PointIndex.h"
+
+namespace eckit {
+class Offset;
+}
+
+namespace metkit {
+namespace pointdb {
+
+class GribFieldInfo;
+
+class GribDataSource : public DataSource {
+public:
+
+    virtual PointResult extract(double lat, double lon) const;
+
+private:
+
+    virtual double value(size_t index) const;
+    virtual std::string geographyHash() const;
+
+
+    virtual eckit::Offset seek(const eckit::Offset&) const = 0;
+    virtual long read(void*, long) const                   = 0;
+    virtual const GribFieldInfo& info() const              = 0;
+
+    friend class GribFieldInfo;
+};
+
+}  // namespace pointdb
+}  // namespace metkit
+
+#endif

metkitlib/include/metkit/pointdb/GribFieldInfo.h

@@ -0,0 +1,76 @@
+/*
+ * (C) Copyright 1996-2013 ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ * In applying this licence, ECMWF does not waive the privileges and immunities
+ * granted to it by virtue of its status as an intergovernmental organisation nor
+ * does it submit to any jurisdiction.
+ */
+
+#ifndef FieldInfoData_H
+#define FieldInfoData_H
+
+#include "eckit/io/Length.h"
+#include "eckit/io/Offset.h"
+#include "eckit/types/FixedString.h"
+
+
+namespace eckit {
+class PathName;
+}
+namespace metkit {
+namespace grib {
+class GribHandle;
+}
+
+
+namespace pointdb {
+
+
+class GribDataSource;
+
+
+class GribFieldInfo {
+public:
+
+    GribFieldInfo();
+
+    std::string geographyHash() const { return geographyHash_; }
+
+    bool ready() const { return numberOfValues_ > 0; }
+
+    void update(const grib::GribHandle& h);
+
+    double value(const GribDataSource&, size_t index) const;
+
+    bool useInterpolation() const { return sphericalHarmonics_ != 0; }
+    double interpolate(GribDataSource&, double& lat, double& lon) const;
+
+private:
+
+    double referenceValue_;
+    long binaryScaleFactor_;
+    long decimalScaleFactor_;
+    unsigned long bitsPerValue_;
+    unsigned long offsetBeforeData_;
+    unsigned long offsetBeforeBitmap_;
+    unsigned long numberOfValues_;
+    unsigned long numberOfDataPoints_;
+    long sphericalHarmonics_;
+
+    eckit::FixedString<32> geographyHash_;
+
+    void print(std::ostream&) const;
+
+    friend std::ostream& operator<<(std::ostream& s, const GribFieldInfo& f) {
+        f.print(s);
+        return s;
+    }
+};
+
+}  // namespace pointdb
+}  // namespace metkit
+
+
+#endif

metkitlib/include/metkit/pointdb/GribHandleDataSource.h

@@ -0,0 +1,58 @@
+/*
+ * (C) Copyright 1996-2013 ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ * In applying this licence, ECMWF does not waive the privileges and immunities
+ * granted to it by virtue of its status as an intergovernmental organisation nor
+ * does it submit to any jurisdiction.
+ */
+
+#ifndef metkit_GribHandleDataSource_H
+#define metkit_GribHandleDataSource_H
+
+#include "metkit/pointdb/GribDataSource.h"
+#include "metkit/pointdb/GribFieldInfo.h"
+
+namespace eckit {
+class DataHandle;
+}
+
+namespace metkit {
+namespace pointdb {
+
+
+class GribHandleDataSource : public GribDataSource {
+public:
+
+    GribHandleDataSource(const eckit::PathName&, const eckit::Offset& = 0);
+    GribHandleDataSource(eckit::DataHandle&, const eckit::Offset& = 0);
+    GribHandleDataSource(eckit::DataHandle*, const eckit::Offset& = 0);
+
+
+    ~GribHandleDataSource();
+
+private:
+
+    mutable eckit::DataHandle* handle_;
+    bool ownsHandle_;
+    mutable bool opened_;
+
+    mutable GribFieldInfo info_;
+    eckit::Offset offset_;
+
+    virtual eckit::Offset seek(const eckit::Offset&) const override;
+    virtual long read(void*, long) const override;
+    virtual const GribFieldInfo& info() const override;
+    virtual void print(std::ostream& s) const override;
+    virtual const std::map<std::string, eckit::Value>& request() const override;
+    virtual std::string groupKey() const override;
+    virtual std::string sortKey() const override;
+
+    void open() const;
+};
+
+}  // namespace pointdb
+}  // namespace metkit
+
+#endif

metkitlib/include/metkit/pointdb/PointIndex.h

@@ -0,0 +1,121 @@
+/*
+ * (C) Copyright 1996-2013 ECMWF.
+ *
+ * This software is licensed under the terms of the Apache Licence Version 2.0
+ * which can be obtained at http://www.apache.org/licenses/LICENSE-2.0.
+ * In applying this licence, ECMWF does not waive the privileges and immunities
+ * granted to it by virtue of its status as an intergovernmental organisation nor
+ * does it submit to any jurisdiction.
+ */
+
+#ifndef PointIndex_H
+#define PointIndex_H
+
+// #include <cmath>
+#include <memory>
+
+// #include "eckit/eckit.h"
+// #include "eckit/filesystem/PathName.h"
+
+
+#include "eckit/container/KDTree.h"
+#include "eckit/geometry/Point3.h"
+
+namespace metkit {
+namespace grib {
+class GribHandle;
+}
+
+namespace pointdb {
+
+
+struct LLPoint2 : public eckit::geometry::Point3 {
+
+    double lat_;
+    double lon_;
+    size_t payload_;
+
+
+    double lat() const { return lat_; }
+    double lon() const { return lon_; }
+
+
+    size_t payload() const { return payload_; }
+
+    const LLPoint2& point() const { return *this; }
+
+
+    LLPoint2() : eckit::geometry::Point3() {}
+
+    LLPoint2(double lat, double lon, size_t index) : eckit::geometry::Point3(), lat_(lat), lon_(lon) {
+
+        // See http://en.wikipedia.org/wiki/Geodetic_system#From_geodetic_to_ECEF
+        payload_  = index;
+        double& X = x_[0];
+        double& Y = x_[1];
+        double& Z = x_[2];
+
+        double h = 0;  // Altitude
+
+        const double earthRadius = 6378137.0;
+        double a                 = earthRadius;  // 6378137.0 ; //  WGS84 semi-major axis
+
+        double e2 = 0;  // 6.69437999014E-3; // WGS84 first numerical eccentricity sqared
+
+        double phi    = lat / 180.0 * M_PI;
+        double lambda = lon / 180.0 * M_PI;
+
+        double cos_phi    = cos(phi);
+        double sin_phi    = sin(phi);
+        double cos_lambda = cos(lambda);
+        double sin_lambda = sin(lambda);
+
+        double N_phi = a / sqrt(1 - e2 * sin_phi * sin_phi);
+
+        X = (N_phi + h) * cos_phi * cos_lambda;
+        Y = (N_phi + h) * cos_phi * sin_lambda;
+        Z = (N_phi * (1 - e2) + h) * sin_phi;
+    }
+
+    friend std::ostream& operator<<(std::ostream& s, const LLPoint2& p) {
+        s << '(' << p.lat_ << "," << p.lon_ << ' ' << p.payload_ << ')';
+        return s;
+    }
+};
+
+
+struct PointIndexTraits {
+    using Point   = LLPoint2;
+    using Payload = size_t;
+};
+
+class PointIndex {
+public:
+
+    typedef eckit::KDTreeMapped<PointIndexTraits> Tree;
+    typedef Tree::Point Point;
+    typedef Tree::NodeInfo NodeInfo;
+
+    NodeInfo nearestNeighbour(double lat, double lon);
+
+    static PointIndex& lookUp(const std::string& md5);
+    static std::string cache(const metkit::grib::GribHandle& h);
+
+    static eckit::PathName cachePath(const std::string& dir, const std::string& name);
+
+private:
+
+    PointIndex(const eckit::PathName&, Tree* tree = 0);
+    ~PointIndex();
+
+    eckit::PathName path_;
+    std::unique_ptr<Tree> tree_;
+
+    std::map<Point, NodeInfo> last_;
+    eckit::Mutex mutex_;
+};
+
+}  // namespace pointdb
+}  // namespace metkit
+
+#endif