XIOS  1.0
Xml I/O Server
 Tout Classes Espaces de nommage Fichiers Fonctions Variables Définitions de type Énumérations Valeurs énumérées Amis Macros
generate_rectilinear_domain.cpp
Aller à la documentation de ce fichier.
1 #include "generate_rectilinear_domain.hpp"
2 #include "type.hpp"
3 
4 namespace xios {
5 
7 
8  CGenerateRectilinearDomain::CGenerateRectilinearDomain(void)
9  : CObjectTemplate<CGenerateRectilinearDomain>(), CGenerateRectilinearDomainAttributes(), CTransformation<CDomain>()
10  { /* Ne rien faire de plus */ }
11 
12  CGenerateRectilinearDomain::CGenerateRectilinearDomain(const StdString & id)
13  : CObjectTemplate<CGenerateRectilinearDomain>(id), CGenerateRectilinearDomainAttributes(), CTransformation<CDomain>()
14  { /* Ne rien faire de plus */ }
15 
16  CGenerateRectilinearDomain::~CGenerateRectilinearDomain(void)
17  {}
18 
19  CTransformation<CDomain>* CGenerateRectilinearDomain::create(const StdString& id, xml::CXMLNode* node)
20  {
21  CGenerateRectilinearDomain* genDomain = CGenerateRectilinearDomainGroup::get("generate_rectilinear_domain_definition")->createChild(id);
22  if (node) genDomain->parse(*node);
23  return static_cast<CTransformation<CDomain>*>(genDomain);
24  }
25 
26  bool CGenerateRectilinearDomain::_dummyRegistered = CGenerateRectilinearDomain::registerTrans();
27  bool CGenerateRectilinearDomain::registerTrans()
28  {
29  registerTransformation(TRANS_GENERATE_RECTILINEAR_DOMAIN, create);
30  }
31 
32  //----------------------------------------------------------------
33 
34  StdString CGenerateRectilinearDomain::GetName(void) { return StdString("generate_rectilinear_domain"); }
35  StdString CGenerateRectilinearDomain::GetDefName(void) { return StdString("generate_rectilinear_domain"); }
36  ENodeType CGenerateRectilinearDomain::GetType(void) { return eGenerateRectilinearDomain; }
37 
38  void CGenerateRectilinearDomain::checkValid(CDomain* domainDst)
39  {
40  const double defaultBndsLonStart = 0;
41  const double defaultBndsLonEnd = 360;
42 
43  const double defaultBndsLatStart = -90;
44  const double defaultBndsLatEnd = 90;
45 
46 
47  int niGlo = domainDst->ni_glo;
48  int njGlo = domainDst->nj_glo;
49 
50 
51  double deltaLon=(defaultBndsLonEnd-defaultBndsLonStart)/niGlo ;
52  if (!lon_start.isEmpty() && !lon_end.isEmpty() ) deltaLon=(lon_end-lon_start)/(niGlo-1) ;
53  if (!bounds_lon_start.isEmpty() && !bounds_lon_end.isEmpty()) deltaLon=(bounds_lon_end-bounds_lon_start)/niGlo ;
54  if (!lon_start.isEmpty() && !bounds_lon_end.isEmpty()) deltaLon=(bounds_lon_end-bounds_lon_start)/(niGlo-0.5) ;
55  if (!bounds_lon_start.isEmpty() && !lon_end.isEmpty()) deltaLon=(bounds_lon_end-bounds_lon_start)/(niGlo-0.5) ;
56 
57  if (lon_start.isEmpty())
58  {
59  if (!bounds_lon_start.isEmpty()) domainDst->lon_start=bounds_lon_start+0.5*deltaLon ;
60  else if (!lon_end.isEmpty()) domainDst->lon_start= lon_end-(niGlo-1.)*deltaLon ;
61  else if (!bounds_lon_end.isEmpty()) domainDst->lon_start=bounds_lon_end-(niGlo-0.5)*deltaLon ;
62  else domainDst->lon_start=defaultBndsLonStart+0.5*deltaLon ;
63  }
64  else domainDst->lon_start=lon_start ;
65 
66  if (bounds_lon_start.isEmpty())
67  {
68  if (!lon_start.isEmpty()) domainDst->bounds_lon_start=lon_start-0.5*deltaLon ;
69  else if (!lon_end.isEmpty()) domainDst->bounds_lon_start= lon_end-(niGlo-0.5)*deltaLon ;
70  else if (!bounds_lon_end.isEmpty()) domainDst->bounds_lon_start=bounds_lon_end-niGlo*deltaLon ;
71  else domainDst->bounds_lon_start=defaultBndsLonStart ;
72  }
73  else domainDst->bounds_lon_start=bounds_lon_start ;
74 
75  if (lon_end.isEmpty())
76  {
77  if (!bounds_lon_end.isEmpty()) domainDst->lon_end=bounds_lon_end-0.5*deltaLon ;
78  else if (!bounds_lon_start.isEmpty()) domainDst->lon_end=bounds_lon_start+(niGlo-0.5)*deltaLon ;
79  else if (!lon_start.isEmpty()) domainDst->lon_end=lon_start+(niGlo-1.)*deltaLon ;
80  else domainDst->lon_end=defaultBndsLonEnd-0.5*deltaLon ;
81  }
82  else domainDst->lon_end=lon_end ;
83 
84  if (bounds_lon_end.isEmpty())
85  {
86  if (!lon_end.isEmpty()) domainDst->bounds_lon_end=lon_end+0.5*deltaLon ;
87  else if (!bounds_lon_start.isEmpty()) domainDst->bounds_lon_end=bounds_lon_start+niGlo*deltaLon ;
88  else if (!lon_start.isEmpty()) domainDst->bounds_lon_end=lon_start+(niGlo-0.5)*deltaLon ;
89  else domainDst->bounds_lon_end=defaultBndsLonEnd ;
90  }
91  else domainDst->bounds_lon_end=bounds_lon_end;
92 
93 
94 
95 
96  double deltaLat=(defaultBndsLatEnd-defaultBndsLatStart)/njGlo ;
97  if (!lat_start.isEmpty() && !lat_end.isEmpty() ) deltaLat=(lat_end-lat_start)/(njGlo-1) ;
98  if (!bounds_lat_start.isEmpty() && !bounds_lat_end.isEmpty()) deltaLat=(bounds_lat_end-bounds_lat_start)/njGlo ;
99  if (!lat_start.isEmpty() && !bounds_lat_end.isEmpty()) deltaLat=(bounds_lat_end-bounds_lat_start)/(njGlo-0.5) ;
100  if (!bounds_lat_start.isEmpty() && !lat_end.isEmpty()) deltaLat=(bounds_lat_end-bounds_lat_start)/(njGlo-0.5) ;
101 
102  if (lat_start.isEmpty())
103  {
104  if (!bounds_lat_start.isEmpty()) domainDst->lat_start=bounds_lat_start+0.5*deltaLat ;
105  else if (!lat_end.isEmpty()) domainDst->lat_start= lat_end-(njGlo-1.)*deltaLat ;
106  else if (!bounds_lat_end.isEmpty()) domainDst->lat_start=bounds_lat_end-(njGlo-0.5)*deltaLat ;
107  else domainDst->lat_start=defaultBndsLatStart+0.5*deltaLat ;
108  }
109  else domainDst->lat_start=lat_start;
110 
111  if (bounds_lat_start.isEmpty())
112  {
113  if (!lat_start.isEmpty()) domainDst->bounds_lat_start=lat_start-0.5*deltaLat ;
114  else if (!lat_end.isEmpty()) domainDst->bounds_lat_start= lat_end-(njGlo-0.5)*deltaLat ;
115  else if (!bounds_lat_end.isEmpty()) domainDst->bounds_lat_start=bounds_lat_end-njGlo*deltaLat ;
116  else domainDst->bounds_lat_start=defaultBndsLatStart ;
117  }
118  else domainDst->bounds_lat_start=bounds_lat_start;
119 
120  if (lat_end.isEmpty())
121  {
122  if (!bounds_lat_end.isEmpty()) domainDst->lat_end=bounds_lat_end-0.5*deltaLat ;
123  else if (!bounds_lat_start.isEmpty()) domainDst->lat_end=bounds_lat_start+(njGlo-0.5)*deltaLat ;
124  else if (!lat_start.isEmpty()) domainDst->lat_end=lat_start+(njGlo-1.)*deltaLat ;
125  else domainDst->lat_end=defaultBndsLatEnd-0.5*deltaLat ;
126  }
127  else domainDst->lat_end=lat_end;
128 
129  if (bounds_lat_end.isEmpty())
130  {
131  if (!lat_end.isEmpty()) domainDst->bounds_lat_end=lat_end+0.5*deltaLat ;
132  else if (!bounds_lat_start.isEmpty()) domainDst->bounds_lat_end=bounds_lat_start+njGlo*deltaLat ;
133  else if (!lat_start.isEmpty()) domainDst->bounds_lat_end=lat_start+(njGlo-0.5)*deltaLat ;
134  else domainDst->bounds_lat_end=defaultBndsLatEnd ;
135  }
136  else domainDst->bounds_lat_end=bounds_lat_end;
137 
138  }
139 
140 }
xios::if
CATCH CDomainAlgorithmReorder::CDomainAlgorithmReorder(CDomain *domainDestination, CDomain *domainSource, CReorderDomain *reorderDomain if)(domainDestination->type!=CDomain::type_attr::rectilinear)
Definition: domain_algorithm_reorder.cpp:46
xios::CGenerateRectilinearDomain::~CGenerateRectilinearDomain
virtual ~CGenerateRectilinearDomain(void)
Destructeur ///.
Definition: generate_rectilinear_domain.cpp:16
StdString
std::string StdString
Definition: xios_spl.hpp:48
xios
#define xios(arg)
Definition: xios_fortran_prefix.hpp:4
xios::CGenerateRectilinearDomain::checkValid
virtual void checkValid(CDomain *axisDest)
Definition: generate_rectilinear_domain.cpp:38
xios::CObjectTemplate
////////////////////// Déclarations ////////////////////// ///
Definition: object_template.hpp:17
xios::CGenerateRectilinearDomain::GetName
static StdString GetName(void)
Accesseurs statiques ///.
Definition: generate_rectilinear_domain.cpp:34
xios::CGenerateRectilinearDomain::CGenerateRectilinearDomain
CGenerateRectilinearDomain(void)
Constructeurs ///.
Definition: generate_rectilinear_domain.cpp:8
xios::CObjectTemplate::parse
virtual void parse(xml::CXMLNode &node)
Definition: object_template_impl.hpp:110
xios::CGenerateRectilinearDomain::create
static CTransformation< CDomain > * create(const StdString &id, xml::CXMLNode *node)
Definition: generate_rectilinear_domain.cpp:19
generate_rectilinear_domain.hpp
xios::ENodeType
enum xios::_node_type ENodeType
////////////////////// Définitions ////////////////////// ///
xios::CTransformation< CDomain >::registerTransformation
static bool registerTransformation(ETranformationType transType, CreateTransformationCallBack createFn)