Dear Fedor,
It seems that I am missing something here. I am running the simulation using the command (mpirun -np 10 ./bstep) which will divide the domain into 10 cuboids. but the above loop only iterates through iC=0.
I apologize for the question but I am genuinely confuse.
Best regards,