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

xbxu at mcs.anl.gov xbxu at mcs.anl.gov
Fri Jun 29 11:39:04 CDT 2012


Author: xbxu
Date: 2012-06-29 11:39:04 -0500 (Fri, 29 Jun 2012)
New Revision: 5615

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

Modified: MOAB/branches/NCparperiodic/src/io/ReadNC.cpp
===================================================================
--- MOAB/branches/NCparperiodic/src/io/ReadNC.cpp	2012-06-29 16:15:07 UTC (rev 5614)
+++ MOAB/branches/NCparperiodic/src/io/ReadNC.cpp	2012-06-29 16:39:04 UTC (rev 5615)
@@ -1794,11 +1794,19 @@
     ilCVals.resize(lDims[3] - lDims[0] + 1);
     lCDims[3] = lDims[3];
   }
-  else if (!locallyPeriodic[0] && globallyPeriodic[0] && lDims[3] > gDims[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;
+      ilVals.resize(lDims[3] - lDims[0]);
+      ilCVals.resize(lDims[3] - lDims[0]);
+      lCDims[3] = lDims[3]-1;
+    }
+    else
+      {
+	ilVals.resize(lDims[3] - lDims[0] + 1);
+	ilCVals.resize(lDims[3] - lDims[0]);
+	lCDims[3] = lDims[3]-1;
+      }
   }
 
   lCDims[0] = lDims[0];



More information about the moab-dev mailing list