mixing cup temperature using posdat.f

Boban

Hi everybody My problem is to calculate an average temperature in a cross section of my channel. It can be done using the command "integrate cset", but in the description of this command is written:"This method is not exact and will provide only a rough estimate of the integrated values". And indeed there are problems in a region where the difference between wall temperature and fluid mean temperature is small, if I want to calculate the heat transfer coefficient (heat flux on the wall has in this case also small value), and I think, that it is due to insufficient accuracy of the values that I get from "integrate cset". I could edit the posdat.f user subroutine to get this mean temperature, but I don´t know, how can I get the coordinates of cells, X, Y and Z in posdat.f. Could someone help me, please? Thank you, Boban