Added params to the list of things passed to the vrm.f90 subroutine (for use of an initial strain rate with plastic materials). Modified to use either an input penalty value or bulk viscosity