source: XIOS/dev/dev_olga/src/node/reduce_domain_to_axis.cpp @ 1158

Last change on this file since 1158 was 980, checked in by mhnguyen, 8 years ago

Minor improvements:

+) Change some type of variable to avoid some error-prone inputs
+) Change some class name to make it meaningful

Test
+) On Curie
+) Okie
File size: 4.0 KB
Line 
1#include "reduce_domain_to_axis.hpp"
2#include "type.hpp"
3#include "axis.hpp"
4#include "domain.hpp"
5
6namespace xios {
7
8  /// ////////////////////// Définitions ////////////////////// ///
9
10  CReduceDomainToAxis::CReduceDomainToAxis(void)
11    : CObjectTemplate<CReduceDomainToAxis>(), CReduceDomainToAxisAttributes(), CTransformation<CAxis>()
12  { /* Ne rien faire de plus */ }
13
14  CReduceDomainToAxis::CReduceDomainToAxis(const StdString & id)
15    : CObjectTemplate<CReduceDomainToAxis>(id), CReduceDomainToAxisAttributes(), CTransformation<CAxis>()
16  { /* Ne rien faire de plus */ }
17
18  CReduceDomainToAxis::~CReduceDomainToAxis(void)
19  {}
20
21  CTransformation<CAxis>* CReduceDomainToAxis::create(const StdString& id, xml::CXMLNode* node)
22  {
23    CReduceDomainToAxis* reduceDomain = CReduceDomainToAxisGroup::get("reduce_domain_to_axis_definition")->createChild(id);
24    if (node) reduceDomain->parse(*node);
25    return static_cast<CTransformation<CAxis>*>(reduceDomain);
26  }
27
28  bool CReduceDomainToAxis::registerTrans()
29  {
30    return registerTransformation(TRANS_REDUCE_DOMAIN_TO_AXIS, CReduceDomainToAxis::create);
31  }
32
33  bool CReduceDomainToAxis::_dummyRegistered = CReduceDomainToAxis::registerTrans();
34
35  //----------------------------------------------------------------
36
37  StdString CReduceDomainToAxis::GetName(void)    { return StdString("reduce_domain_to_axis"); }
38  StdString CReduceDomainToAxis::GetDefName(void) { return StdString("reduce_domain_to_axis"); }
39  ENodeType CReduceDomainToAxis::GetType(void)    { return eReduceDomainToAxis; }
40
41  void CReduceDomainToAxis::checkValid(CAxis* axisDst, CDomain* domainSrc)
42  {
43    if (CDomain::type_attr::unstructured == domainSrc->type)
44      ERROR("CReduceDomainToAxis::checkValid(CAxis* axisDst, CDomain* domainSrc)",
45       << "Domain reduction is only supported for rectilinear or curvillinear grid."
46       << "Domain source " <<domainSrc->getId() << std::endl
47       << "Axis destination " << axisDst->getId());
48
49    int axis_n_glo = axisDst->n_glo;
50    int domain_ni_glo = domainSrc->ni_glo;
51    int domain_nj_glo = domainSrc->nj_glo;
52
53    if (this->operation.isEmpty())
54      ERROR("CReduceDomainToAxis::checkValid(CAxis* axisDst, CDomain* domainSrc)",
55             << "An operation must be defined."
56             << "Domain source " <<domainSrc->getId() << std::endl
57             << "Axis destination " << axisDst->getId());
58
59    if (this->direction.isEmpty())
60      ERROR("CReduceDomainToAxis::checkValid(CAxis* axisDst, CDomain* domainSrc)",
61             << "A direction to apply the operation must be defined. It should be: 'iDir' or 'jDir'"
62             << "Domain source " <<domainSrc->getId() << std::endl
63             << "Axis destination " << axisDst->getId());
64
65    if (this->direction.isEmpty())
66      ERROR("CReduceDomainToAxis::checkValid(CAxis* axisDst, CDomain* domainSrc)",
67             << "A direction to apply the operation must be defined. It should be: 'iDir' or 'jDir'"
68             << "Domain source " <<domainSrc->getId() << std::endl
69             << "Axis destination " << axisDst->getId());
70
71   
72    switch (direction)
73    {
74      case direction_attr::jDir:
75        if (axis_n_glo != domain_ni_glo)
76          ERROR("CReduceDomainToAxis::checkValid(CAxis* axisDst, CDomain* domainSrc)",
77            << "Extract domain along j, axis destination should have n_glo equal to ni_glo of domain source"
78            << "Domain source " <<domainSrc->getId() << " has nj_glo " << domain_ni_glo << std::endl
79            << "Axis destination " << axisDst->getId() << " has n_glo " << axis_n_glo);
80         break;
81
82      case direction_attr::iDir:
83        if (axis_n_glo != domain_nj_glo)
84          ERROR("CReduceDomainToAxis::checkValid(CAxis* axisDst, CDomain* domainSrc)",
85            << "Extract domain along i, axis destination should have n_glo equal to nj_glo of domain source"
86            << "Domain source " <<domainSrc->getId() << " has nj_glo " << domain_nj_glo << std::endl
87            << "Axis destination " << axisDst->getId() << " has n_glo " << axis_n_glo);
88        break;
89
90      default:
91        break;
92    }
93  }
94
95}
Note: See TracBrowser for help on using the repository browser.