$ OMP_NUM_THREADS=128 ./miluph 
Usage ./miluph [options]
                        sph program.
Available options:
-h, --help               This message.
-V, --version            Print out version number.
-v, --verbose            Be talkative (stdout).
-g, --grav               Include self gravitation (default is no).
-a, --theta              Value for the Barnes-Hut theta criterion (allowed values 0 <= theta <= 0.99). No default.  If
not set, pp-gravity is used.
-n, --num                Number of simulation steps.
-I, --integrator         Type of time integration scheme: implemented are euler, rk2, rk2_adaptive, pc, defaults to
rk2_adaptive.
-T, --starttime          Time of simulation start.
-t, --timeperstep        Time for one simulation step.
-F, --firsttimestep      First time step size before adaptive integrator takes over (default 1).
-D, --dt                 Fixed time step size. (default: not used, CFL or adaptive instead, depending on integrator)
-m, --materials          Name of material constants file (default: materialconstants.data).
-s, --scenario           Name of the scenario constants file (default: materialscenario.data).
-f, --filename           Name of input data file (default: disk.0000).
                         Input data file name format is something like 'string'.XXXX, where
                         XXXX means runlevel and zeros.
-o, --threadio           Use extra thread to write output file (needs more mem).
-H, --hdf5_output                Use hdf5 for output (default is FALSE). Implies --threadio.
-X, --hdf5_input                 Use hdf5 for input (default is FALSE).
-Y, --format             Print information about input and output format of the data files.
-R, --restart            Input file is an old miluph output file, which has a different format.
-r, --rdump              Specify dump file name if run is to be continued from a dump file.
                 The values of -t, -D, -g are read from the dump file, but can be
                 overridden on the command line. If -r is specified, the file name
                 given in the -f option ('string'.XXXX) will be interpreted as the output file
                 name with XXXX replaced by the next runlevel from the dump file.
-b, --blevel             Specify 'binary level' (default value is 0):
                 0: do not write binary files
                 1: write binary dump files at the start and after the last simulation step.
                    Backups are kept to ensure a valid dump in case the program crashes.
                 2: write binary dump files at the start and after each simulation step.
                 3: write binary dump files at the start and after each step, suppress text files.
                 If 'string'.XXXX is specified in -f, the binary files will be called
                 'string'.bin.XXXX.start (the start file) and 'string'.bin.XXXX, respectively.

File format is 'x y [z] v_x v_y [v_z] mass [rho] [e] [material_type] [number_of_flaws] [damage] [S_xx [S_xy S_xz]]
[[flaw_1] ... [flaw_number_of_flaws]]'.
File output format is 'x y [z] v_x v_y [v_z] mass [rho] [e] [hsml] [noi] [material_type] [number_of_flaws]
[number_of_activated_flaws] [local_strain] [damage] [p] [S_xx [S_xy S_xz]] [[flaw_1] ... [flaw_number_of_flaws]]'.


authors: Christoph Schaefer & Thomas Maindl.

$ ./miluph --format
Data file format for ./miluph
dimension = 2
./miluph version (hydro or solid): solid

input file format (space seperated values, integer: material type, number of flaws, doubles: all the remaining):
1:x[0] 2:x[1] 3:v[0] 4:v[1] 5:mass 6:density 7:energy 8:material type 9:number of flaws 10:damage 11:S[0] 12:S[1] 
13:S[2] 14:S[3] 15->15+number of flaws:activation thresholds for this particle

output file format: (non-hdf5, for hdf5 use h5ls):
1:x[0] 2:x[1] 3:v[0] 4:v[1] 5:mass 6:density 7:energy 8:smoothing length 9:number of interaction partners 
10:material type 11:number of flaws 12:number of activated flaws 13:local_strain 14:damage 15:pressure 
16:S[0] 17:S[1] 18:S[2] 19:S[3] 20->20+number of flaws:activation thresholds for this particle

#
# This is the master file of material constants for use with miluph
# TIM 2012, 2013
#
#
# DO NOT CHANGE THE DESCRIPTORS!!
#
# Ice
#
MAT = Ice
# Weibull parameters start (cf. Lange, Ahrens & Boslough 1984)
W_M = 9.1
W_K = 1e46
# Weibull parameters end
shear_modulus = 2.8e9
bulk_modulus = 9.47e9
yield_stress = 1e9
# Tillotson EOS parameters
till_rho_0 = 0.917e3
till_A = 9.47e9
till_B = 9.47e9
till_E_0 = 10e6
till_E_iv = 0.773e6
till_E_cv = 3.04e6
till_a = 0.3
till_b = 0.1
till_alpha = 10.0
till_beta = 5.0
# Murnaghan EOS parameters
# from Melosh book
K_0 = 1.502e9
n = 5.23
rho_0 = 0.91e3
#
#
#
# re-normed rubber
MAT = TestRubber
# No Weibull parameters as we do not want the material to break
shear_modulus = 0.22
bulk_modulus = 1.0
yield_stress = 0.3
# # Murnaghan EOS parameters (Melosh refers to Kieffer and Simonds 1980)
K_0 = 1.0
n = 1.0
rho_0 = 1.0
#
# polytropic gas
MAT = Polytropic_Gas
polytropic_K = 1.0
polytropic_gamma = 1.6

#
# These are material-dependent scenario parameters for use with miluph
# TIM 2012
#
# DO NOT CHANGE THE DESCRIPTORS!!
#
MTYPE = 1
MAT = Polytropic_Gas
EOS = P
hsml = 2.01
alpha = 1.0
beta = 2.0