Changeset 999 for XIOS/trunk

Timestamp:

11/21/16 15:31:49 (7 years ago)

Author:

mhnguyen

Message:

Correcting a minor bug on writing a field in file

+) Make sure each element of grid exist before using them

Test
+) On Curie
+) Work

File:

• XIOS/trunk/src/io/nc4_data_output.cpp (modified) (3 diffs)

XIOS/trunk/src/io/nc4_data_output.cpp

@@ -1470 +1470 @@
         std::vector<StdString> domainList = grid->getDomainList();
         std::vector<StdString> axisList   = grid->getAxisList();
-        std::vector<StdString> scalarList = grid->getScalarList();
-
-        CDomain* domain = CDomain::get(domainList[idxDomain]);
+        std::vector<StdString> scalarList = grid->getScalarList();        
 
         StdString timeid  = getTimeCounterName();
@@ -1485 +1483 @@
           if (2 == axisDomainOrder(i))
           {
+            CDomain* domain = CDomain::get(domainList[idxDomain]);
             StdString domId = domain->getDomainOutputName();
             StdString appendDomId  = singleDomain ? "" : "_" + domId ;
@@ -1654 +1653 @@
            if (!SuperClassWriter::useCFConvention)
            {
-             StdString mesh = domain->name;
-             SuperClassWriter::addAttribute("mesh", mesh, &fieldid);
-             StdString location;
-             if (domain->nvertex == 1)
-               location = "node";
-             else if (domain->nvertex == 2)
-               location = "edge";
-             else if (domain->nvertex > 2)
-               location = "face";
-             SuperClassWriter::addAttribute("location", location, &fieldid);
+            if (!domainList.empty())
+            {
+              CDomain* domain = CDomain::get(domainList[0]); // Suppose that we have only domain
+              StdString mesh = domain->name;
+              SuperClassWriter::addAttribute("mesh", mesh, &fieldid);
+              StdString location;
+              if (domain->nvertex == 1)
+                location = "node";
+              else if (domain->nvertex == 2)
+                location = "edge";
+              else if (domain->nvertex > 2)
+                location = "face";
+              SuperClassWriter::addAttribute("location", location, &fieldid);
+            }
+
            }