function jrad=deg2rad(deg)
degd=fix(deg);
degf=(deg-degd)*100;
degm=fix(degf);
degs=(degf-degm)/36.0;
jrad=(degd+degm/60.00+degs)* pi / 180.0;