Update harmonize.R