This works fine in a single core run, but when I do a parallel run, the data extraction fails. I have read the forum and someone mentioned the formation of the (Temp_Fluid1) array and having the first element of the array as a cuboid rank. Is it due to that? if so, how can I obtain the cuboid rank using certain lattice coordinates? and then implement that digit into position? I hope the problem is well identified.
you better work with global coordinates, like in the examples cylinder2d and cylinder3d it is done for the pressure drop. You can get the cubpoid number by methods provided by the cuboidGeometry2d class. Note that you have to differ betrween local and global coordinates as well as local and global cuboidNumbers.