ps_spells function

Create "spells" by cross-sectional unit, even more generally