Added rate-dependence effects (Duvaut-Lions formulation, see Simo & Hughes, p. 217) to model. .in files should add a rate-dependence parameter just before debug parameter. This should be 0 to recover the rate-independent version of the model.
Also cleaned up and added some checks to reading and writing of input parameters.
