Changeset 123

Timestamp:

09/27/10 14:01:47 (14 years ago)

Author:

hozdoba

Message:

Début de gestion des opérations sur les champs + quelques modifications

Location:

XMLIO_V2/dev/dev_rv/src/XMLIO

Files:

• NetCDF4_data_output.hpp (modified) (6 diffs)
• abstract_calendar.hpp (modified) (1 diff)
• abstract_data_treatment.hpp (modified) (1 diff)
• c_interface.cpp (modified) (2 diffs)
• context.hpp (modified) (2 diffs)
• date.hpp (modified) (1 diff)
• duration.hpp (modified) (2 diffs)
• field.hpp (modified) (4 diffs)
• field_operation.hpp (added)
• file_attribut.hpp (modified) (1 diff)
• grid.hpp (modified) (3 diffs)
• xml_parser.hpp (modified) (1 diff)
• xmlio.hpp (modified) (5 diffs)

XMLIO_V2/dev/dev_rv/src/XMLIO/NetCDF4_data_output.hpp

@@ -51 +51 @@
               timeinfo = localtime ( &rawtime );
 
-              strftime (buffer, 100,"%Y-%b-%d %H:%M:%S %Z",timeinfo);
+              strftime (buffer, 100, "%Y-%b-%d %H:%M:%S %Z", timeinfo);
 
               return (string(buffer));
