function rad = dms2rad( dms )
% rad = dms2rad( [d m s] )

d = dms(:,1);
m = dms(:,2);
s = dms(:,3);

rad = (d + m./60 + s./3600) *pi/180;
