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)
Enter choice: 6

Data Processing Options
	individual pore surface visualization/area (1)
	throat/body statistics (2)
	pore principal diameters (3)
	throat principal diameters (4)
	pore surface curvature (under constr.) (5)
	pore/throat subset info (6)
Enter choice: 2

This routine will histogram the np2th and th2np information contained
	in either single or multiple pairs of files.

Enter number of pairs of files: 1
Enter np2th and th2np file basename: ../throats/sw_dij_cln
Are np2th and th2np info files compressed? [y,n]: y

Ignore nodal pores that touch the exterior? (y,n(dflt)): y

Enter voxel length in physical units: 4.1
Enter basename for throat/nodal pore statistics plots files: ../results/np_thr_net
Enter title for plots: Berea sample

Set log axis for nodal pore volume plots
	x (y,n(dflt)): y
	y (y,n(dflt)): n

Set log axis for effective pore radii plots
	x (y,n(dflt)): n
	y (y,n(dflt)): n

Set binning for nodal pore volume histogram
	number of bins (c)
	bin width      (w)
Enter choice: w

Specify bin width to be used: 0.2

Set binning for effective pore radii histogram
	number of bins (c)
	bin width      (w)
Enter choice: w

Specify bin width to be used: 2.5

Set log axis for throat surface area plots
	x (y,n(dflt)): y
	y (y,n(dflt)): y

Set log axis for effective throat radii plots
	x (y,n(dflt)): n
	y (y,n(dflt)): n

Set binning for throat surface area histogram
	number of bins (c)
	bin width      (w)
Enter choice: w
Specify bin width to be used: 0.08

Set binning for effective throat radii histogram
	number of bins (c)
	bin width      (w)
Enter choice: w
Specify bin width to be used: 2.5

Set binning for thr/np size ratio histogram
	number of bins (c)
	bin width      (w)
Enter choice: w

Specify bin width to be used: 0.05