Fortran for (l, b, Vlsr) to Vhelio

	program v
	c1=-10.27
	c2=15.32
	c3=7.74
	print 10
10	format(' l, b, vlsr :',$)
	read*,xl,b,vlsr
	ut=-cosd(b)*cosd(xl)
	vt=cosd(b)*sind(xl)
	wt=sind(b)
	sum=ut*c1+vt*c2+wt*c3
	vhelio=vlsr-sum
	print*,'vhelio= ',vhelio
	end