source: XIOS/dev/dev_trunk_graph/src/transformation/domain_algorithm/domain_algorithm_interpolate.hpp @ 2019

Last change on this file since 2019 was 2019, checked in by yushan, 9 months ago

Graph intermedia commit to a tmp branch

  • Property svn:executable set to *
File size: 4.1 KB
Line 
1/*!
2   \file domain_algorithm_interpolate_from_file.hpp
3   \author Ha NGUYEN
4   \since 09 July 2015
5   \date 09 Sep 2015
6
7   \brief Algorithm for interpolation on a domain.
8 */
9#ifndef __XIOS_DOMAIN_ALGORITHM_INTERPOLATE_HPP__
10#define __XIOS_DOMAIN_ALGORITHM_INTERPOLATE_HPP__
11
12#include "algorithm_transformation_weight.hpp"
13#include "transformation.hpp"
14#include "nc4_data_output.hpp"
15
16namespace xios {
17
18class CDomain;
19class CInterpolateDomain;
20
21
22/*!
23  \class CDomainAlgorithmInterpolate
24  Reading interpolation from file then apply on a domain
25*/
26class CDomainAlgorithmInterpolate : public CAlgorithmTransformationWeight
27{
28public:
29  CDomainAlgorithmInterpolate(bool isSource, CDomain* domainDestination, CDomain* domainSource, CInterpolateDomain* interpDomain);
30
31  virtual ~CDomainAlgorithmInterpolate() {}
32
33  virtual void apply(const std::vector<std::pair<int,double> >& localIndex,
34                     const double* dataInput,
35                     CArray<double,1>& dataOut,
36                     std::vector<bool>& flagInitial,                     
37                     bool ignoreMissingValue, bool firstPass);
38  virtual void updateData(CArray<double,1>& dataOut);
39
40  static bool registerTrans();
41
42private:
43  void readInterpolationInfo(std::string& filename, std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
44  void writeInterpolationInfo(std::string& filename, std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
45  void processPole(std::map<int,std::vector<std::pair<int,double> > >& interMapValuePole,
46                   int nbGlobalPointOnPole);
47  void computeRemap();
48  void readRemapInfo();
49  void writeRemapInfo(std::map<int,std::vector<std::pair<int,double> > >&);
50  void exchangeRemapInfo(std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
51  void convertRemapInfo(std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
52
53private:
54  CArray<double,1> renormalizationFactor ;
55  CArray<bool,1> allMissing ;
56  bool detectMissingValue ;
57  bool renormalize ;
58  bool quantity ;
59 
60  CInterpolateDomain* interpDomain_;
61  bool writeToFile_;
62  bool readFromFile_;
63  StdString fileToReadWrite_;
64  bool fortranConvention ;
65
66  // class WriteNetCdf;
67  class WriteNetCdf : public CNc4DataOutput
68  {
69  public:
70    WriteNetCdf(const StdString& filename, const MPI_Comm comm);
71    int addDimensionWrite(const StdString& name, const StdSize size = UNLIMITED_DIM);
72    int addVariableWrite(const StdString& name, nc_type type,
73                         const std::vector<StdString>& dim);
74    void endDefinition();
75    void writeDataIndex(const CArray<int,1>& data, const StdString& name,
76                        bool collective, StdSize record,
77                        const std::vector<StdSize>* start = NULL,
78                        const std::vector<StdSize>* count = NULL);
79    void writeDataIndex(const CArray<double,1>& data, const StdString& name,
80                        bool collective, StdSize record,
81                        const std::vector<StdSize>* start = NULL,
82                        const std::vector<StdSize>* count = NULL);
83  };
84
85
86public:
87
88  static CGenericAlgorithmTransformation* create(bool isSource, CGrid* gridDst, CGrid* gridSrc,
89                                                CTransformation<CDomain>* transformation,
90                                                int elementPositionInGrid,
91                                                std::map<int, int>& elementPositionInGridSrc2ScalarPosition,
92                                                std::map<int, int>& elementPositionInGridSrc2AxisPosition,
93                                                std::map<int, int>& elementPositionInGridSrc2DomainPosition,
94                                                std::map<int, int>& elementPositionInGridDst2ScalarPosition,
95                                                std::map<int, int>& elementPositionInGridDst2AxisPosition,
96                                                std::map<int, int>& elementPositionInGridDst2DomainPosition);
97private: 
98  CDomain* domainSrc_ ;
99  CDomain* domainDest_ ;
100
101public:
102  static bool dummyRegistered_;
103
104
105};
106
107}
108#endif // __XIOS_DOMAIN_ALGORITHM_INTERPOLATE_HPP__
Note: See TracBrowser for help on using the repository browser.