pro x_date,date,X,Cl00,crn,Cl01

    if n_params( ) lt 1 then begin
       print,'x_date,date,X,Cl0'
       return
    endif
    year_sec_date,date,year,month,day,hh,mm,ss
    ny=year-1995+1
    yeart=INDGEN(ny)+1995
    X=365
    for yn=0,ny-1 do begin
        if yn lt (ny-1) then begin
           res=yeart(yn) MOD 4
           if res gt 0 then X=X+365 else X=X+366
        endif else begin
           dmy2doy,year,month,day,doy
           X=X+doy
        endelse
    endfor
    Cl00=349.03 - (360.* X / 27.2753)
    crn=-FIX(Cl00/360)+1891
    Cl00=-Cl00-crn*360
    a=1.91787 & b=-0.13067 & c=-8.25278E-2 & d=-0.17505
    e=365.27116 & f=0.26318  & g=-26379.45   & h=5.20448E-3
    i=-5.56336E-3  & j=-1.22842E-2
    pi=!PI
    Fit = f + X/g + a*SIN(2*pi*X/e) + b*SIN(4*pi*X/e) + h*SIN(6*pi*X/e) $
              + c*COS(2*pi*X/e) + d*COS(4*pi*X/e) + i*COS(6*pi*X/e) 
    Cl01=Cl00+Fit
    print,Cl00,-crn,-Cl01

end