Dear Fedor,
What I want to do is to avoid using (getBlock). So I am thinking to obtain concrete cell number using getCellId(iX,iY)) and then use Cell (BlockLattice< T, TDESCRIPTOR > &block, std::size_t CellId) to set rho using [olb::Cell< T, TDESCRIPTOR >::defineRho(scalar)]. All of these are in the setTemperature function that I have written. However, it seems that I am missing something. What is the BlockLattice here?
Your help is really needed and appreciated.
Thanks