popsub function

Subset data by population