nope -- velocity BC still goofy trying again

This commit is contained in:
James E McClure 2015-07-23 15:44:59 -04:00
parent 90cbc71a45
commit 7e572d4936

View File

@ -447,6 +447,11 @@ int main(int argc, char **argv)
for (i=0; i<Dm.Nx*Dm.Ny*Dm.Nz; i++) Dm.id[i] = id[i];
Dm.CommInit(MPI_COMM_WORLD);
//...........................................................................
if (rank==0) printf ("Create ScaLBL_Communicator \n");
// Create a communicator for the device
ScaLBL_Communicator ScaLBL_Comm(Dm);
// set reservoirs
if (BoundaryCondition > 0){
for ( k=0;k<Nz;k++){
@ -463,11 +468,6 @@ int main(int argc, char **argv)
}
}
//...........................................................................
if (rank==0) printf ("Create ScaLBL_Communicator \n");
// Create a communicator for the device
ScaLBL_Communicator ScaLBL_Comm(Dm);
//...........device phase ID.................................................
if (rank==0) printf ("Copy phase ID to device \n");
char *ID;