added back if(nrow(mutated.populations > 0) { total.mutations = total.mutations + 1, which had been combined into subsequent if statements and I think was the source of error for why species were not diversifying in latest version of code
replacing a lost % in %dopar%
deleted print statement in the extinct.pops.output.times loop