fixed contact rate calculation, game stopifnot commands to keep effect of infection on death rate positivemade r, d, beta_d, c, and phi arguments to mpi_f