pro skymap
readcol,'coords',name,r1,r2,r3,d1,d2,d3,format='a,f,f,f,f,f,f'

r=r1+r2/60.+r3/3600.
d=d1+d2/60.+d3/3600.

rt=r/24.*360.

euler,rt,d,a,b,1
print,b

n1=10.
n2=1.
n3=.5
gx1=indgen(n1*360)/n1
gx2=indgen(n2*360)/n2
gx3=indgen(n3*360)/n3
gy1=0.
gy2=20.
gy3=40.
gy4=-20.
gy5=-40.
gy6=60.
gy7=-60.

euler,gx1,gy1,rg1,dg1,2
euler,gx2,gy2,rg2,dg2,2
euler,gx2,gy3,rg3,dg3,2
euler,gx2,gy4,rg4,dg4,2
euler,gx2,gy5,rg5,dg5,2
euler,gx3,gy6,rg6,dg6,2
euler,gx3,gy7,rg7,dg7,2
euler,0,0,cx,cy,2

rg1=rg1/360.*24.
rg2=rg2/360.*24.
rg3=rg3/360.*24.
rg4=rg4/360.*24.
rg5=rg5/360.*24.
rg6=rg6/360.*24.
rg7=rg7/360.*24.
cx=cx/360.*24.

plot,rg1,dg1,psym=3,xrange=[24,0],xstyle=1,yrange=[-90,90],ystyle=1,xtitle='RA (hours)',ytitle='Dec (degrees)'
oplot,rg2,dg2,psym=3
oplot,rg3,dg3,psym=3
oplot,rg4,dg4,psym=3
oplot,rg5,dg5,psym=3
oplot,rg6,dg6,psym=3
oplot,rg7,dg7,psym=3
oplot,[cx,cx],[cy,cy],psym=1

oplot,r,d,psym=2
;read,temp

;plot,a,b,psym=2


stop
end