fracdeg2ast

 [ra,dec]=fracdeg2ast([ra],[dec])

 Translate simple fractional degrees.
 to the stupid ra=hh:mm:ss dec=dd:mm:ss
 format which astronomers often use.

 Use strings to match ast2fracdeg function
 which is forced to by +/- zero degrees dec
 problem.

0001 function [rao,deco]=fracdeg2ast(ra,dec)
0002 % [ra,dec]=fracdeg2ast([ra],[dec])
0003 %
0004 % Translate simple fractional degrees.
0005 % to the stupid ra=hh:mm:ss dec=dd:mm:ss
0006 % format which astronomers often use.
0007 %
0008 % Use strings to match ast2fracdeg function
0009 % which is forced to by +/- zero degrees dec
0010 % problem.
0011 
0012 ra=ra/15;
0013 raoh=fix(ra);
0014 ra=60*abs(ra-raoh);
0015 raom=floor(ra);
0016 raos=60*(ra-raom);
0017 
0018 decod=fix(dec);
0019 dec=60*abs(dec-decod);
0020 decom=floor(dec);
0021 decos=60*(dec-decom);
0022 
0023 for i=1:length(ra)
0024   rao(i,:)=sprintf('%02.0f:%02.0f:%08.5f',raoh(i),raom(i),raos(i));
0025   deco(i,:)=sprintf('%+03.0f:%02.0f:%07.4f',decod(i),decom(i),decos(i));
0026 end