[MOAB-dev] r5612 - MOAB/branches/NCparperiodic/src/io

xbxu at mcs.anl.gov xbxu at mcs.anl.gov
Fri Jun 29 10:57:35 CDT 2012


Author: xbxu
Date: 2012-06-29 10:57:34 -0500 (Fri, 29 Jun 2012)
New Revision: 5612

Modified:
   MOAB/branches/NCparperiodic/src/io/ReadNC.cpp
Log:
set the local dimension for cell centers properly

Modified: MOAB/branches/NCparperiodic/src/io/ReadNC.cpp
===================================================================
--- MOAB/branches/NCparperiodic/src/io/ReadNC.cpp	2012-06-29 15:07:29 UTC (rev 5611)
+++ MOAB/branches/NCparperiodic/src/io/ReadNC.cpp	2012-06-29 15:57:34 UTC (rev 5612)
@@ -1792,14 +1792,19 @@
       // if locally periodic, doesn't matter what global periodicity is, # vertex coords = # elem coords
     ilVals.resize(lDims[3] - lDims[0] + 1);
     ilCVals.resize(lDims[3] - lDims[0] + 1);
+    lCDims[3] = lDims[3];
   }
   else if (!locallyPeriodic[0] && globallyPeriodic[0] && lDims[3] > gDims[3]) {
       // globally periodic and I'm the last proc, get fewer vertex coords than vertices in i
     ilVals.resize(lDims[3] - lDims[0]);
     ilCVals.resize(lDims[3] - lDims[0]);
+    lCDims[3] = lDims[3]-1;
   }
-  
 
+  lCDims[0] = lDims[0];
+  lCDims[4] = lDims[4]-1;
+  lCDims[1] = lDims[1];
+
   if (-1 != lDims[1]) {
     jlVals.resize(lDims[4] - lDims[1] + 1);
     jlCVals.resize(lCDims[4] - lCDims[1] + 1);



More information about the moab-dev mailing list