pro test14,crn,i1,i2,di,dt121,dt122,dt123,dt124

    if n_params( ) lt 1 then begin
       print,'test14,crn,i1,i2,di,dt121,dt122,dt123,dt124'
       print,'crn=2075 & crl1=150.550 & crl2=190.505'
       return
    endif
  ; crn=2075 
    crnt=INDGEN(15)+crn
    crl1=150.550 & crl2=190.505
    dt121=fltarr(15)
    dt122=fltarr(15)
    dt123=fltarr(15)
    dt124=fltarr(15)
;for ii=0,14 do begin
for ii=i1,i2,di do begin
    crni=crnt(ii)
; 1
print,'crn,crl:',crni,crl1
    date_crncrlz1,crni,crl1,date1
print,'crn,crl:',crni,crl2
    date_crncrlz1,crni,crl2,date2
    fdoy_date,date1,fdoy1,year1
    fdoy_date,date2,fdoy2,year2
    dt=fdoy2-fdoy1
    CASE 1 OF
       dt gt 54: begin
          md=year2 MOD 4
          if md eq 0 then fdoy1=fdoy1+366 else fdoy1=fdoy1+365
                 end
       dt lt -54: begin
          md=year1 MOD 4
          if md eq 0 then fdoy2=fdoy2+366 else fdoy2=fdoy2+365
                  end
       else: print,' '
    ENDCASE
    dt1=fdoy2-fdoy1
; print,'z1-dt1;',dt1
    dt121(ii)=dt1
  ; print,date1,' ',date2
  ; print,'z1-dt1:',dt1
; 2
    CASE 1 OF
     crl1 lt 10: scrl1='00'+STRTRIM(crl1,2)
     crl1 lt 100: scrl1='0'+STRTRIM(crl1,2)
     else: scrl1=STRTRIM(crl1,2)
    ENDCASE
    scrncl0=STRTRIM(crni,2)+':'+scrl1
    time_crncl0,scrncl0,date1
    CASE 1 OF
     crl2 lt 10: scrl2='00'+STRTRIM(crl2,2)
     crl2 lt 100: scrl2='0'+STRTRIM(crl2,2)
     else: scrl2=STRTRIM(crl2,2)
    ENDCASE
    scrncl0=STRTRIM(crni,2)+':'+scrl2
    time_crncl0,scrncl0,date2
    fdoy_date,date1,fdoy1,year1
    fdoy_date,date2,fdoy2,year2
    dt=fdoy2-fdoy1
    CASE 1 OF
       dt gt 54: begin
          md=year2 MOD 4
          if md eq 0 then fdoy1=fdoy1+366 else fdoy1=fdoy1+365
                 end
       dt lt -54: begin
          md=year1 MOD 4
          if md eq 0 then fdoy2=fdoy2+366 else fdoy2=fdoy2+365
                  end
       else: print,' '
    ENDCASE
    dt2=fdoy2-fdoy1
;print,dt2
    dt122(ii)=dt2
  ; print,date1,' ',date2
  ; print,'z2-dt2:',dt2
; 3 
if ii eq 0 then begin    
      dct=(1-crl2/360.) - (1-crl1/360.)
      dt4=dct*27.2753
      for k=0,14 do dt124(k)=dt4
endif 
 
 print,'z3-crn,crl:',crni,crl1
    date_crncrlz3,crni,crl1,date1
 print,'z3-crn,crl:',crni,crl2
    date_crncrlz3,crni,crl2,date2
    fdoy_date,date1,fdoy1,year1
    fdoy_date,date2,fdoy2,year2
    dt=fdoy2-fdoy1
    CASE 1 OF
       dt gt 54: begin
          md=year2 MOD 4
          if md eq 0 then fdoy1=fdoy1+366 else fdoy1=fdoy1+365
                 end
       dt lt -54: begin
          md=year1 MOD 4
          if md eq 0 then fdoy2=fdoy2+366 else fdoy2=fdoy2+365
                  end
       else: print, ' '
    ENDCASE
  ; print,date1,' ',date2
    dt3=fdoy2-fdoy1
    dt123(ii)=dt3
  ; print,'z3-dt3:',dt3
endfor
plt,dt121,dt122,dt123,dt124

end