* Fixing bug in RT load balancing when the number of grids with

radiation is less than number of processors.

* StarParticleFinalize should be called before setting BCs.

--HG--
branch : week-of-code

src/enzo/EvolveLevel.C

@@ -538,6 +538,11 @@ int EvolveLevel(TopGridData *MetaData, LevelHierarchyEntry *LevelArray[],
 
     }  // end loop over grids
 
+    /* Finalize (accretion, feedback, etc.) star particles */
+
+    StarParticleFinalize(Grids, MetaData, NumberOfGrids, LevelArray,
+			 level, AllStars, TotalStarParticleCountPrevious);
+
     /* For each grid: a) interpolate boundaries from the parent grid.
                       b) copy any overlapping zones from siblings. */
 
@@ -549,11 +554,6 @@ int EvolveLevel(TopGridData *MetaData, LevelHierarchyEntry *LevelArray[],
 			  Exterior, LevelArray[level]);
 #endif
 
-    /* Finalize (accretion, feedback, etc.) star particles */
-
-    StarParticleFinalize(Grids, MetaData, NumberOfGrids, LevelArray,
-			 level, AllStars, TotalStarParticleCountPrevious);
-
     /* If cosmology, then compute grav. potential for output if needed. */
 
     //dcc cut second potential cut: Duplicate?

src/enzo/LoadBalanceHilbertCurve.C

@@ -671,7 +671,7 @@ int LoadBalanceHilbertCurve(grid *GridPointers[], int NumberOfGrids,
     if (WorkLeft == 0) {
       BlockDivisions[i] = grid_num-1;
       ProcessorWork[i] = 0;
-      break;
+      continue;
     }
     WorkThisProcessor = 0;
     GridsThisProcessor = 0;