HDFSPArrayGeoField.h

// This file is part of the hdf4 data handler for the OPeNDAP data server.
// It retrieves the latitude and longitude fields for some special HDF4 data products. 
// The products include TRMML2_V6,TRMML3B_V6,CER_AVG,CER_ES4,CER_CDAY,CER_CGEO,CER_SRB,CER_SYN,CER_ZAVG,OBPGL2,OBPGL3
// To know more information about these products,check HDFSP.h.
// Each product stores lat/lon in different way, so we have to retrieve them differently. 
//  Authors:   Kent Yang <myang6@hdfgroup.org>  
// Copyright (c) The HDF Group
#ifndef HDFSPARRAYGeoField_H
#define HDFSPARRAYGeoField_H
 
#include <libdap/Array.h>
#include "hdf.h"
#include "mfhdf.h"
#include "HDFSPEnumType.h"
 
class HDFSPArrayGeoField:public libdap::Array
{
    public:
        HDFSPArrayGeoField (int32 rank, const std::string& filename, const int sdfd, int32 fieldref, int32 dtype, SPType sptype, int fieldtype, const std::string & fieldname, const std::string & n = "", libdap::BaseType * v = nullptr):
            libdap::Array (n, v),
            rank (rank),
            filename(filename),
            sdfd(sdfd),
            fieldref (fieldref),
            dtype (dtype),
            sptype (sptype), 
            fieldtype (fieldtype), 
            name (fieldname) {
        }
        ~ HDFSPArrayGeoField () override = default;
 
        // Standard way of DAP handlers to pass the coordinates of the subsetted region to the handlers
        // Return the number of elements to read. 
        int format_constraint (int *cor, int *step, int *edg);
 
        libdap::BaseType *ptr_duplicate () override
        {
            return new HDFSPArrayGeoField (*this);
        }
 
        bool read () override;
 
 
    private:

        int32 rank;

        std::string filename;
        int sdfd;

        int32 fieldref;

        int32 dtype;

        SPType sptype;

        int fieldtype;

        std::string name;
 
        // Read TRMM level 2 version 6 lat/lon
        void readtrmml3a_v6 (const int32 *, const int32 *, const int32 *, int);
 
        // Read TRMM level 2 version 6 lat/lon
        void readtrmml3c_v6 (const int32 *, const int32 *, const int32 *, int);
        // Read TRMM level 2 version 6 lat/lon
        void readtrmml2_v6 (const int32 *, const int32 *, const int32 *, int);
 
 
        // Read OBPG level 2 lat/lon
        void readobpgl2 (int32 *, int32 *, int32 *, int);
 
        // Read OBPG level 3 lat/lon
        void readobpgl3 (const int *, const int *,  int);
 
        // Read TRMM level 3 version 6 lat/lon
        void readtrmml3b_v6 (const int32 *, const int32 *, const int32 *, int);
 
        // Read TRMM level 3 version 7 lat/lon
        void readtrmml3_v7 (const int32 *,  const int32 *, int);
 
 
        // Read CERES SAVG and CERES ICCP_DAYLIKE lat/lon
        void readcersavgid1 (const int *, const int *, const int *, int);
 
        // Read CERES SAVG and ICCP_DAYLIKE lat/lon
        void readcersavgid2 (const int *, const int *, const int *, int);
 
        // Read CERES ZAVG lat/lon
        void readcerzavg (const int32 *, const int32 *, const int32 *, int);
 
        // Read CERES AVG and SYN lat/lon
        void readceravgsyn (int32 *, int32 *, int32 *, int);
 
        // Read CERES ES4 and ICCP_GEO lat/lon
        void readceres4ig (const int32 *, const int32 *, const int32 *, int);
 
        template <typename T>  void LatLon2DSubset (T* outlatlon, int xdim, T* latlon, const int32 * offset, const int32 * count, const int32 * step);
 
 
 
};
 
 
#endif