; Procedure that plots roster pattern for a reduced map.
; Input info is given as a structure.
; Useful: strsec converts seconds after midnight into time series (hh:mm:ss)
;
; Modified to accept stacks of maps, April 4 2001.

pro roster_plot, indata
!P.CHARSIZE=1.5
!P.CHARTHICK=1.8

data=indata
pom=size(data)
if (pom(0) eq 3) then nz=1
if (pom(0) eq 4) then nz=pom(4)
nx=pom(2)
ny=pom(3)

srcname=string(data[0,0,0].h.proc.srcname)

ramax=max(data[0,*,*,*].rahr)
ramin=min(data[0,*,*,*].rahr)
decmin=min(data[0,*,*,*].decdeg)
decmax=max(data[0,*,*,*].decdeg)
;a=strsec((ramax+0.003)*3600.)
;b=strsec((ramin-0.003)*3600.)
plot,data[0,*,1,0].rahr,data[0,*,1,0].decdeg,/nodata,xstyle=1,ystyle=1,$
xrange=[ramax+0.003,ramin-0.003],yrange=[decmin-0.03,decmax+0.03],$
xthick=3,ythick=3,xtitle='Right Ascention (hours)', $
ytitle='Declination (degrees)', title='Roster Pattern for map '+srcname

for k=0,nz-1 do begin
for j=0,ny-1 do oplot,data[*,*,j,k].rahr,data[*,*,j,k].decdeg,psym=7

endfor

end