pro test14_0624,crn,i1,i2,di,dt121,dt122,dt123,dt124,dt1240,CS3=cs3

    if n_params( ) lt 1 then begin
       print,'test14_0624,crn,i1,i2,di,dt121,dt122,dt123,dt124,dt1240'
       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)
    dt1240=fltarr(15)
    snp=[27.3341,27.3185,27.2931,27.2658,27.2416,27.2307,27.2203,27.2275,27.2452,27.2702,27.2973,27.4060,27.3353]
;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
; 4 
    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
    year_sec_date,date1,year1,month1,day1,hh,mm,ss
    CASE 1 OF
      month1 eq 1 AND day1 gt 15: snr1=27.3400
      month1 eq 2: snr1=27.3400
      month1 eq 3: snr1=27.3300
      month1 eq 4: snr1=27.2800
      month1 eq 5: snr1=27.2400
      month1 eq 6: snr1=27.2100
      month1 eq 7: snr1=27.2000
      month1 eq 8 AND day1 gt 15: snr1=27.2100
      month1 eq 8 AND day1 le 15: snr1=27.2300
      month1 eq 9: snr1=27.2700
      month1 eq 10: snr1=27.2900
      month1 eq 11: snr1=27.3000
      month1 eq 12: snr1=27.3200
      month1 eq 1 AND day1 le 15: snr1=27.3300
    ENDCASE
print,'year1,month1,day1:',year1,month1,day1
    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
    year_sec_date,date2,year2,month2,day2,hh,mm,ss
    CASE 1 OF
      month2 eq 1 AND day2 gt 15: snr2=27.3400
      month2 eq 2: snr2=27.3400
      month2 eq 3: snr2=27.3300
      month2 eq 4: snr2=27.2800
      month2 eq 5: snr2=27.2400
      month2 eq 6: snr2=27.2100
      month2 eq 7: snr2=27.2000
      month2 eq 8 AND day2 gt 15: snr2=27.2100
      month2 eq 8 AND day2 le 15: snr2=27.2300
      month2 eq 9: snr2=27.2700
      month2 eq 10: snr2=27.2900
      month2 eq 11: snr2=27.3000
      month2 eq 12: snr2=27.3200
      month2 eq 1 AND day2 le 15: snr2=27.3300
    ENDCASE
print,'year2,month2,day2:',year2,month2,day2
      if month1 ne month2 then snr=0.5*(snr1+snr2) else snr=snr1
      dct=(1-crl2/360.)-(1-crl1/360.)
    ; dct4=(1-crl2/360.)*snr2 - (1-crl1/360.)*snr1
      dt124(ii)=dct*snr
      dt1240(ii)=dct*27.2753
if keyword_set(cs3) then begin
; 3 
 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
endif
endfor
;dt1240(0)=!values.f_nan
plt0624,dt121,dt122,dt123,dt124,dt1240

end