Fortran for (RA, DEC, Vlsr) to Vhelio

	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