source: XIOS/trunk/src/io/nc4_data_input.hpp @ 1638

Last change on this file since 1638 was 1638, checked in by yushan, 5 years ago

dev on ADA
  • Property copyright set to
    Software name : XIOS (Xml I/O Server)
    http://forge.ipsl.jussieu.fr/ioserver
    Creation date : January 2009
    Licence : CeCCIL version2
    see license file in root directory : Licence_CeCILL_V2-en.txt
    or http://www.cecill.info/licences/Licence_CeCILL_V2-en.html
    Holder : CEA/LSCE (Laboratoire des Sciences du CLimat et de l'Environnement)
    CNRS/IPSL (Institut Pierre Simon Laplace)
    Project Manager : Yann Meurdesoif
    yann.meurdesoif@cea.fr
File size: 3.9 KB
Line 
1#ifndef __XIOS_NC4_DATA_INPUT__
2#define __XIOS_NC4_DATA_INPUT__
3
4/// XIOS headers ///
5#include "xios_spl.hpp"
6#include "data_input.hpp"
7#include "inetcdf4.hpp"
8
9namespace xios
10{
11  class CDomain;
12  class CAxis;
13  class CScalar;
14
15  class CNc4DataInput
16    : protected CINetCDF4
17    , public virtual CDataInput
18  {
19  public:
20    /// Type definitions ///
21    typedef CINetCDF4  SuperClassWriter;
22    typedef CDataInput SuperClass;
23
24    /// Constructors ///
25    CNc4DataInput(const StdString& filename, ep_lib::MPI_Comm comm_file, bool multifile, bool isCollective = true,
26                  bool readMetaDataPar = false, bool ugridConvention = false, const StdString& timeCounterName = "time_counter");
27    CNc4DataInput(const CNc4DataInput& dataInput);       // Not implemented.
28    CNc4DataInput(const CNc4DataInput* const dataInput); // Not implemented.
29
30    /// Destructor ///
31    virtual ~CNc4DataInput(void);
32
33    /// Getters ///
34    const StdString& getFileName(void) const;
35
36  protected:
37    // Read methods
38    virtual StdSize getFieldNbRecords_(CField* field);
39    virtual void readFieldData_(CField* field);
40    virtual void readFieldAttributes_(CField* field, bool readAttributeValues);
41    virtual void closeFile_(void);
42
43  private:
44    void readDomainAttributesFromFile(CDomain* domain, std::list<std::pair<StdString, StdSize> >& dimSizeMap,
45                                      int elementPosition, const StdString& fieldId);
46    void readDomainAttributeValueFromFile(CDomain* domain, std::list<std::pair<StdString, StdSize> >& dimSizeMap,
47                                          int elementPosition, const StdString& fieldId);
48
49    void readAxisAttributesFromFile(CAxis* axis, std::list<std::pair<StdString, StdSize> >& dimSizeMap,
50                                    int elementPosition, const StdString& fieldId);
51    void readAxisAttributeValueFromFile(CAxis* axis, std::list<std::pair<StdString, StdSize> >& dimSizeMap,
52                                        int elementPosition, const StdString& fieldId);
53
54    void readScalarAttributesFromFile(CScalar* scalar, std::list<std::pair<StdString, StdSize> >& dimSizeMap,
55                                      int elementPosition, const StdString& fieldId);
56    void readScalarAttributeValueFromFile(CScalar* scalar, std::list<std::pair<StdString, StdSize> >& dimSizeMap,
57                                          int elementPosition, const StdString& fieldId);
58
59    template <typename T, int Ndim>
60    void readFieldVariableValue(CArray<T, Ndim>& var, const StdString& varId,
61                                const std::vector<StdSize>& nBegin,
62                                const std::vector<StdSize>& nSize,
63                                bool forceIndependent = false);
64
65  private:
66    std::set<StdString> readMetaDataDomains_, readValueDomains_;
67    std::set<StdString> readMetaDataAxis_, readValueAxis_;
68    std::set<StdString> readMetaDataScalar_, readValueScalar_;
69
70  private:
71    /// Private attributes ///
72    ep_lib::MPI_Comm comm_file;
73    const StdString filename;
74    bool isCollective;
75    bool ugridConvention;
76  }; // class CNc4DataInput
77
78template <typename T, int Ndim>
79void CNc4DataInput::readFieldVariableValue(CArray<T,Ndim>& var, const StdString& varId,
80                                           const std::vector<StdSize>& nBegin,
81                                           const std::vector<StdSize>& nSize,
82                                           bool forceIndependent)
83{
84  if (SuperClass::type==MULTI_FILE || !isCollective) return;
85
86  bool openCollective = isCollective;
87  if (forceIndependent) openCollective = !isCollective;
88  switch (SuperClass::type)
89  {
90    case MULTI_FILE:
91      ERROR("CINetCDF4::readFieldVariableValue(...)",
92            << "Only support attributes reading with one_file mode");
93      break;
94    case ONE_FILE:
95    {
96      SuperClassWriter::getData(var, varId, openCollective, 0, &nBegin, &nSize);
97      break;
98    }
99  }
100}
101} // namespace xios
102
103#endif //__XIOS_NC4_DATA_INPUT__
Note: See TracBrowser for help on using the repository browser.