looking at bug in writer

analysis/SubPhase.cpp

@@ -190,6 +190,7 @@ void SubPhase::Basic(){
 				if (Phi(n) != Phi(n)){
 					// check for NaN
 					Phi(n) = 0.0;
+					//printf("Nan at %i %i %i \n",i,j,k);
 				}
 			}
 		}
@@ -344,7 +345,7 @@ void SubPhase::Basic(){
 			}
 		}
 	}
-	printf("wetting interaction = %f, count = %f\n",total_wetting_interaction,count_wetting_interaction);
+	//printf("wetting interaction = %f, count = %f\n",total_wetting_interaction,count_wetting_interaction);
 	total_wetting_interaction_global=Dm->Comm.sumReduce(  total_wetting_interaction);
 	count_wetting_interaction_global=Dm->Comm.sumReduce(  count_wetting_interaction);
 	/* normalize wetting interactions */

analysis/runAnalysis.cpp

@@ -978,8 +978,7 @@ void runAnalysis::basic(int timestep, std::shared_ptr<Database> input_db, SubPha
         /*if (d_regular)
             d_ScaLBL_Comm->RegularLayout(d_Map,Phi,Averages.Phi);
         else */
-        ScaLBL_CopyToHost(Averages.Phi.data(),Phi,N*sizeof(double));
-        
+        	ScaLBL_CopyToHost(Averages.Phi.data(),Phi,N*sizeof(double));
         // copy other variables        
         d_ScaLBL_Comm->RegularLayout(d_Map,Pressure,Averages.Pressure);
         d_ScaLBL_Comm->RegularLayout(d_Map,&Den[0],Averages.Rho_n);

example/Plates/input.db

@@ -7,10 +7,10 @@ Color {
     beta  = 0.95;
     F = 0, 0, 0
     Restart = false
-    timestepMax = 3000
+    timestepMax = 500
     flux = 0.0
     ComponentLabels = -2, -1
-    ComponentAffinity = -1.0, -0.5;
+    ComponentAffinity = 1.0, 1.0;
 }
 
 Domain {
@@ -26,13 +26,14 @@ Domain {
 }
 
 Analysis {
-    blobid_interval = 1000      // Frequency to perform blob identification
-    analysis_interval = 1000    // Frequency to perform analysis
-    restart_interval = 1000    // Frequency to write restart data
-    visualization_interval = 1000        // Frequency to write visualization data
+    analysis_interval = 100    // Frequency to perform analysis
+    visualization_interval = 500
+    subphase_analysis_interval = 100
+    restart_interval = 100000
     restart_file = "Restart"    // Filename to use for restart file (will append rank)
     N_threads    = 4            // Number of threads to use
     load_balance = "independent" // Load balance method to use: "none", "default", "independent"
 }
 
-
+Visualization {
+}

tests/lbpm_color_simulator.cpp

@@ -79,5 +79,3 @@ int main(int argc, char **argv)
 
   Utilities::shutdown();
 }
-
-