source: XIOS/dev/XIOS_DEV_CMIP6/src/transformation/domain_algorithm_interpolate.hpp @ 1269

Last change on this file since 1269 was 1269, checked in by ymipsl, 3 years ago

interpolate_domain : add new attribute "read_write_convention" which msut be "c" or fortran. The default remain fortran (global indices begin at 1)

YM

File size: 4.0 KB
RevLine 
[689]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 "domain_algorithm_transformation.hpp"
[933]13#include "transformation.hpp"
[1021]14#include "nc4_data_output.hpp"
[689]15
16namespace xios {
[933]17
18class CDomain;
19class CInterpolateDomain;
20
[1021]21
[689]22/*!
23  \class CDomainAlgorithmInterpolate
24  Reading interpolation from file then apply on a domain
25*/
26class CDomainAlgorithmInterpolate : public CDomainAlgorithmTransformation
27{
28public:
29  CDomainAlgorithmInterpolate(CDomain* domainDestination, CDomain* domainSource, CInterpolateDomain* interpDomain);
30
31  virtual ~CDomainAlgorithmInterpolate() {}
32
[1264]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
[933]40  static bool registerTrans();
[827]41protected:
42  void computeIndexSourceMapping_(const std::vector<CArray<double,1>* >& dataAuxInputs);
[689]43
44private:
45  void readInterpolationInfo(std::string& filename, std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
[1021]46  void writeInterpolationInfo(std::string& filename, std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
[743]47  void processPole(std::map<int,std::vector<std::pair<int,double> > >& interMapValuePole,
48                   int nbGlobalPointOnPole);
[689]49  void computeRemap();
50  void readRemapInfo();
[1021]51  void writeRemapInfo(std::map<int,std::vector<std::pair<int,double> > >&);
[856]52  void exchangeRemapInfo(std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
[1173]53  void convertRemapInfo(std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
[1021]54
[689]55private:
[1264]56  CArray<double,1> renormalizationFactor ;
57  bool detectMissingValue ;
58  bool renormalize ;
59  bool quantity ;
60 
[689]61  CInterpolateDomain* interpDomain_;
[1021]62  bool writeToFile_;
63  bool readFromFile_;
64  StdString fileToReadWrite_;
[1269]65  bool fortranConvention ;
[689]66
[1021]67  // class WriteNetCdf;
68  class WriteNetCdf : public CNc4DataOutput
69  {
70  public:
71    WriteNetCdf(const StdString& filename, const MPI_Comm comm);
72    int addDimensionWrite(const StdString& name, const StdSize size = UNLIMITED_DIM);
73    int addVariableWrite(const StdString& name, nc_type type,
74                         const std::vector<StdString>& dim);
[1158]75    void endDefinition();
[1021]76    void writeDataIndex(const CArray<int,1>& data, const StdString& name,
77                        bool collective, StdSize record,
78                        const std::vector<StdSize>* start = NULL,
79                        const std::vector<StdSize>* count = NULL);
80    void writeDataIndex(const CArray<double,1>& data, const StdString& name,
81                        bool collective, StdSize record,
82                        const std::vector<StdSize>* start = NULL,
83                        const std::vector<StdSize>* count = NULL);
84  };
85
86
[933]87private:
88
89  static CGenericAlgorithmTransformation* create(CGrid* gridDst, CGrid* gridSrc,
90                                                CTransformation<CDomain>* transformation,
91                                                int elementPositionInGrid,
92                                                std::map<int, int>& elementPositionInGridSrc2ScalarPosition,
93                                                std::map<int, int>& elementPositionInGridSrc2AxisPosition,
94                                                std::map<int, int>& elementPositionInGridSrc2DomainPosition,
95                                                std::map<int, int>& elementPositionInGridDst2ScalarPosition,
96                                                std::map<int, int>& elementPositionInGridDst2AxisPosition,
97                                                std::map<int, int>& elementPositionInGridDst2DomainPosition);
[689]98};
99
100}
101#endif // __XIOS_DOMAIN_ALGORITHM_INTERPOLATE_HPP__
Note: See TracBrowser for help on using the repository browser.