This program propagates muons in Ice to/through the detector
Available options are: -length=[LENGTH of the detector volume in meters]
                       -radius=[RADIUS of the detector volume in meters]
                       -vcut=[value of vcut used for the 1st region]
                       -ecut=[ecut in  MeV  used for the 2nd region]
                       -medi=[medium name]
                       -mediadef=[file with media definitions]
                       -tau  propagate taus instead of muons
                       -sdec enable stopped muon decay treatment
                       -user     enable the mmc_en user line
                       -user=[z] same, but record energy at z, not CPD
                       -rdmc     enforce compliance with rdmc
                       -amasim   turn on workarounds for amasim
                       -time precise time of flight calculation
                       -cont enable continuous loss randomization
                       -scat enable Moliere scattering
                       -lpm  enable lpm treatment
                       -allm enable allm photonuclear parametrization
                       -phnu=[1-4] choose another photonuclear formula
                       -elow=[muon energy in GeV below which it is lost]
                       -ebig=[upper bound in GeV of the paramet. tables]
                       -surf=[h in meters]  propagate to the plane z=[h]
                       -face  only if detector is on opposite side of it
                       -maxbuf=[max number of lines per event] (10000)
                       -romb=[number of interpolation points]
                       -seed=[integer] sets random number generator seed
                       -raw  save tables in raw format
                       -tdir=[dir] specify directory for paramet. tables