OilfieldToolsNet/LBPM
4
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

print voxel length when specified

Browse Source
This commit is contained in:
JamesEMcclure
2019-05-01 13:15:53 -04:00
parent cd904c9c91
commit 18d0affaaa
1 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

12
common/Domain.cpp
View File

@@ -118,12 +118,6 @@ void Domain::initialize( std::shared_ptr<Database> db )
d_db = db;
auto nproc = d_db->getVector<int>("nproc");
auto n = d_db->getVector<int>("n");
voxel_length = 1.0;
if (d_db->keyExists( "voxel_length" )){
auto voxel_length = d_db->getScalar<double>("voxel_length");
if (rank==0) printf("voxel length = %f micron \n", voxel_length);
}
if (d_db->keyExists( "InletLayers" )){
auto InletCount = d_db->getVector<int>( "InletLayers" );
@@ -164,6 +158,12 @@ void Domain::initialize( std::shared_ptr<Database> db )
// Fill remaining variables
N = Nx*Ny*Nz;
Volume = nx*ny*nx*nproc[0]*nproc[1]*nproc[2]*1.0;
voxel_length = 1.0;
if (d_db->keyExists( "voxel_length" )){
auto voxel_length = d_db->getScalar<double>("voxel_length");
if (myrank==0) printf("voxel length = %f micron \n", voxel_length);
}
id = new signed char[N];
memset(id,0,N);
BoundaryCondition = d_db->getScalar<int>("BC");