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 limitation of the is the format in which it reads the files: either 999 or 9999 vertices (see comment below).
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.
In the file cell-vertex.dat, iv1, iv2, ..., iv14 = 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 let me know and I will modify the code to accommodate you.
Comment: to input 9999 vertices, reduce the number of neighbours in the cell-vertex.dat file to at most 13 and change the line format to (14(I4,1x)).
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