We present a new method for finding principal null directions (PNDs).
Because our method assumes as an input the intrinsic metric and extrin
sic curvature of a space-like hypersurface, we expect it will be usefu
l to numerical relativists; We illustrate our method by finding the PN
Ds of the Kastor-Traschen spacetimes, which contain arbitrarily many Q
= M black holes in a de Sitter background.