source: XIOS/trunk/src/node/domain.hpp @ 676

Last change on this file since 676 was 676, checked in by rlacroix, 6 years ago

Add support for indexed output.

If the new field attribute "indexed_output" is set to true and a mask is defined (either at grid, domain or axis level), the indexed data will be outputed instead of the full data with missing values.

See http://cfconventions.org/Data/cf-conventions/cf-conventions-1.5/build/cf-conventions.html#compression-by-gathering for more information.

  • 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: 6.9 KB
Line 
1#ifndef __XIOS_CDomain__
2#define __XIOS_CDomain__
3
4/// XIOS headers ///
5#include "xios_spl.hpp"
6#include "group_factory.hpp"
7
8#include "declare_group.hpp"
9#include "declare_ref_func.hpp"
10#include "event_client.hpp"
11#include "event_server.hpp"
12#include "buffer_in.hpp"
13#include "array_new.hpp"
14#include "attribute_array.hpp"
15#include "attribute_enum.hpp"
16#include "transformation.hpp"
17#include "transformation_enum.hpp"
18
19namespace xios {
20
21   /// ////////////////////// Déclarations ////////////////////// ///
22
23   class CDomainGroup;
24   class CDomainAttributes;
25   class CDomain;
26
27   ///--------------------------------------------------------------
28
29   // Declare/Define CDomainAttribute
30   BEGIN_DECLARE_ATTRIBUTE_MAP(CDomain)
31#  include "domain_attribute.conf"
32   END_DECLARE_ATTRIBUTE_MAP(CDomain)
33
34   ///--------------------------------------------------------------
35
36   class CDomain
37      : public CObjectTemplate<CDomain>
38      , public CDomainAttributes
39   {
40         enum EEventId
41         {
42           EVENT_ID_SERVER_ATTRIBUT, EVENT_ID_INDEX, EVENT_ID_LON, EVENT_ID_LAT, EVENT_ID_AREA
43         } ;
44
45         /// typedef ///
46         typedef CObjectTemplate<CDomain>   SuperClass;
47         typedef CDomainAttributes SuperClassAttribute;
48
49      public:
50
51         typedef CDomainAttributes RelAttributes;
52         typedef CDomainGroup      RelGroup;
53         typedef CTransformation<CDomain>::TransformationMapTypes TransMapTypes;
54
55         /// Constructeurs ///
56         CDomain(void);
57         explicit CDomain(const StdString & id);
58         CDomain(const CDomain & domain);       // Not implemented yet.
59         CDomain(const CDomain * const domain); // Not implemented yet.
60
61         static CDomain* createDomain();
62         void duplicateAttributes(CDomain* domain);
63
64         virtual void parse(xml::CXMLNode & node);
65
66         /// Vérifications ///
67         void checkAttributes(void);
68
69         void checkAttributesOnClient();
70         void checkAttributesOnClientAfterTransformation();
71
72         void checkEligibilityForCompressedOutput(void);
73
74         void sendCheckedAttributes();
75
76         bool hasTransformation();
77         void solveInheritanceTransformation();
78         TransMapTypes getAllTransformations();
79
80      public:
81         const std::set<StdString> & getRelFiles(void) const;
82         bool IsWritten(const StdString & filename) const;
83         bool isWrittenCompressed(const StdString& filename) const;
84
85         const std::vector<int>& getIndexesToWrite(void) const;
86         int getNumberWrittenIndexes() const;
87         int getTotalNumberWrittenIndexes() const;
88         int getOffsetWrittenIndexes() const;
89
90         bool isEmpty(void) const;
91         bool isDistributed(void) const;
92         bool isCompressible(void) const;
93
94         int global_zoom_ibegin, global_zoom_ni;
95         int global_zoom_jbegin, global_zoom_nj;
96
97         int ni_srv,ibegin_srv,iend_srv ;
98         int zoom_ni_srv,zoom_ibegin_srv,zoom_iend_srv ;
99
100         int nj_srv,jbegin_srv,jend_srv ;
101         int zoom_nj_srv,zoom_jbegin_srv,zoom_jend_srv ;
102
103         CArray<double, 1> lonvalue_srv, latvalue_srv ;
104         CArray<double, 2> bounds_lon_srv, bounds_lat_srv ;
105         CArray<double, 1> lonvalue_client, latvalue_client;
106         CArray<double, 2> bounds_lon_client, bounds_lat_client;
107         CArray<double, 1> area_srv;
108
109        vector<int> connectedServer ; // list of connected server
110        vector<int> nbSenders ; // for each communication with a server, number of communicating client
111        vector<int> nbDataSrv ; // size of data to send to each server
112        vector< vector<int> > i_indSrv ; // for each server, i global index to send
113        vector< vector<int> > j_indSrv ; // for each server, j global index to send
114
115      public:
116         /// Mutateur ///
117         void addRelFile(const StdString & filename);
118         void addRelFileCompressed(const StdString& filename);
119         void completeLonLatClient(void);
120         void sendServerAttribut(void) ;
121         void sendLonLatArea(void);
122         void computeConnectedServer(void) ;
123
124         static bool dispatchEvent(CEventServer& event);
125         static void recvServerAttribut(CEventServer& event);
126         static void recvIndex(CEventServer& event);
127         static void recvLon(CEventServer& event);
128         static void recvLat(CEventServer& event);
129         static void recvArea(CEventServer& event);
130         void recvServerAttribut(CBufferIn& buffer);
131         void recvIndex(int rank, CBufferIn& buffer);
132         void recvLon(int rank, CBufferIn& buffer);
133         void recvLat(int rank, CBufferIn& buffer);
134         void recvArea(int rank, CBufferIn& buffer);
135
136         /// Destructeur ///
137         virtual ~CDomain(void);
138
139         /// Accesseurs statiques ///
140         static StdString GetName(void);
141         static StdString GetDefName(void);
142
143         static ENodeType GetType(void);
144         const std::map<int, vector<size_t> >& getIndexServer() const;
145         CArray<int, 2> local_mask;
146         bool isCurvilinear ;
147         bool hasBounds ;
148         bool hasArea;
149         bool hasLonLat;
150      private:
151         void checkDomain(void);
152         void checkLocalIDomain(void);
153         void checkLocalJDomain(void);
154
155         void checkMask(void);
156         void checkDomainData(void);
157         void checkCompression(void);
158
159         void checkBounds(void);
160         void checkArea(void);
161         void checkLonLat();
162
163         void checkTransformations();
164         void setTransformations(const TransMapTypes&);
165         void computeNGlobDomain();
166
167         void sendIndex();
168         void sendArea();
169         void sendLonLat();
170
171       private:
172         bool isChecked;
173         std::set<StdString> relFiles, relFilesCompressed;
174         bool isClientChecked; // Verify whether all attributes of domain on the client side are good
175         bool isClientAfterTransformationChecked;
176         std::map<int, CArray<int,1> > indiSrv, indjSrv;
177         std::map<int,int> nbConnectedClients_; // Mapping of number of communicating client to a server
178         std::map<int, vector<size_t> > indSrv_; // Global index of each client sent to server
179         std::map<int, vector<int> > indWrittenSrv_; // Global written index of each client sent to server
180         std::vector<int> indexesToWrite;
181         int numberWrittenIndexes_, totalNumberWrittenIndexes_, offsetWrittenIndexes_;
182         std::vector<int> connectedServerRank_;
183         bool isDistributed_;
184         //! True if and only if the data defined on the domain can be outputted in a compressed way
185         bool isCompressible_;
186         TransMapTypes transformationMap_;
187         std::vector<int> nGlobDomain_;
188         bool isUnstructed_;
189
190         DECLARE_REF_FUNC(Domain,domain)
191
192   }; // class CDomain
193
194   ///--------------------------------------------------------------
195
196   // Declare/Define CDomainGroup and CDomainDefinition
197   DECLARE_GROUP(CDomain);
198
199   ///--------------------------------------------------------------
200
201} // namespace xios
202
203#endif // __XIOS_CDomain__
Note: See TracBrowser for help on using the repository browser.