Changeset 2200 for XIOS/trunk/src/node/file.cpp

Timestamp:

08/19/21 13:28:53 (3 years ago)

Author:

ymipsl

Message:

Bugs fix for UGRID convention output

• global indexation was not taking into account
• coherence problem in connectivity for node, edge and face mesh
• add UGRID testcase based on sphere cube

YM

File:

• XIOS/trunk/src/node/file.cpp (modified) (3 diffs)

XIOS/trunk/src/node/file.cpp

@@ -867 +867 @@
    \brief Sorting domains with the same name (= describing the same mesh) in the decreasing order of nvertex for UGRID files.
    This insures that the domain with the highest nvertex is written first and thus all known mesh connectivity is generated at once by this domain.
+   Secondly, we associate the UGRID mesh to the domain in incresing order, i.e. nvertex=1 first, nvertex=2 and then nvertex>2.
+   In this case the connectivity of each component of the mesh (node, edge and face) are fully coherent.
    */
    void CFile::sortEnabledFieldsForUgrid()
@@ -874 +876 @@
      std::vector<int> domainNvertices;
      std::vector<StdString> domainNames;
+     std::map<string, tuple<CDomain*,CDomain*,CDomain*>> registeredDomains ;
 
      for (int i = 0; i < size; ++i)
@@ -916 +919 @@
          domainNvertices.push_back(nvertex);
        }
+
+       if (nvertex==1)  std::get<0>(registeredDomains[domainName])=domain[0] ;
+       else if (nvertex==2) std::get<1>(registeredDomains[domainName])=domain[0] ;
+       else  std::get<2>(registeredDomains[domainName])=domain[0] ;
+     }
+
+     for(auto& it:registeredDomains)
+     {
+       list<CDomain*> domains ;
+       string domainName=it.first ;
+
+       if (std::get<0>(it.second)!=nullptr) domains.push_back(std::get<0>(it.second)) ;
+       if (std::get<1>(it.second)!=nullptr) domains.push_back(std::get<1>(it.second)) ; 
+       if (std::get<2>(it.second)!=nullptr) domains.push_back(std::get<2>(it.second)) ; 
+       
+       // for each component of a given mesh (i.e. domains with same name but different number of vertices)
+       // associate the UGRID mesh in increasing order
+       for(auto& domain : domains )
+       {
+         domain-> computeWrittenIndex();
+         CArray<int, 1>& indexToWrite = domain->localIndexToWriteOnServer;
+         int nbWritten = indexToWrite.numElements();
+         CArray<double,1> writtenLat, writtenLon;
+         CArray<double,2> writtenBndsLat, writtenBndsLon;
+
+         writtenLat.resize(nbWritten);
+         writtenLon.resize(nbWritten);
+         for (int idx = 0; idx < nbWritten; ++idx)
+         {
+           writtenLat(idx) = domain->latvalue(indexToWrite(idx));
+           writtenLon(idx) = domain->lonvalue(indexToWrite(idx));
+         }
+    
+         int nvertex = domain->nvertex, idx;
+        if (nvertex>1)
+         {
+           writtenBndsLat.resize(nvertex, nbWritten);
+           writtenBndsLon.resize(nvertex, nbWritten);
+           CArray<double,2>& boundslat = domain->bounds_latvalue;
+           CArray<double,2>& boundslon = domain->bounds_lonvalue;
+           for (idx = 0; idx < nbWritten; ++idx)
+             for (int nv = 0; nv < nvertex; ++nv)
+             {
+                writtenBndsLat(nv, idx) = boundslat(nv, int(indexToWrite(idx)));
+                writtenBndsLon(nv, idx) = boundslon(nv, int(indexToWrite(idx)));
+             }
+         }
+         domain->assignMesh(domainName, domain->nvertex);
+         CContextServer* server=CContext::getCurrent()->server ;
+         domain->mesh->createMeshEpsilon(server->intraComm, writtenLon, writtenLat, writtenBndsLon, writtenBndsLat);
+       }
+
      }
    }