IDL routine for Gagné et al. 2014b, Figure 1