move error calculation to data dir
script to make kathryn's suggested error plots
separate out proper motion plot
return a quantity object
use functions from streams.data.gaia