pro date0_crncrl,crn,crl,date0,PLT=plt

    if n_params( ) lt 1 then begin
       print,'date0_crncrl,crn,crl,date0,/plt'
       return
    endif
    degs=crl-(crn-1891)*360.           ;Total degrees from 19950101_000000
    fday0=(349.03-degs)*27.2753/360    ;Total fractional days from 19950101
   ;dds=FIX(fday0)
    ny0=FIX(fday0/365)                 ;Number of years if all are 365 days
    if ny0 eq 0 then begin
       year1=1995
       fdoy=fday0+1
    endif else begin
       yeart=INDGEN(ny0)+1995
       fday=-1
       for yn=0,ny0-1 do begin
           yeari=yeart(yn)
           res=yeari MOD 4
           if res gt 0 then fday=fday+365 else fday=fday+366 
           dfday=fday-fday0
           if fday gt fday0 then goto,lbl
           year1=yeari+1
           fdoy=fday0-fday
       endfor
    endelse
lbl: year=year1
     doy=FIX(fdoy)
     doy2dmy,year,doy,day,month
     dfday=fdoy-doy
     fhh=dfday*24
     hh=FIX(fhh)
     fmm=(fhh-hh)*60
     mm=FIX(fmm)
     ss=FIX((fmm-mm)*60)
     if hh eq 23 then begin
        day=day+1
        hh=0
     endif
     date_year_sec,year,month,day,hh,mm,ss,date0     
     if keyword_set(plt) then print,'date0:',date0

end