NetCDFCFCollection.h

Go to the documentation of this file.

//
// $Id$
//
 
#ifndef _NetCDFCFCollection_h_
#define _NetCDFCFCollection_h_
 
#include <vector>
#include <map>
#include <algorithm>
 
#include <sstream>
#include <vapor/MyBase.h>
#include <vapor/UDUnitsClass.h>
#include <vapor/NetCDFCollection.h>
 
union ut_unit;
struct ut_system;
 
namespace VAPoR {
 
//
//
 
class VDF_API NetCDFCFCollection : public NetCDFCollection {
public:
    NetCDFCFCollection();
    virtual ~NetCDFCFCollection();
 
    virtual int Initialize(const std::vector<string> &files);
 
    //
    virtual bool IsCoordVarCF(string var) const { return (std::find(_coordinateVars.begin(), _coordinateVars.end(), var) != _coordinateVars.end()); }
 
    //
    virtual bool IsAuxCoordVarCF(string var) const { return (std::find(_auxCoordinateVars.begin(), _auxCoordinateVars.end(), var) != _auxCoordinateVars.end()); }
 
 
    //
    virtual bool IsCoordinateVar(string varName) const { return (IsCoordVarCF(varName) || IsAuxCoordVarCF(varName)); }
 
 
 
    //
    virtual bool IsLatCoordVar(string var) const { return (std::find(_latCoordVars.begin(), _latCoordVars.end(), var) != _latCoordVars.end()); }
 
    //
    virtual bool IsLonCoordVar(string var) const { return (std::find(_lonCoordVars.begin(), _lonCoordVars.end(), var) != _lonCoordVars.end()); }
 
    //
    virtual bool IsTimeCoordVar(string var) const { return (std::find(_timeCoordVars.begin(), _timeCoordVars.end(), var) != _timeCoordVars.end()); }
 
    //
    virtual bool IsVertCoordVar(string var) const { return (std::find(_vertCoordVars.begin(), _vertCoordVars.end(), var) != _vertCoordVars.end()); }
 
    //
    virtual bool IsVertCoordVarPressure(string var) const;
 
    //
    virtual bool IsVertCoordVarLength(string var) const;
 
    virtual bool IsVertCoordVarUp(string var) const;
 
    virtual std::vector<string> GetLatCoordVars() const { return (_latCoordVars); };
 
    virtual std::vector<string> GetLonCoordVars() const { return (_lonCoordVars); };
 
    virtual std::vector<string> GetTimeCoordVars() const { return (_timeCoordVars); };
 
    virtual std::vector<string> GetVertCoordVars() const { return (_vertCoordVars); };
 
    //
    virtual std::vector<string> GetDataVariableNames(int ndim, bool spatial) const;
 
    //
    virtual int GetVarCoordVarNames(string var, std::vector<string> &cvars) const;
 
    virtual int GetVarUnits(string var, string &units) const;
 
    const UDUnits *GetUDUnits() const { return (_udunit); };
 
    virtual int Convert(const string from, const string to, const double *src, double *dst, size_t n) const;
    virtual int Convert(const string from, const string to, const float *src, float *dst, size_t n) const;
 
    virtual bool GetMissingValue(string varname, double &mv) const;
 
    virtual int OpenRead(size_t ts, string varname);
 
    virtual bool IsVertDimensionless(string cvar) const;
 
    //
    bool GetMapProjectionProj4(string varname, string &proj4string) const;
 
    void FormatTimeStr(double time, string &str) const;
 
    friend std::ostream &operator<<(std::ostream &o, const NetCDFCFCollection &ncdfc);
 
private:
    std::vector<std::string>       _coordinateVars;
    std::vector<std::string>       _auxCoordinateVars;
    std::vector<std::string>       _lonCoordVars;
    std::vector<std::string>       _latCoordVars;
    std::vector<std::string>       _vertCoordVars;
    std::vector<std::string>       _timeCoordVars;
 
    UDUnits *_udunit;
 
    //
    // Map a variable name to it's missing value (if any)
    //
    std::map<string, double> _missingValueMap;
 
    std::vector<std::string> _GetCoordAttrs(const NetCDFSimple::Variable &varinfo) const;
 
    int _Initialize(const std::vector<string> &files);
 
    //
    bool _IsCoordinateVar(const NetCDFSimple::Variable &varinfo) const;
 
 
    bool _IsLonCoordVar(const NetCDFSimple::Variable &varinfo) const;
 
    bool _IsLatCoordVar(const NetCDFSimple::Variable &varinfo) const;
 
    bool _IsVertCoordVar(const NetCDFSimple::Variable &varinfo) const;
 
    bool _IsTimeCoordVar(const NetCDFSimple::Variable &varinfo) const;
 
    bool _GetMissingValue(string varname, string attname, double &mv) const;
    void _GetMissingValueMap(map<string, double> &missingValueMap) const;
 
};
};    // namespace VAPoR
 
#endif