pro azelp, ra, dec, year, month, day

lat=-22.97166667
long=292.2969444


;######### Time zone ; summer -3, winter -4 ############




ut=(dindgen(100)*24/99.0)
ut1=[0, 3, 6 , 9 , 12 , 15, 18, 21, 24]
;CT2LST, gst1, 0., 0., 0.0, day, month, year
CT2LST,gst, 0.,0,ut1,day, month, year
; converting from Greenwich Sidereal Time to Local Sidereal Time
; the formula is LST = GST + long
;lst1=float(gst1)+long/15.  ; the long. has a "+" because is measured towards East
lst2=float(gst)+long/15
lst2= lst2 mod 24
;lst=lst*2.*!pi/24.
print,lst2

radec2azel,az,el,ra,dec,year,month,day,ut,lat,long

set_plot,'x'
!p.multi=[0,1,2]

plot, ut, el,xrange=[0,24],xstyle=9,yrange=[0,90],ystyle=1,xtitle='!6 UT ',$
 ytitle='!6 Elevation (deg) ';, title=' Az - El for Pampa la Bola '
axis,0,90, xaxis=1,xticks=8,xtickname=string(lst2),xtitle=' LST '

plot, ut, az,xrange=[0,24],xstyle=9,yrange=[-180,180.],ystyle=1,xtitle='!6 UT ',$
ytitle='!6 Azimuth (deg)'
axis,0,180,xaxis=1,xticks=8,xtickname=string(lst2)

!p.multi=0

set_plot,'ps'

device, filename='azelout.ps'

!p.multi=[0,1,2]

plot, ut, el,xrange=[0,24],xstyle=9,yrange=[0,90],ystyle=1,xtitle='!6 UT ',$
 ytitle='!6 Elevation (deg) ';, title=' Az - El for Pampa la Bola '
axis,0,90, xaxis=1,xticks=8,xtickname=string(lst2),xtitle=' LST '

plot, ut, az,xrange=[0,24],xstyle=9,yrange=[-180,180.],ystyle=1,xtitle='!6 UT ',$
ytitle='!6 Azimuth (deg)'
axis,0,180,xaxis=1,xticks=8,xtickname=string(lst2)

!p.multi=0

device,/close

end