@@ -141 +141 @@
          virtual void setVars(void)
          {
+            NcType tvar;
             NcVar *var  = NULL;
             Poco::HashMap<string, string> hm;
@@ -156 +157 @@
                   if (string(field->operation).compare("once") == 0) wtime = false;
 
+               tvar = ncFloat;
+               if (field->prec.hasValue())
+                  if (field->prec == 8) tvar = ncDouble;
+
                if (wtime)
                {
                   if (lonlat) // 2D spatio + temps
                   {
-                     if (!(var = dataFile->add_var(bfield->getId().c_str(), ncFloat,
+                     if (!(var = dataFile->add_var(bfield->getId().c_str(), tvar,
                         dataFile->get_dim("time"), dataFile->get_dim("lon"), dataFile->get_dim("lat"))))
                         throw XMLIOUndefinedValueException("Impossible d'ajouter le champ "+ field->getId() +" !");
@@ -166 +171 @@
                   else // 3D spatio + temps
                   {
-                     if (!(var = dataFile->add_var(bfield->getId().c_str(), ncFloat,
+                     if (!(var = dataFile->add_var(bfield->getId().c_str(), tvar,
                         dataFile->get_dim("time"), dataFile->get_dim("lon"), dataFile->get_dim("lat"),
                         dataFile->get_dim(field->getGrid()->getRelAxis()->getId().c_str()))))
@@ -176 +181 @@
                   if (lonlat) // 2D spatio sans temps
                   {
-                     if (!(var = dataFile->add_var(bfield->getId().c_str(), ncFloat,
+                     if (!(var = dataFile->add_var(bfield->getId().c_str(), tvar,
                         dataFile->get_dim("lon"), dataFile->get_dim("lat"))))
                         throw XMLIOUndefinedValueException("Impossible d'ajouter le champ "+ field->getId() +" !");
@@ -182 +187 @@
                   else // 3D spatio sans temps
                   {
-                     if (!(var = dataFile->add_var(bfield->getId().c_str(), ncFloat,
+                     if (!(var = dataFile->add_var(bfield->getId().c_str(), tvar,
                         dataFile->get_dim("lon"), dataFile->get_dim("lat"),
                         dataFile->get_dim(field->getGrid()->getRelAxis()->getId().c_str()))))

XMLIO_V2/dev/dev_rv/src/XMLIO/abstract_calendar.hpp

@@ -17 +17 @@
             return (NoLeapMonthLength[d.getMonth()-1]);
          }
+
+         friend std::ostream& operator<<(std::ostream& out, AbstractCalendar& c)
+         { out << "[type: " << c.getId() << ", start: " << c.getInitDate() << ", current: " << c.getCurrentDate() << "]"; return (out); }
 
          virtual int getYearTotalLength(const Date& d) const { return (365 * 86400); } // Retourne la durée d'une année en seconde.

XMLIO_V2/dev/dev_rv/src/XMLIO/abstract_data_treatment.hpp

@@ -36 +36 @@
          AbstractDataTreatment(Context* const _ctxt) : currentContext(_ctxt), enabledFiles()
          {
-            std::cout << "(Message temporaire) Résolution des héritages ..." << std::endl;
+            //std::cout << "(Message temporaire) Résolution des héritages ..." << std::endl;
             // Résolution des héritages pour le context actuel (Context::ResolveInheritance() résoud les héritages pour tous les contextes).
             solveAllInheritance();
 
             // Initialisation du vecteur 'enabledFiles' contenant la liste des fichiers à sortir.
-            std::cout << "(Message temporaire) Initialisation du vecteur enabledFiles ..." << std::endl;
+            //std::cout << "(Message temporaire) Initialisation du vecteur enabledFiles ..." << std::endl;
             findEnabledFiles();
 
             // Recherche des champs à sortir (enable à true + niveau de sortie correct) pour chaque fichier précédemment listé.
-            std::cout << "(Message temporaire) Recherche des champs à sortir ..." << std::endl;
+            //std::cout << "(Message temporaire) Recherche des champs à sortir ..." << std::endl;
             findAllEnabledFields();
 
             // Résolution des références de grilles pour chacun des champs.
-            std::cout << "(Message temporaire) Résolution des références de grilles ..." << std::endl;
+            //std::cout << "(Message temporaire) Résolution des références de grilles ..." << std::endl;
             solveAllGridRef();
 
-            std::cout << "(Message temporaire) fin traitement sorties ..." << std::endl;
+            //std::cout << "(Message temporaire) fin traitement sorties ..." << std::endl;
          }
 

XMLIO_V2/dev/dev_rv/src/XMLIO/c_interface.cpp

@@ -8 +8 @@
 extern "C" void main_c_(void) ;
 
+// Cette fonction remplit un tableau à N dimensions de valeurs de type T (T = double ou float) de la maniÚre suivante :
+// 0.0 1.0 2.0 3.0, puis 0.1 1.1 2.1 3.1 ... au second appel de la fonction, etc.
+template <class T>
+   static void updateDataTest(T begin, T end)
+{
+   static int called = 0; int i = 0;
+   for (T it = begin; it != end; it++, i++)
+      (*it) = i + called * 0.1;
+   called++;
+}
+
 void main_c_ (void)
 {
    try
    {
-      //string file("/local/XMLIOSERVER_DEV/dev_rv/test/iodef_test.xml");
-      //string file("/local/XMLIOSERVER_DEV/dev_rv/iodef_test.xml");
-      string file("/local/XMLIOSERVER_DEV/dev_rv/test/iodef_simple_test.xml");
-
+      Array<double, 3> arr(3, 3, 3) ; // Les données.
+      string file("/local/XMLIOSERVER_DEV/dev_rv/test/iodef_simple_test.xml"); // Le fichier de définition.
       ifstream istr( file.c_str() , ifstream::in );
 
       // On commence la lecture du flux de donnée xml qui doit posséder pour racine un unique noeud nommé "simulation".
       XMLNode node = XMLNode::CreateNode(istr, Context::GetRootName());
+
       // On parse le fichier xml noeud par noeud (ie on construit dynamiquement notre arbre d'objets).
       XMLParser::Parse(node);
-
-      // On résoud les héritages descendants ainsi que les héritages par référence.
-      //Context::ResolveInheritance(); // Plus nécessaire!!
 
       // On se place dans le contexte 'context1'.
       Context::SetCurrentContext("context1") ;
 
-      std::cout << " * ----------- <<<TOP>>> Début du traitement ----------- * " << std::endl;
+      std::cout << " * ----------- Début du traitement ----------- * " << std::endl;
+
+      //Affichage de la date initiale (pour vérification).
+      std::cout << *Context::GetObject("context1")->getCalendar() << std::endl;
+
       // On crée une instance de traitement de données dans le format NetCDF4.
       AbstractDataTreatment *dtreat = new NetCDF4DataTreatment(Context::GetObject("context1"));
-
-      // Affichage de la liste des fichiers à sortir (par id).
-      //dtreat->printEnabledFilesId();
 
       NetCDF4DataOutput dop(CFile::GetObject("fichier1"));
       dop.createFileAndHeader();
 
-      // On écrit l'arborescence résultante du traitement sur la sortie.
+      for(int tcourant = 0, tfinal = 9; tcourant <= tfinal; tcourant++)
+      {
+         // Mise à jour de la date.
+
+
+         // Mise à jour des données.
+         updateDataTest<Array<double, 3>::iterator>(arr.begin(), arr.end());
+
+         // Ecriture des données.
+         // dtreat->writeData("Premier_champ", arr);
+
+      }
+
+      // On écrit l'arborescence résultante du traitement sur la sortie de log.
       //Context::ShowTree(std::clog);
 
@@ -52 +73 @@
    }
 
-// return (0);
+   // return (0);
 }
 
 
-      /* CODE TEST POUR CALENDRIER
-      AbstractCalendar* calendar = new GregorianCalendar();
-
-      Date dd (*calendar, 2000, 05, 03, 00, 00, 00) ;
-      Date ddd(*calendar, 2030, 11, 02, 10, 33, 40) ;
-
-      Date d1 (*calendar, 2007, 11, 12, 00, 00, 00) ;
-      Date d2 (*calendar, 2008, 11, 12, 00, 00, 00) ;
-
-      string sttr("12/07/2005 - 10:11:30");
-      string stttr("111s 14.36 mi 7y 3.9y 3.4 mg");
-      Date d3 = Date::FromString(sttr, *calendar);
-
-      // Les mois sont considérés comme des entiers mais stockés comme des flottants.
-      Duration sc0 = Second * -65.3 + Minute * 117.5   + Hour * 65.9  + Day * 417.8 + Month * 16.0  - Year * 17.5;
-      Duration sc1 = Second * 175.3 - Minute * 13109.6 + Hour * 135.1 - Day * 616.4 + Month * 236.0 + Year * 33.9;
-      Duration sc2 = 60 * Day;
-
-      Duration sc3 = Duration::FromString(stttr);
-
-      std::cout << calendar->getInitDate()    << std::endl;
-      std::cout << calendar->getCurrentDate() << std::endl;
-
-      std::cout << boolalpha << (calendar->getInitDate() != calendar->getCurrentDate()) << std::endl;
-
-      std::cout << (Time)d2 - (Time)(d1 + 365 * Day) << std::endl;
-      std::cout << sc0 << std::endl << sc1 << std::endl << sc1 - sc0 << std::endl;
-      std::cout << ddd - dd << std::endl;
-      std::cout << dd << std::endl;
-      std::cout << (Time)dd << std::endl;
-      std::cout << calendar->getNbSecond(dd) << std::endl;
-
-      std::cout << (dd+sc2) << std::endl;
-      std::cout << (d3)  << std::endl;
-
-      std::cout << (sc3) << std::endl;
-
-      delete calendar;
-      */
-

XMLIO_V2/dev/dev_rv/src/XMLIO/context.hpp

@@ -94 +94 @@
                   if (attributes.end() != attributes.find("id"))
                   { WARNING("Le noeud de définition possÚde un identifiant, ce dernier ne sera pas pris en compte lors du traitement !"); }
-
-                  if (attributes.end() != attributes.find("calendar_type") && attributes.end() != attributes.find("start_date"))
-                     setCalendar(attributes["calendar_type"], attributes["start_date"]);
-
 
                   if (name.compare(FieldDefinition::GetDefName())  == 0) // Parsing pour la définition des champs.
@@ -153 +149 @@
          GridDefinition   * getGridDefinition(void)   const { return (this->gridDef  ); }
          DomainDefinition * getDomainDefinition(void) const { return (this->domainDef); }
+
+         void updateCurrentDate(const Duration* const dr = NULL)
+         {
+            static Duration dur = Hour;
+            if (dr != NULL) dur = *dr;
+            Date newCurrentDate = ccalendar->getCurrentDate() + dur;
+            ccalendar->getCurrentDate() = newCurrentDate;
+         }
 
          AbstractCalendar * getCalendar(void) const { return (this->ccalendar ); }

XMLIO_V2/dev/dev_rv/src/XMLIO/date.hpp

@@ -15 +15 @@
          }
 
-         Date(const Date& d) : relCalendar(d.relCalendar)
+         Date(const Date& d) : relCalendar(d.getRelCalendar())
             , year(d.year), month(d.month), day(d.day), hour(d.hour), minute(d.minute), second(d.second)
          {
             if(!checkDate())
                WARNING("La date initialisée a été modifiée car elle était incorrecte par rapport au calendrier souhaité.");
+         }
+
+         Date& operator=(const Date& d)
+         {
+            // relCalendar = d.getRelCalendar(); << inutile si fonction bien utilisée
+            year = d.year; month = d.month; day = d.day;
+            hour = d.hour; minute = d.minute; second = d.second;
+            return (*this);
          }
 

XMLIO_V2/dev/dev_rv/src/XMLIO/duration.hpp

@@ -26 +26 @@
          friend std::istream& operator>>(std::istream& in, struct _duration& d)
          {
+            d.year = d.month = d.day = d.hour = d.minute = d.second = 0.0;
             double v = 1.0;
             char   c = '/';
@@ -77 +78 @@
                          Hour   = {0.0, 0.0, 0.0, 1.0, 0.0, 0.0},
                          Minute = {0.0, 0.0, 0.0, 0.0, 1.0, 0.0},
-                         Second = {0.0, 0.0, 0.0, 0.0, 0.0, 1.0};
+                         Second = {0.0, 0.0, 0.0, 0.0, 0.0, 1.0},
+                         NoneDu = {0.0, 0.0, 0.0, 0.0, 0.0, 0.0};
 
 } // namespace XMLIOSERVER

XMLIO_V2/dev/dev_rv/src/XMLIO/field.hpp

@@ -15 +15 @@
          CField(void) : ObjectTemplate<CField>(), FieldAttribut(), grid(NULL)
          {/* Ne rien faire de plus */}
-         CField(const string& _id) : ObjectTemplate<CField>(_id), FieldAttribut(), grid(NULL)
+         CField(const string& _id) : ObjectTemplate<CField>(_id), FieldAttribut(), grid(NULL), foper(NULL)
          {/* Ne rien faire de plus */}
 
@@ -23 +23 @@
 
          const CGrid* getGrid(void) const { return (grid); }
+         const Array<double, 1>& getData(void) const { return (data); }
+         const FieldOperation* getFieldOperation(void) const { return (foper); }
 
          virtual CField* getReference(void) const
@@ -28 +30 @@
             if(!field_ref.hasValue()) return (NULL);
             if (!CField::HasObject(field_ref))
-            { WARNING("Référence invalide sur l'objet "+GetName()+" nommé \""+((string)field_ref)+"\""); return (NULL);}
+            { WARNING("Référence invalide sur l'objet "+GetName()+" nommé \""+((string)field_ref)+"\""); return (NULL); }
 
             return (CField::GetObject(field_ref));
@@ -39 +41 @@
 
          CGrid* grid ;
+         FieldOperation* foper;
+         Array<double, 1> data;
 
    }; // class CField

XMLIO_V2/dev/dev_rv/src/XMLIO/file_attribut.hpp

@@ -10 +10 @@
          DECLARE_ATTR(name, string) ;
          DECLARE_ATTR(description, string) ;
-         DECLARE_ATTR(output_freq, int) ;
+         DECLARE_ATTR(output_freq, Duration) ;
          DECLARE_ATTR(output_level, int) ;
          DECLARE_ATTR(enabled, bool);

XMLIO_V2/dev/dev_rv/src/XMLIO/grid.hpp

@@ -44 +44 @@
          { /* Ne rien faire de plus */ }
 
-
       private:
 
@@ -52 +51 @@
          CDomain* domain ;
 
-         Array<int,1> storeIndex ;
-         Array<int,1> out_i_index ;
-         Array<int,1> out_j_index ;
-         Array<int,1> out_l_index ;
+         Array<int, 1> storeIndex ;
+         Array<int, 1> out_i_index ;
+         Array<int, 1> out_j_index ;
+         Array<int, 1> out_l_index ;
 
    }; // class CGrid
-
 
    CGrid* CGrid::CreateObject(const CDomain* const a_domain, const CAxis* const a_axis)
@@ -66 +64 @@
 
       grid->domain_ref = a_domain->getId() ;
-      grid->axis_ref = a_axis->getId() ;
+      grid->axis_ref   = a_axis  ->getId() ;
 
       return (grid);

XMLIO_V2/dev/dev_rv/src/XMLIO/xml_parser.hpp

@@ -37 +37 @@
                   Context& context = *Context::CreateObject(attributes["id"]);
 
+                  if (attributes.end() != attributes.find("calendar_type") && attributes.end() != attributes.find("start_date"))
+                     context.setCalendar(attributes["calendar_type"], attributes["start_date"]);
+
                   context.parse(_node);
                   attributes.clear();

XMLIO_V2/dev/dev_rv/src/XMLIO/xmlio.hpp

@@ -40 +40 @@
 #include "abstract_object.hpp"
 
-// Gestion du calendrier
+// Gestion du calendrier.
 #include "duration.hpp"
 #include "date.hpp"
@@ -47 +47 @@
 #include "calendar_utils.hpp"
 
-// Gestion de l'écriture des données
+// Gestion de l'écriture des données.
 #include "abstract_data_output.hpp"
 
+//....
 #include "xml_node.hpp"
 #include "object_template.hpp"
@@ -59 +60 @@
 #include "field_attribut.hpp"
 #include "file_attribut.hpp"
+
+// Traitement des opérations sur les champs
+#include "field_operation.hpp"
 
 #include "group_template.hpp"
@@ -70 +74 @@
 #include "file.hpp"
 
+// Résolution des références aux grilles pour les champs.
 #include "field_gridrefsolver.hpp"
 
+//....
 #include "context.hpp"
 #include "xml_parser.hpp"
@@ -79 +85 @@
 #include "NetCDF4_data_output.hpp"
 
-// Traitement des données
+// Traitement des données (abstract puis selon types de sortie)
 #include "abstract_data_treatment.hpp"
 #include "NetCDF4_data_treatment.hpp"