SPH-Hydrodynamic Accretion and Remnant Debris (SHARD)

This repository contains a Python driver for SHARD: a collision-resolving late-stage planet-formation workflow built on REBOUND/mercurius. The script detects N-body collisions, maps non-special embryo--embryo impacts onto an SPH outcome catalogue, interpolates the largest remnants, reconstructs small debris with mass-weighted KMeans clustering in velocity space, and writes collision/restart diagnostics.