diff --git a/common/Domain.cpp b/common/Domain.cpp index 32e13501..e355310f 100644 --- a/common/Domain.cpp +++ b/common/Domain.cpp @@ -390,14 +390,11 @@ void Domain::Decomp( const std::string& Filename ) } } } - if (RANK==0){ - for (size_t idx=0; idx 0){ // use checkerboard pattern @@ -516,31 +513,31 @@ void Domain::Decomp( const std::string& Filename ) } } } - } - else { - if (inlet_layers_z > 0){ - printf("Mixed reflection pattern at z inlet for %i layers, saturated with phase label=%i \n",inlet_layers_z,inlet_layers_phase); - for (int k = zStart; k < zStart+inlet_layers_z; k++){ - for (int j = 0; j 0){ - SegData[k*global_Nx*global_Ny+j*global_Nx+i] = reflection_id; + else { + if (inlet_layers_z > 0){ + printf("Mixed reflection pattern at z inlet for %i layers, saturated with phase label=%i \n",inlet_layers_z,inlet_layers_phase); + for (int k = zStart; k < zStart+inlet_layers_z; k++){ + for (int j = 0; j 0){ + SegData[k*global_Nx*global_Ny+j*global_Nx+i] = reflection_id; + } } } } } - } - if (outlet_layers_z > 0){ - printf("Mixed reflection pattern at z outlet for %i layers, saturated with phase label=%i \n",outlet_layers_z,outlet_layers_phase); - for (int k = zStart + nz*nprocz - outlet_layers_z; k < zStart + nz*nprocz; k++){ - for (int j = 0; j 0){ - SegData[k*global_Nx*global_Ny+j*global_Nx+i] = reflection_id; + if (outlet_layers_z > 0){ + printf("Mixed reflection pattern at z outlet for %i layers, saturated with phase label=%i \n",outlet_layers_z,outlet_layers_phase); + for (int k = zStart + nz*nprocz - outlet_layers_z; k < zStart + nz*nprocz; k++){ + for (int j = 0; j 0){ + SegData[k*global_Nx*global_Ny+j*global_Nx+i] = reflection_id; + } } } }