pro annulv22010,crl21,crl22,dt122x1,dt1240r1,dt122x2,dt1240r2

    scrl21=strtrim(crl21,2)
    scrl22=strtrim(crl22,2)
  ; ttl1='crl1=180, crl2='+scrl21
  ; ttl2='crl1=180, crl2='+scrl22
    ttl1='crl1=200, crl2=180
    ttl2='crl1=200, crl2=25
    ofl='annualvariation2010.ps'
    dt122x1=dt122x1*24
    dt122x2=dt122x2*24
;   dt1220=dt1220*24
;   dt123=dt123*24
;   dt124=dt124*24
    dt1240r1=dt1240r1*24
    dt1240r2=dt1240r2*24
    
    tek_color
    psout,021,ofl
    !p.multi=[0,1,2]
    mx=MAX(dt122x1) & mn=MIN(dt122x1)
    plot,dt122x1,xrange=[0,14],yrange=[mn,mx],xsty=1,ysty=1,title=ttl1,xtitle='CR Number from 1928 to 1941 (1997:10:10--1998:10:08)',ytitle='Time difference from crl1 to crl2 in hrs',/nodata
    dy=(mx-mn)/18.
    yy=indgen(18)*dy+mn
    oplot,dt1240r1,line=0,color=4,thick=3
    oplot,[7.5,8.5],[yy(15),yy(15)],line=0,color=4,thick=3
    xyouts,9.0,yy(15),'Negleting eccentricity'
    oplot,dt122x1,color=2,thick=3
    oplot,[7.5,8.5],[yy(13),yy(13)],color=2,thick=3
    xyouts,9.0,yy(13),'Including eccentricity'
    ddt=(dt122x1-dt1240r1)
    mx=STRMID(STRTRIM(MAX(ddt),2),0,5) & mn=STRMID(STRTRIM(MIN(ddt),2),0,6)
    xyouts,1.0,yy(2),'Difference between red and blue lines:'
    xyouts,1.0,yy(1),'Maximum='+mx+' hrs, Minimum='+mn+' hrs'
;
    mx=MAX(dt122x2) & mn=MIN(dt122x2)
    plot,dt122x2,xrange=[0,14],yrange=[mn,mx],xsty=1,ysty=1,title=ttl2,xtitle='CR Number from 1928 to 1941 (1997:10:10--1998:10:08)',ytitle='Time difference from crl1 to crl2 in hrs',/nodata
    dy=(mx-mn)/18.
    yy=indgen(18)*dy+mn
    oplot,dt1240r2,line=0,color=4,thick=3
    oplot,[7.5,8.5],[yy(15),yy(15)],line=0,color=4,thick=3
    xyouts,9.0,yy(15),'Negleting eccentricity'
    oplot,dt122x2,color=2,thick=3
    oplot,[7.5,8.5],[yy(13),yy(13)],color=2,thick=3
    xyouts,9.0,yy(13),'Including eccentricity'
    ddt=(dt122x2-dt1240r2)
    mx=STRMID(STRTRIM(MAX(ddt),2),0,5) & mn=STRMID(STRTRIM(MIN(ddt),2),0,6)
    xyouts,1.0,yy(2),'Difference between red and blue lines:'
    xyouts,1.0,yy(1),'Maximum='+mx+' hrs, Minimum='+mn+' hrs'
    !p.multi=0
    psdone
    print,ofl

end