PRO ohtimser,ps=ps

r1=ddread('\data\zenith\93123103.txt')
HELP,r1
q1=ddread('\data\zenith\93123150.txt')
HELP,q1
r2=ddread('\data\zenith\94010102.txt')
HELP,r2
r2(0,*)=r2(0,*)+86400.
q2=ddread('\data\zenith\94010152.txt')
HELP,q2
q2(0,*)=q2(0,*)+86400.
r=[transpose(r1),transpose(r2)]
q=[transpose(q1),transpose(q2)]
tr=r(*,0)
ir=r(*,1)
tq=q(*,0)
iq=q(*,1)
ir=ir/0.7356
iq=iq/0.6581
tr=[tr(0:2),tq(3),tr(3:*)]
ir=[ir(0:2),(ir(2)+ir(3))/2.,ir(3:*)]
rqratio=ir/iq
trot=FINDGEN(7)*20.+165.
rqvst=[0.1449/0.1857,0.1552/0.1691,0.1643/0.1542 $
      ,0.1720/0.1433,0.1794/0.1330,0.1859/0.1242 $
      ,0.1927/0.1169]

IF KEYWORD_SET(ps) $
    THEN psopen,'\csa\ohtimser.ps',/long $
    ELSE window,0,xs=600,ys=750

!p.multi=[0,0,2,0,0]
plot,tr/3600.,ir,psym=-1,xrange=[0,42],xticks=7,xtickv=INDGEN(8)*6 $
    ,xtickname=['!1700','06','12','18','00','06','12','18'],xminor=3 $
    ,xtitle='!1731 Dec 1993          Hours (UT)          1 Jan 1994' $
    ,ytitle='!17OH (6-2) R / Q Branch Brightness (R)' $
    ,title='!17Zenith OH Brightness, 31 Dec 1993 - 1 Jan 1994'
oplot,tq/3600.,iq,psym=-2
xyouts,2,470,'!17R branch (829 nm)'
xyouts,2,160,'!17Q branch (835 nm)'

plot,tq/3600.,rqratio,yran=[1.,1.5],xrange=[0,42],xticks=7 $
    ,xtickv=INDGEN(8)*6,xminor=3 $
    ,xtickname=['!1700','06','12','18','00','06','12','18'] $
    ,xtitle='!1731 Dec 1993          Hours (UT)          1 Jan 1994' $
    ,ytitle='!17R / Q Branch Brightness Ratio' $
    ,title='!17Mesopause Rotational Temperature Inferred from I(R)/I(Q)'
for i=2,5 do begin
    plots,[0,42],replicate(rqvst(i),2)
    xyouts,36.,rqvst(i)-0.02,STRING(fix(trot(i)),format='(I3," K")')
    endfor

ggs=TRANSPOSE(ddread('\csa\ohtemp.933'))
tiggs=ggs(*,0)
teggs=ggs(*,1)
OPLOT,tiggs,teggs,line=1

IF KEYWORD_SET(ps) THEN psclose

return
end