  We have developed the anisotropic heat transport equation for rotating
neutron stars. With a simple model of neutron star, we also model the
propagation of heat pulses result from transient energy releases inside the
star. Even in slow rotation limit, the results with rotational effects involved
could differ significantly from those obtained with a spherically symmetric
metric in the time scale of the thermal afterglow.
