1: function
2: Array containing property values
3:
4: parallel for each node or cell in GZ do
5: x ¬ Either grid node or cell barycentre
6: Map point in R to point in D
7:
8: Interpolate Φ in x'
9:
10:
11: end parallel for
12: return a
13: end function