program v c1=-.14 c2=-17.32 c3=10.06 print 10 10 format(' ra, dec (in degrees), vlsr :',$) read*,ra,dec,vlsr xt=-cosd(dec)*cosd(ra) yt=cosd(dec)*sind(ra) zt=sind(dec) sum=xt*c1+yt*c2+zt*c3 vhelio=vlsr-sum print*,'vhelio= ',vhelio end