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: 8

Data Processing Options
	LBGK basic compressible model (1)
	LBGK incompressible model computing pressure (2)
	LBGK incompressible model (3)
	LBGK linearized basic model (4)
	MRT multiple-relaxation-time model (5)
	2D LBM (6)
Enter choice: 1

Options
	1) Poiseulle flow in an axis aligned tube
	2) Poiseulle flow in an inclined tube
	3) Hyperboloid flow
	4) Flow in general porous rock volume
	5) Flow in a specific pore body
	6) Flow through a throat connecting two pores
Enter option: 6

Solid wall boundary condition implementations:
	1) Standard bounce-back(dflt)
	2) Half-way wall Bounce-back(straight walls)
	3) Curved boundary(Mei/Shyy/Luo 2000.)
Enter option: 3

Are initial conditions stored in a file?(y/n): n
Store final iteration state to a file?(y/n): n

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

Enter first and last slice of data to use: 1 256
Slice size of data, (nx,ny): 256 256

Enter throat data file name: ../../sw/throats/sw_dij_cln
Are throat data file(s) compressed? [y,n]: y

Enter basename for nodal pore/throat network (*.np2th) file(s): ../../sw/throats/sw_dij_cln
Are file(s) compressed? [y,n]: y


Ignore throats not used in pore/throat network?(y(dflt),n): n


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

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

1135 relevant throats used out of total 1135

Num throat voxels that are - grain 0 shared 75

	Throat indices start from 1 and end with 1135
Enter range of throats through which you wish to simulate flow(start,end): 10 10

Produce rasterplot of default direction velocity? (y,n(dlft)): n

Plot velocity in Geomview(3D)?(y/n(dflt)): y
Enter velocity plot file basename: ../plots/veloc
	Rainbow coloring of velocity direction vectors shows velocity vector 
	length - red = the smallest through violet = the largest

	Computation is done in lattice units. To study flow qualitatively
	there is no need to input physical information.
Use physical units? (y/n(dflt)): y
Enter voxel length in microns(1e-6m): 4.93

	Time step in physical units will be computed from sound speed.
Enter sound speed in the fluid in m/s: 1482

time step in phys. units 1.92061e-09s lattice speed 2566.9m/s

Enter fluid density  in kg*m^(-3): 1000
Reference mass chosen as 1.19823e-13
	Knudsen number depends on temperature, pressure, diameter
	of the fluid molecules and the voxel length.
Check Knudsen number for the flow? (y,n(dflt)): n

	Physical LBE relaxation parameter is computed from fluid viscosity.
Compute physical relaxation parameter?(y,n(dflt)): n
Enter pressure gradient in physical units(Pa/m): 3.35466e+06
Enter relaxation parameter for LBE: 0.750000

	Resolution can be refined, i.e.each voxel can be replaced by 
	M x M x M voxels.
Enter M (enter 1 for no refining): 1

Enter tolerance for relative velocity error at each node: 1e-06

	Iterations for slow-converging runs should be limited.
Enter max. number of iterations allowed: 5000