add output
add 'modelfit' output and adjust the default for 'dontsave'; add 'dosave' input; add totalnumvxs output; allow <dontsave> to refer to the auxiliary variables
new function