Add calculation of planetary albedo