Changeset 127 for XMLIO_V2/dev/dev_rv/src/XMLIO/NetCDF4_data_output.hpp

Timestamp:

10/07/10 16:17:04 (14 years ago)

Author:

hozdoba

Message:

suite du précédent commit

File:

• XMLIO_V2/dev/dev_rv/src/XMLIO/NetCDF4_data_output.hpp (modified) (11 diffs)

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

@@ -6 +6 @@
 namespace XMLIOSERVER
 {
+
+   static const char* TimeName = "time";
+
    class NetCDF4DataOutput : public AbstractDataOutput
    {
@@ -74 +77 @@
          }
 
-         void writeCoords(const string& id, const Array<float, 1>& cdata)
+         void writeCoords(const string& id, const Array<double, 1>& cdata)
          {
             NcVar *cVar = dataFile->get_var(id.c_str());
@@ -115 +118 @@
             // Ajout de la dimension temporelle non limitée.
             if (withTime)
-               if (NULL == dataFile->add_dim("time"))
+               if (NULL == dataFile->add_dim(TimeName))
                   throw XMLIOUndefinedValueException("Impossible d'ajouter la dimension temporelle !");
          }
@@ -139 +142 @@
                if (!(latVar = dataFile->add_var(latid.c_str(), ncFloat, dataFile->get_dim(latid.c_str()))))
                   throw XMLIOUndefinedValueException("Impossible d'ajouter la variable de latitude !");
+
                if (!(lonVar = dataFile->add_var(lonid.c_str(), ncFloat, dataFile->get_dim(lonid.c_str()))))
                   throw XMLIOUndefinedValueException("Impossible d'ajouter la variable de longitude !");
@@ -147 +151 @@
                hm["units"]         = "degrees_north";
                hm["long_name"]     = "Latitude" ;
-               addStringAttributesToVar(latVar, hm);
+               this->addStringAttributesToVar(latVar, hm);
                hm.clear();
 
@@ -155 +159 @@
                hm["units"]          = "degrees_east";
                hm["long_name"]      = "Longitude" ;
-               addStringAttributesToVar(lonVar, hm);
-               hm.clear();
+               this->addStringAttributesToVar(lonVar, hm);
+               hm.clear();
+
+               if ((*itt)->yvalue.hasValue())
+                  this->writeCoords(lonid, (*itt)->xvalue);
+               else throw XMLIOUndefinedValueException("Les coordonnées de longitude (xvalue) ne sont pas définies pour le domaine \""+domid+"\".");
+
+               if ((*itt)->xvalue.hasValue())
+                  this->writeCoords(latid, (*itt)->yvalue);
+               else throw XMLIOUndefinedValueException("Les coordonnées de latitude (yvalue) ne sont pas définies pour le domaine \""+domid+"\".");
             }
 
@@ -171 +183 @@
                if ((*it)->long_name.hasValue())    hm["long_name"]     = (*it)->long_name ;
                if ((*it)->unit.hasValue())         hm["units"]         = (*it)->unit;
-               addStringAttributesToVar(othvar, hm);
-               hm.clear();
+               this->addStringAttributesToVar(othvar, hm);
+               hm.clear();
+
+               if ((*it)->value.hasValue())
+                  this->writeCoords(axisid, (*it)->value);
+               else throw XMLIOUndefinedValueException("Les coordonnées de l'axe \""+axisid+"\" (value) ne sont pas définies.");
             }
          }
@@ -228 +244 @@
                   {
                      if (!(var = dataFile->add_var(fieldid.c_str(), tvar,
-                        dataFile->get_dim("time"), dataFile->get_dim(latid.c_str()), dataFile->get_dim(lonid.c_str()))))
+                        dataFile->get_dim(TimeName), dataFile->get_dim(latid.c_str()), dataFile->get_dim(lonid.c_str()))))
                         throw XMLIOUndefinedValueException("Impossible d'ajouter le champ "+ field->getId() +" !");
                   }
@@ -234 +250 @@
                   {
                      if (!(var = dataFile->add_var(fieldid.c_str(), tvar,
-                        dataFile->get_dim("time"), dataFile->get_dim(axisid.c_str()),
+                        dataFile->get_dim(TimeName), dataFile->get_dim(axisid.c_str()),
                         dataFile->get_dim(latid.c_str()), dataFile->get_dim(lonid.c_str()))))
                         throw XMLIOUndefinedValueException("Impossible d'ajouter le champ "+ field->getId() +" !");
@@ -258 +274 @@
                if (field->long_name.hasValue())    hm["long_name"]     = field->long_name ;
                if (field->unit.hasValue())         hm["units"]         = field->unit;
-               addStringAttributesToVar(var, hm);
+               this->addStringAttributesToVar(var, hm);
                hm.clear();
             }
@@ -265 +281 @@
       private :
 
-         void addStringAttributesToVar(NcVar * var, const Poco::HashMap<string, string>& attr)
+         void addStringAttributesToVar(NcVar * const var, const Poco::HashMap<string, string>& attr)
          {
             Poco::HashMap<string, string>::ConstIterator it;