Added dates and dev_stages vector to the results object; this was needed to analyse the pheno_opt simulations directly; no need to analyse the rmout for those outputs during the optimisation