RF code incorporating methods of Murphy et al 2010 and  Evans et al 2010
outputs