A calculational technique is presented for predicting thermal neutron fluxes in the primary shields of reactor systems which eliminates reliance on mock-up experimental data. A multigroup P1 approach is employed with the spatial dependence of the neutron attenuation adjusted through use of a point source attenuation kernel for a homogeneous hydrogenous medium. Comparison of calculation with experiment is presented.