From 5f85b767d6f5fd01463506d832d44f6be02b2ba5 Mon Sep 17 00:00:00 2001
From: Rex Zhe Li <zhe.rex.li@gmail.com>
Date: Mon, 9 Dec 2019 22:30:36 -0500
Subject: [PATCH] add debugging for greyscale lbm

---
 models/GreyscaleModel.cpp          | 56 +++++++++++++++---------------
 tests/lbpm_greyscale_simulator.cpp |  1 +
 2 files changed, 29 insertions(+), 28 deletions(-)

diff --git a/models/GreyscaleModel.cpp b/models/GreyscaleModel.cpp
index 5af10205..2ec3b85e 100644
--- a/models/GreyscaleModel.cpp
+++ b/models/GreyscaleModel.cpp
@@ -643,34 +643,36 @@ void ScaLBL_GreyscaleModel::VelocityField(){
 
 void ScaLBL_GreyscaleModel::WriteDebug(){
 	// Copy back final phase indicator field and convert to regular layout
-/*	ScaLBL_CopyToHost(Porosity.data(), Poros, sizeof(double)*N);
+	DoubleArray PhaseField(Nx,Ny,Nz);
 
-	FILE *OUTFILE;
-	sprintf(LocalRankFilename,"Phase.%05i.raw",rank);
-	OUTFILE = fopen(LocalRankFilename,"wb");
-	fwrite(PhaseField.data(),8,N,OUTFILE);
-	fclose(OUTFILE);
+	//ScaLBL_CopyToHost(Porosity.data(), Poros, sizeof(double)*N);
 
-    ScaLBL_Comm->RegularLayout(Map,&Den[0],PhaseField);
-	FILE *AFILE;
-	sprintf(LocalRankFilename,"A.%05i.raw",rank);
-	AFILE = fopen(LocalRankFilename,"wb");
-	fwrite(PhaseField.data(),8,N,AFILE);
-	fclose(AFILE);
-
-	ScaLBL_Comm->RegularLayout(Map,&Den[Np],PhaseField);
-	FILE *BFILE;
-	sprintf(LocalRankFilename,"B.%05i.raw",rank);
-	BFILE = fopen(LocalRankFilename,"wb");
-	fwrite(PhaseField.data(),8,N,BFILE);
-	fclose(BFILE);
-
-	ScaLBL_Comm->RegularLayout(Map,Pressure,PhaseField);
-	FILE *PFILE;
-	sprintf(LocalRankFilename,"Pressure.%05i.raw",rank);
-	PFILE = fopen(LocalRankFilename,"wb");
-	fwrite(PhaseField.data(),8,N,PFILE);
-	fclose(PFILE);
+//	FILE *OUTFILE;
+//	sprintf(LocalRankFilename,"Phase.%05i.raw",rank);
+//	OUTFILE = fopen(LocalRankFilename,"wb");
+//	fwrite(PhaseField.data(),8,N,OUTFILE);
+//	fclose(OUTFILE);
+//
+//    ScaLBL_Comm->RegularLayout(Map,&Den[0],PhaseField);
+//	FILE *AFILE;
+//	sprintf(LocalRankFilename,"A.%05i.raw",rank);
+//	AFILE = fopen(LocalRankFilename,"wb");
+//	fwrite(PhaseField.data(),8,N,AFILE);
+//	fclose(AFILE);
+//
+//	ScaLBL_Comm->RegularLayout(Map,&Den[Np],PhaseField);
+//	FILE *BFILE;
+//	sprintf(LocalRankFilename,"B.%05i.raw",rank);
+//	BFILE = fopen(LocalRankFilename,"wb");
+//	fwrite(PhaseField.data(),8,N,BFILE);
+//	fclose(BFILE);
+//
+//	ScaLBL_Comm->RegularLayout(Map,Pressure,PhaseField);
+//	FILE *PFILE;
+//	sprintf(LocalRankFilename,"Pressure.%05i.raw",rank);
+//	PFILE = fopen(LocalRankFilename,"wb");
+//	fwrite(PhaseField.data(),8,N,PFILE);
+//	fclose(PFILE);
 
 	ScaLBL_Comm->RegularLayout(Map,&Velocity[0],PhaseField);
 	FILE *VELX_FILE;
@@ -693,7 +695,5 @@ void ScaLBL_GreyscaleModel::WriteDebug(){
 	fwrite(PhaseField.data(),8,N,VELZ_FILE);
 	fclose(VELZ_FILE);
 
- * 
- */
 
 }
diff --git a/tests/lbpm_greyscale_simulator.cpp b/tests/lbpm_greyscale_simulator.cpp
index 9f910a32..61322d6d 100644
--- a/tests/lbpm_greyscale_simulator.cpp
+++ b/tests/lbpm_greyscale_simulator.cpp
@@ -55,6 +55,7 @@ int main(int argc, char **argv)
 		Greyscale.Initialize();   // initializing the model will set initial conditions for variables
 		Greyscale.Run();	 
 		Greyscale.VelocityField();
+		Greyscale.WriteDebug();
 	}
 	// ****************************************************
 	MPI_Barrier(comm);