add plot option prior to adding back 2 axles
add plotfile to function signature

if pass a plotfile argument, then you get amelia diagnostic plots, which
take a while to generate