  We investigate the influence of different magnetic field configurations on
the temperature distribution in neutron star crusts. We consider axisymmetric
dipolar fields which are either restricted to the stellar crust, ``crustal
fields'', or allowed to penetrate the core, ``core fields''. By integrating the
two-dimensional heat transport equation in the crust, taking into account the
classical (Larmor) anisotropy of the heat conductivity, we obtain the crustal
temperature distribution, assuming an isothermal core. Including quantum
magnetic field effects in the envelope as a boundary condition, we deduce the
corresponding surface temperature distributions. We find that core fields
result in practically isothermal crusts unless the surface field strength is
well above $10^{15}$ G while for crustal fields with surface strength above a
few times $10^{12}$ G significant deviations from isothermality occur at core
temperatures inferior or equal to $10^8$ K. At the stellar surface, the cold
equatorial region produced by the quantum suppression of heat transport
perpendicular to the field in the envelope, present for both core and crustal
fields, is significantly extended by the classical suppression at higher
densities in the case of crustal fields. This can result, for crustal fields,
in two small warm polar regions which will have observational consequences: the
neutron star has a small effective thermally emitting area and the X-ray pulse
profiles are expected to have a distinctively different shape compared to the
case of a neutron star with a core field. These features, when compared with
X-ray data on thermal emission of young cooling neutron stars, will open a way
to provide observational evidence in favor, or against, the two radically
different configurations of crustal or core magnetic fields.
