Implementation of parcel method for surface reactions in DSMC