I have read your response on setting the value of rho to each cell. That is exactly what I need and I am thankful for your help but it seems that this function does not work for parallel running. Does it? The value of iC is always zero and does not increase further. Any help on that would really be appreciated.