Calculations of accreting magnetized neutron star atmospheres heated by the gradual deceleration of protons via Coulomb collisions are presented. Self consistent determinations of the temperature and density structure for different accretion rates are made by assuming hydrostatic equilibrium and energy balance, coupled with radiative transfer. The full radiative transfer in two polarizations, using magnetic cross sections but with cyclotron resonance effects treated approximately, is carried out in the inhomogeneous atmospheres.