### Pretvorba GLON,GLAT v RA,DEC - S53MV 31.10.2016 ###

import math

def gllvrd(glon,glat):				#galakticne koordinate v rektascenzijo,deklinacijo
	glon=glon-32.9				#galakticna dolzina dviznega vozla galaksije 32.9
	glon=math.radians(glon)			#pretvorba v radiane
	glat=math.radians(glat)
	x=math.cos(glon)*math.cos(glat)		#pretvorba v kartezicne koordinate
	y=math.sin(glon)*math.cos(glat)
	z=math.sin(glat)
	nagib=math.radians(62.8)		#nagib ravnine galaksije 62.8
	yy=y*math.cos(nagib)-z*math.sin(nagib)
	z=y*math.sin(nagib)+z*math.cos(nagib)
	ra=math.degrees(math.atan2(yy,x))	#pretvorba v krogelne koordinate
	dec=math.degrees(math.asin(z))
	ra=ra+282.9				#rektascenzija dviznega vozla galaksije 282.9
	if ra>360:				#omeji obmocje rektascenzije 0...360
		ra=ra-360
	return ra,dec

while True:
	glon=input('Vstavi galakticno dolzino: ')
	glat=input('Vstavi galakticno sirino: ')

	ra,dec=gllvrd(glon,glat)		#galakticne koordinate v rektascenzijo,deklinacijo

	print 'Rektascenzija:',ra,'Deklinacija',dec
	print 'RA:',int(ra/15.0),'h',(ra/15.0-int(ra/15.0))*60 

#konec