iv = an index running over all vertices from 1 to Nv. In both the vertex coordinates file and in the vertex-vertex
connectivity file the vertices should be ordered by increasing number without skipping. Currently, the code takes at most 999 vertices.
ic = an index running over all cells from 1 to Nc. In both the cell coordinates file and in the cell-vertex
connectivity file the cells should be ordered by increasing number without skipping. Currently, the code takes at most 999 cells.
iv1, iv2, iv3 = the indices of the vertices that connect to vertex iv. If a vertex does not have three neighbours
(eg. a boundary vertex) please fill the index with a 0 (zero). Eg., a line that reads 17 252 69 0
means that vertex number 17 is connected to vertices 252 and 69 but has no more neighbours.
ic1, ic2, ..., ic14 = the indices of the vertices that surround cell ic. If the cell does not have 14 vertices please
fill the index with a 0 (zero). If you have in your system cells surrounded by more than 14 vertices please let me know and
I will modify the code to accommodate you.
At present the code applies only to two dimensional foams. A modification to planar granular assemblies is
straightforward (in the sense that a couple of the subroutines need to be bypassed). If If you are keen to analyse a
structure of the latter please let me know and if there is enough demand to run the code on granular systems I will modify
it.