1: Fix a given dipolar + quadrupolar configuration and compute consistently the thermal map of the surface
 We can fix 7 parameters and ³see² the rotation of the thermal surface:
                  b quadi = Bquadi /Bdip      i=0,Š4
                  c = angle between LOS and spin axis
                  x = angle between magnetic and spin axis