source: XIOS/trunk/src/transformation/domain_algorithm_interpolate.hpp @ 773

Last change on this file since 773 was 743, checked in by mhnguyen, 9 years ago

Processing pole of rectangular grid with latitude = 90

+) Use average value of all points with latitude 90 (-90)

Test
+) On Curie
+) test_remap pass

File size: 1.3 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 "domain_algorithm_transformation.hpp"
13#include "domain.hpp"
14#include "interpolate_domain.hpp"
15
16namespace xios {
17/*!
18  \class CDomainAlgorithmInterpolate
19  Reading interpolation from file then apply on a domain
20*/
21class CDomainAlgorithmInterpolate : public CDomainAlgorithmTransformation
22{
23public:
24  CDomainAlgorithmInterpolate(CDomain* domainDestination, CDomain* domainSource, CInterpolateDomain* interpDomain);
25
26  virtual ~CDomainAlgorithmInterpolate() {}
27
28  virtual void computeIndexSourceMapping();
29
30private:
31  void readInterpolationInfo(std::string& filename, std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
32  void processPole(std::map<int,std::vector<std::pair<int,double> > >& interMapValuePole,
33                   int nbGlobalPointOnPole);
34  void computeRemap();
35  void readRemapInfo();
36  void exchangeRemapInfo(const std::map<int,std::vector<std::pair<int,double> > >& interpMapValue);
37private:
38  CInterpolateDomain* interpDomain_;
39
40};
41
42}
43#endif // __XIOS_DOMAIN_ALGORITHM_INTERPOLATE_HPP__
Note: See TracBrowser for help on using the repository browser.