pro crtsnp_crncrl,crn,crl,snp366,crt,snp crt=crn+(1-crl/360.) CASE 1 OF crl lt 10: scrl='00'+STRTRIM(crl,2) crl lt 100: scrl='0'+STRTRIM(crl,2) else: scrl1=STRTRIM(crl1,2) ENDCASE scrncl0=STRTRIM(crn,2)+':'+scrl time_crncl0,scrncl0,date year_sec_date,date,year,month,day,hh,mm,ss dmy2doy,year,month,day,doy snp=snp366(doy) end