Ben's script + plots