diff --git a/tests/lbpm_BGK_simulator.cpp b/tests/lbpm_BGK_simulator.cpp
index 9053ada9..f16e7d93 100644
--- a/tests/lbpm_BGK_simulator.cpp
+++ b/tests/lbpm_BGK_simulator.cpp
@@ -275,7 +275,7 @@ int main(int argc, char **argv)
 		ScaLBL_Communicator ScaLBL_Comm(Mask);
 		int *neighborList;
 		IntArray Map(Nx,Ny,Nz);
-		neighborList= new int[18*(Np)];
+		neighborList= new int[18*Npad];
 		Np = ScaLBL_Comm.MemoryOptimizedLayoutAA(Map,neighborList,Mask.id,Np);
 		MPI_Barrier(comm);
 		
diff --git a/tests/lbpm_permeability_simulator.cpp b/tests/lbpm_permeability_simulator.cpp
index 10f62309..1a9f5418 100644
--- a/tests/lbpm_permeability_simulator.cpp
+++ b/tests/lbpm_permeability_simulator.cpp
@@ -256,8 +256,7 @@ int main(int argc, char **argv)
 		MPI_Allreduce(&sum_local,&sum,1,MPI_DOUBLE,MPI_SUM,comm);
 		porosity = sum*iVol_global;
 		if (rank==0) printf("Media porosity = %f \n",porosity);
-		//int Npad=(Np/32 + 1)*32;
-		//Np=Npad;
+
 		//.........................................................
 		// don't perform computations at the eight corners
 		id[0] = id[Nx-1] = id[(Ny-1)*Nx] = id[(Ny-1)*Nx + Nx-1] = 0;