Wcr fixed back (was averaged strangely and wrong)
iph, ith outputs to shock.dat set tom min/max FS velocity (for the complex morphologies)
with iphmin/iphmax, itmin/ithmax