Fix OMP bug in AddParticles

1 files changed, 9 insertions, 3 deletions

diff --git a/Source/PhysicalParticleContainer.cpp b/Source/PhysicalParticleContainer.cpp
index e8196eedc..c221ddd84 100644
--- a/Source/PhysicalParticleContainer.cpp
+++ b/Source/PhysicalParticleContainer.cpp
@@ -103,12 +103,18 @@ PhysicalParticleContainer::AddParticles (int lev, Box part_box)
 #elif BL_SPACEDIM==2
     scale_fac = dx[0]*dx[2]/num_ppc;
 #endif
-
+    
+    // First touch all tiles in the map in serial
+    for (MFIter mfi = MakeMFIter(lev); mfi.isValid(); ++mfi) {
+        const int grid_id = mfi.index();
+        const int tile_id = mfi.LocalTileIndex();        
+        GetParticles(lev)[std::make_pair(grid_id, tile_id)];
+    }
+    
 #ifdef _OPENMP
 #pragma omp parallel
 #endif
-    {
-
+    {        
         std::array<Real,PIdx::nattribs> attribs;
         attribs.fill(0.0);