“…where S , the solar flux at 1 AU, is equal to 1370 W m −2 , A H is the Bond albedo, taken to be 0.0108, from a standard geometric value of 0.04 and phase integral value of 0.27 (Lamy et al 2007), is the infrared emissivity, taken to be 0.95 (Birkebak 1972), σ is the Stefan-Boltzmann constant, r h is the heliocentric distance, and i is the angle of incidence of the sun to the local surface. As a first step, the illumination geometry and projected shadows are calculated for each facet of the shape model over one complete nucleus rotation.…”