Input Data Options
	tomographic data (1)
	segmented data (2)
	burn data (3)
	medial axis data (4)
	throat data (5)
	pore-throat network data (6)
	fluid data (7)
	lattice Boltzmann simulation (8)
Enter choice: 4

Data Processing Options
	plot MA (1)
	plot MA with MS (2)
	convert MS files to segmented file type (3)
	histogram and fit pore MA burn data (4)
	embed search tree(s) (5)
	shortest paths (6)
	min cut/max flow (7)
	convert data to/from ascii format (8)
	convert MA list to/from node/path formats (9)
	modify MA (10)
	histogram MA info (11)
	throat calculation (12)
	convert to graph format (13)
	reconstruct burn layers (14)
	MA path curvature and torsion (15)
	maximal balls algorithm (16)
	unpad MA cluster/path files (17)
Enter choice: 12

Throat computation algorithm options
	Dijkstra shortest paths algorithm (1)
	wedge finding algorithm (2)
	subsectional Dijkstra search (3)
Enter choice: 1

    Files assumed to be labelled  basename.ext,
    basename  is limited to 255 characters maximum,
    ext       is a numerical designator lying between 000 and 999
    It is assumed that the files have consecutively numbered extensions
        ie.  000 -> 056,  or   021 -> 049


Enter basename for burn file(s): ../burn/sw
Are file(s) compressed? [y,n]: y


Resize information selection.
NOTE, row and column numbers begin from 0

For each slice,
   enter first and last column of data to use (x-dir): 0 255
   enter first and last row of data to use (y-dir): 0 255

enter first and last slice of data to use (z-dir): 1 256



Enter basename for singlet loc/burn file(s): ../ma_t/cp_sngl

Enter basename for cluster/path loc/burn file(s): ../ma_t/cp_loc

Enter basename for cluster/path structure file(s): ../ma_t/cp_struct
Are file(s) compressed? [y,n]: y


Enter basename for throat file(s): ../throats/sw_dij
Are throat file(s) to be compressed? [y,n]: y

Is the data two dimensional (y,n(dflt): n

Request expert menu? (y,n(dflt)): n

By default 13 (maximum burn number) burn steps are done before
abandoning a search for throats along a medial axis path.
You however have a choice of
	changing this value (c)
	opting for a local criteria based on each path (p)
	leaving it unchanged (u).
Enter choice (c,p,u(dflt)): c

Enter maximum number of burn steps to be done before
abandoning a search for throats along a medial axis path: 40

Average burn number on Medial Axis is 2

By default throat calculation is attempted for paths
which are at least 2 (average burn number) long.
Do you wish to specify a new length (y,n(dflt)): n
Compute throats for
	branch-branch paths[n,y(dft)]: y
	branch-leaf paths[n,y(dft)]: y
	leaf-leaf paths[n,y(dft)]: n

There are 2560 ma paths. Compute throats for all? (y,n): y

Subdivide non monotonic paths to find multiple throats on the
	same path (y(dflt),n): n

In addition to the path length control option specified,
ignore branch-branch and branch-leaf paths if their
length is less than the maximum burn number on the
clusters to which they attach? (n(dflt),y): n