Let be a structural model

Let be the regions from to be gridded and populated

Let be the property functions Φ for each L containing an R in

Let be a set of grids GZ

Let be a set of mappings f

Let be a set of arrays a for storing grid-based property values

parallel for each do

Select a property function

PR ¬ Construct the polygonal boundary of R from

Retrieval of GZ and are independent processes

Select Z (subgrid resolution and quality)

(An existing GZ may be reused)

end parallel for