pro time_line, guide=guide
;
;
;Day light savings is last sunday in Apr through last sunday in Oct
;
clearplot
nhours = 36
n = 7
plot, fltarr(2), fltarr(2), xrange=[0,nhours],yrange=[0,26], /nodata, xstyle=1+4+8, ystyle=1+4+8
t = anytim2ints('2-jan-93', off=findgen(n)/(n-1)*nhours*60.*60)
xtickv = findgen(n)/(n-1)*nhours			;tickv and tickname is one greater than ticks
day0 = gt_day('1-jan-93')
ref_day= ['Sat', 'Sun', 'Mon']
off_arr = [0,6,12,18,24]
;
if (!d.name eq 'X') then wshow
erase
;
week_arr = ['Saturday', 'Sunday', 'Monday']
;
xtickname = ref_day(t.day-day0) + ' ' + strmid(gt_time(t, /str, /nolead), 0, 5)
axis, 0, off_arr(0), xrange = [0,nhours], xtitle = 'Universal Time / England (Nov-Dec, Jan-Apr)', $
	xticks = n-1, xtickv=xtickv, xtickname=xtickname, xstyle=1	;, ymargin=[4,0]
off = int2secarr(t(0), '2-jan-93 8:00')/60./60.
for i=-1,1 do oplot,indgen(11)-off+i*24, intarr(11)+off_arr(0)+1.5, thick=3*(i+1)^1.6
;for i=-1,1 do xyouts, -off+i*24+5.5, off_arr(0)+2, week_arr(i+1), align=0.5
for i=0,1 do xyouts, -off+i*24+5.5, off_arr(0)+2, week_arr(i+1), align=0.5
;
t2 = anytim2ints(t, off=1*60.*60)	;BST is 1 hour ahead of UT
xtickname = ref_day(t2.day-day0) + ' ' + strmid(gt_time(t2, /str, /nolead), 0, 5)
axis, 0, off_arr(1), xrange = [0,nhours], xtitle = 'British Summer Time (May thru Oct)', $
	xticks = n-1, xtickv=xtickv, xtickname=xtickname, xstyle=1	;, ymargin=[4,0]
off = int2secarr(t2(0), '2-jan-93 8:00')/60./60.
for i=-1,1 do oplot,indgen(11)-off+i*24, intarr(11)+off_arr(1)+1.5, thick=3*(i+1)^1.6
;for i=-1,1 do xyouts, -off+i*24+5.5, off_arr(1)+2, week_arr(i+1), align=0.5
for i=0,1 do xyouts, -off+i*24+5.5, off_arr(1)+2, week_arr(i+1), align=0.5
;
t2 = anytim2ints(t, off=9*60.*60)	;JST is 9 hrs ahead of UT
xtickname = ref_day(t2.day-day0) + ' ' + strmid(gt_time(t2, /str, /nolead), 0, 5)
axis, 0, off_arr(2), xrange = [0,nhours], xtitle = 'Japan Standard Time', $
	xticks = n-1, xtickv=xtickv, xtickname=xtickname, xstyle=1	;, ymargin=[4,0]
off = int2secarr(t2(0), '2-jan-93 8:00')/60./60.
for i=-1,1 do oplot,indgen(11)-off+i*24, intarr(11)+off_arr(2)+1.5, thick=3*(i+1)^1.6
for i=-1,1 do xyouts, -off+i*24+5.5, off_arr(2)+2, week_arr(i+1), align=0.5
;
t2 = anytim2ints(t, off=-7*60.*60)	;PST is 7 hrs behind of UT
xtickname = ref_day(t2.day-day0) + ' ' + strmid(gt_time(t2, /str, /nolead), 0, 5)
axis, 0, off_arr(3), xrange = [0,nhours], xtitle = 'Pacific Daylight Savings Time (May thru Oct)', $
	xticks = n-1, xtickv=xtickv, xtickname=xtickname, xstyle=1	;, ymargin=[4,0]
off = int2secarr(t2(0), '2-jan-93 8:00')/60./60.
for i=-1,1 do oplot,indgen(11)-off+i*24, intarr(11)+off_arr(3)+1.5, thick=3*(i+1)^1.6
;for i=-1,1 do xyouts, -off+i*24+5.5, off_arr(3)+2, week_arr(i+1), align=0.5
for i=-1,0 do xyouts, -off+i*24+5.5, off_arr(3)+2, week_arr(i+1), align=0.5
;
t2 = anytim2ints(t, off=-8*60.*60)	;PDT  is 8 hrs behind of UT
xtickname = ref_day(t2.day-day0) + ' ' + strmid(gt_time(t2, /str, /nolead), 0, 5)
axis, 0, off_arr(4), xrange = [0,nhours], xtitle = 'Pacific Standard Time (Nov-Dec, Jan-Apr)', $
	xticks = n-1, xtickv=xtickv, xtickname=xtickname, xstyle=1	;, ymargin=[4,0]
off = int2secarr(t2(0), '2-jan-93 8:00')/60./60.
for i=-1,1 do oplot,indgen(11)-off+i*24, intarr(11)+off_arr(4)+1.5, thick=3*(i+1)^1.6
;for i=-1,1 do xyouts, -off+i*24+5.5, off_arr(4)+2, week_arr(i+1), align=0.5
for i=-1,0 do xyouts, -off+i*24+5.5, off_arr(4)+2, week_arr(i+1), align=0.5
;
;JST is +16 from PST
;
if (not keyword_set(guide)) then begin
    plottime, 0, 0, 'TIME_LINE.PRO
    plottime
end
;
end