pro test0619,i1,i2,di

    crnt=INDGEN(15)+2065
    y2008t=INDGEN(13) & for i=0,12 do y2008t(i)=2008
    yeart=[2007,y2008t,2009]
    montht=[12,1,2,3,4,5,6,7,8,8,9,10,11,12,1]
    fdayt=[29.05,25.39,21.73,20.06,16.34,13.58,9.79,6.99,3.20,30.43,26.70,23.99,20.29,17.61,13.94]
    dayt=FIX(fdayt)
    fdoyt=fltarr(15)
    fdoytz=fltarr(15)
    crr=27.2753
    for i=i1,i2,di do begin
        year=yeart(i)
        month=montht(i)
        day=dayt(i)
        fhh=fdayt(i)-day
        dmy2doy,year,month,day,doy
        fdoyt(i)=doy+fhh
    ;   crni=crnt(i) & crl=360.00
    ;   date_crncrlz1,crni,crl,date
    ;   fdoy_date,date,fdoyz,year
    ;   fdoytz(i)=fdoyz
    endfor
    print,'fdoyt:',fdoyt
    print,'fdoytz:',fdoytz
    fdoynt=fdoyt
    dfdoyt=fltarr(15)
    fdoyntz=fdoytz
    dfdoytz=fltarr(15)
    for i=1,14 do begin
        year0=yeart(i-1)
        year1=yeart(i)
        fdoy0=fdoyt(i-1)
        fdoy1=fdoyt(i)
        dfdoy=fdoy1-fdoy0
; print,'year0,fdoy0,fdoy1,dfdoy:',year0,fdoy0,fdoy1,dfdoy
        if dfdoy lt -30 then begin
           md=year0 MOD 4
           if md eq 0 then dfdoyt(i)=fdoy1-(fdoy0-366) $
           else dfdoyt(i)=fdoy1-(fdoy0-365)
        endif else dfdoyt(i)=dfdoy
; Calculated result
        fdoy0z=fdoytz(i-1)
        fdoy1z=fdoytz(i)
        dfdoyz=fdoy1z-fdoy0z
    ;   if dfdoyz lt -30 then begin
    ;      md=year0 MOD 4
    ;      if md eq 0 then dfdoytz(i)=fdoy1z-(fdoy0z-366) $
    ;      else dfdoytz(i)=fdoy1z-(fdoy0z-365)
    ;   endif else dfdoytz(i)=dfdoyz
    endfor
 ;  print,'dfdoyt:',dfdoyt
    dfdoyt(0)=!values.f_nan
    print,'dfdoyt:',dfdoyt
  ; dfdoytz(0)=!values.f_nan
  ; print,'dfdoytz:',dfdoytz
    mx=MAX(dfdoyt) & mn=MIN(dfdoyt)
    tek_color
    ofl='test0619.ps'
    psout,011,ofl
    xtt='Carrington Rotation Numbers from 2007:12:25 to 2009:01:13'
    ytt='Synodic rotation period in days'
    avg=TOTAL(dfdoyt(1:14))/14.
    print,'Avg:',avg
    plot,crnt,dfdoyt,xrange=[2065,2080],yrange=[27.00,27.50],xsty=1,ysty=1,$
         xtitle=xtt,ytitle=ytt,title='Data from Astronomical Almanac 2008'
  ; oplot,crnt,dfdoytz,color=2
    oplot,[2065,2080],[27.2753,27.2753],color=2
    oplot,[2065,2080],[27.2779,27.2779]
    xyouts,2065,27.47,'08:01:25.39'
    xyouts,2070.5,27.47,'2008:07:06.99'
    xyouts,2078,27.47,'09:01:13.94'
    xyouts,2069.5,27.254,'Carrington rotation rate=27.2753'
    xyouts,2069.5,27.276,'Averaged synodic rate='+STRTRIM(avg,2)
  ; xyouts,2069.5,27.224,'Averaged synodic rate='+STRTRIM(avg,2)
    psdone
    print,ofl

end