function ang = angdiffLocal(alpha, beta)

   ang = rem(  (beta-alpha), 2*pi  );
   
   if ( ang >= pi ) ang = ang - 2*pi; end

   if ( ang < -pi ) ang = ang + 2*pi; end

end   

