import time, math	#uporabljene funkcije
pomlad=0.275	#poenostavljeno pomladisce za Julijanski koledar

maxracun=9999	#max stevilo racunov
stracun=0	#stevec racunov
casracun=5	#interval racunanja (sekunde)

sirina=45.948	#zemljepisna sirina opazovalca v stopinjah (+sever)
dolzina=13.634	#zemljepisna dolzina opazovalca v stopinjah (+vzhod)

cs=math.cos(math.radians(sirina))	#smerni vektorji opazovalca
ss=math.sin(math.radians(sirina))
cd=math.cos(math.radians(dolzina))
sd=math.sin(math.radians(dolzina))
ovzhod=[-sd,cd,0]
osever=[-ss*cd,-ss*sd,cs]
ogor=[cs*cd,cs*sd,ss]

def skalar(a,b):
	return (a[0]*b[0])+(a[1]*b[1])+(a[2]*b[2])	#skalarni produkt dveh vektorjev

#zvezda Sirius
rektascenzija=101.284	#rektascenzija v stopinjah
deklinacija=-16.743	#deklinacija v stopinjah

while stracun<maxracun:		#zanka racunanja

	t=time.gmtime()		#precitaj datum in uro realnega casa racunalnika (UTC)
	leto=t[0]
	mesec=t[1]
	dan=t[2]
	ura=t[3]
	minuta=t[4]
	sekunda=t[5]

	print str(leto)+'-'+str(mesec)+'-'+str(dan),		#izpisi datum
	print str(ura)+':'+str(minuta)+':'+str(sekunda),	#izpisi uro

	ddan=(((sekunda/60.0)+minuta)/60.0+ura)/24.0		#decimalni del dneva

	pleto=leto-int((leto/4)*4)				#dan v stiriletju
	sdan=dan
	dmesec=[31,28,31,30,31,30,31,31,30,31,30,31]	#dodaj stevilo dni v mesecih
	n=1	
	while n<mesec:
		sdan=sdan+dmesec[n-1]
		if pleto==0:				#popravek za prestopno leto
			if n==2:
				sdan=sdan+1
		n=n+1
	dleto=[366,365,365,365]				#dodaj stevilo dni v letih
	n=0
	while n<pleto:
		sdan=sdan+dleto[n]
		n=n+1

	zdolzina=(ddan*(1+(1/365.25)))+(sdan/365.25)		#izracunaj zemljepisno dolzino
	zdolzina=zdolzina+pomlad				#dodaj zacetni polozaj pomladisca
	zdolzina=(1+int(zdolzina)-zdolzina)*360.0		#rezultat v stopinje (+vzhod)

	zdolzina=zdolzina+rektascenzija				#izracunaj smerni vektor zvezde
	cr=math.cos(math.radians(zdolzina))
	sr=math.sin(math.radians(zdolzina))
	cd=math.cos(math.radians(deklinacija))
	sd=math.sin(math.radians(deklinacija))
	zvezda=[cd*cr,cd*sr,sd]

	azimut=math.degrees(math.atan2(skalar(zvezda,ovzhod),skalar(zvezda,osever)))
	if azimut<0:
		azimut=azimut+360
	elevacija=math.degrees(math.asin(skalar(zvezda,ogor)))

	print 'Azimut %5.1f'%(azimut), 'Elevacija %4.1f'%(elevacija)

	stracun=stracun+1	#zakljucek zanke racunanja
	time.sleep(casracun)