Subgroup method for taking account of the spatial distribution of unscattered and singly scattered neutrons in multigroup shielding calculations