13 labelList neighbourCell = neighbour[celli];
16 scalar neighboursEx = 0.0;
18 forAll (neighbourCell, cellj)
20 labelList neighboursNeighbour = neighbour[neighbourCell[cellj]];
21 scalar neighboursNeighbourCellVolumes = 0.0;
23 forAll (neighboursNeighbour, cellk)
25 neighboursNeighbourCellVolumes +=
26 cellVolumes[neighboursNeighbour[cellk]];
30 alphaEx[neighbourCell[cellj]]*cellVolumes[celli]
31 /neighboursNeighbourCellVolumes;
34 alpha[celli] += neighboursEx - alphaEx[celli];