function jdeg=rtod(rad)
rad=rad*180/pi;
deg=fix(rad);
xsbf=(rad-deg)*60;
degf=fix(xsbf);
degm=(xsbf-degf)*60;
jdeg=deg+degf/100.0+degm/10000.0;