def get_atherr(dset):      
      """Returns array with uncertainty in ash top height from VOLCAT."""
      """Default units are km above sea-level"""
      height_err = dset.ash_cth_uncertainty
      height_err = height_err.where(height_err != height_err._FillValue, drop=True)
      return height_err