Fix.

python/UnstructuredMesh.py

@@ -1025,8 +1025,13 @@ def _computeOFConn(self):
         return  faceSizes, faceConn, pts
 
     def _readOFGrid(self, caseDir):
-        """Read in the mesh points and connectivity from the polyMesh
-        directory
+        """
+        Read in the mesh points and connectivity from the polyMesh
+        directory.
+
+        NOTE: To ensure a consistant starting point in DAFoam. One needs to copy 
+        points_orig to points in pyDAFoam before initializing USMesh. This function 
+        will no longer perform this copying!
 
         Parameters
         ----------
@@ -1042,12 +1047,6 @@ def _readOFGrid(self, caseDir):
         # Copy the reference points file to points to ensure
         # consistant starting point
         self.OFData = ofm.getFileNames(caseDir,comm=self.comm)
-        # NOTE: comment out the copying of points_orig to points here because it should 
-        # be done in the pyDAFoam layer
-        #try:
-        #    shutil.copyfile(self.OFData['refPointsFile'], self.OFData['pointsFile'])
-        #except:
-        #    raise Error('USMesh: Unable to copy %s to %s.'%(self.OFData['refPointsFile'], self.OFData['pointsFile']))
 
         # Read in the volume points
         self.OFData['x0'] = ofm.readVolumeMeshPoints()