As well as calculating D and v_drift, also calculate errors