pro heatercomp ; Created March 2000 by Chuck See ; This program compares the data from the in flight CO2 heater test datasets. ; The test logs are hardcoded below. ; This program plots the ratio of the temperature rise profiles, F3/F1, F5/F1 & F5/F3. ; This program uses the conversion equation for DISR 3 (Deg K = 0.09934*p(1:9,*)+23.973) ; Fortunately the constants are the same for the strap and the auxboard thermistors. ; The heater test samples every 15 seconds for 2 minutes. ; The data is in the pixel portion of the file in two rows, p=intarr(10,2) ; P(0,0) is 1, P(0,1) is 2, P(1:9,0) is the strap data, P(1:9,1) is the aux brd data. cd,'',current=cdir & cd,cdir ;establishes cdir as current directory F1='/df3/23Oct97/Log/L+8day_checkout_23Oct97.3' F3='/df3/28Dec98/Log/F3_28Dec98.1' F5='/df3/03Feb00/Log/F5_03Feb00.1' time=[0,.25,.5,.75,1.,1.25,1.5,1.75,2.] ;the times when the thermistor is read (minutes). ;print,'Log = ',log ;print,'Today is: ',systime(0) ;__________________________________________________________________ ;F1... cd,F1+'/DB/HeatTest' F1files=findfile('*', count = nfiles) ; p is an intarr(10,2) containing the strap and aux board data. ;for i=0,nfiles-1 do begin d_read,F1files,h,p ;Date=strmid(d_value(h,15),4,6)+' '+strmid(d_value(h,15),20,4) Date=strmid(d_value(h,15),4,21) print,'F1 - Date: ',date,string(10b) print,'Raw Data:' print,p print,'' f1data=fltarr(10,2) f1data(0,0)=1 & f1data(0,1)=2 f1data(1:9,*)=0.09934*p(1:9,*)+23.973 ;conversion from DN to deg K print,'Converted Data:' print,f1data, format='(i8,9f8.2)' print,'' ;rate1 is the CCD rate... fpheat=4.91733 ;the standard strap heater rate from F1 f1rate1= 0.09934*(p(9,0)-p(1,0))/2 ; Scale factor = 0.09934 deg/dn print,f1rate1,format='("Focal Plane Heater Rate:",T30,f7.3," deg/min")' print,fpheat,format='("Typical FP Heater Rate:",T30,f7.3," deg/min")' print,100*(f1rate1-fpheat)/fpheat,format='("Difference:",T30,f7.3," %",/)' ;rate2 is the aux board rate... auxheat=1.88746 ;the standard aux. board heater rate from F1 f1rate2= 0.09934*(p(9,1)-p(1,1))/2 ; Scale factor = 0.09934 deg/dn print,f1rate2,format='("Aux Board Heater Rate:",T30,f7.3," deg/min")' print,auxheat,format='("Typical Aux Board Rate:",T30,f7.3," deg/min")' print,100*(f1rate2-auxheat)/auxheat,format='("Difference:",T30,f7.3," %",/)' ;__________________________________________________________________ ;F3... cd,f3+'/DB/HeatTest' f3files=findfile('*', count = nfiles) ; p is an intarr(10,2) containing the strap and aux board data. ;for i=0,nfiles-1 do begin d_read,f3files,h,p ;Date=strmid(d_value(h,15),4,6)+' '+strmid(d_value(h,15),20,4) Date=strmid(d_value(h,15),4,21) print,'f3 - Date: ',date,string(10b) print,'Raw Data:' print,p print,'' f3data=fltarr(10,2) f3data(0,0)=1 & f3data(0,1)=2 f3data(1:9,*)=0.09934*p(1:9,*)+23.973 ;conversion from DN to deg K print,'Converted Data:' print,f3data, format='(i8,9f8.2)' print,'' ;rate1 is the CCD rate... fpheat=4.91733 ;the standard strap heater rate from F1 f3rate1= 0.09934*(p(9,0)-p(1,0))/2 ; Scale factor = 0.09934 deg/dn print,f3rate1,format='("Focal Plane Heater Rate:",T30,f7.3," deg/min")' print,fpheat,format='("Typical FP Heater Rate:",T30,f7.3," deg/min")' print,100*(f3rate1-fpheat)/fpheat,format='("Difference:",T30,f7.3," %",/)' ;rate2 is the aux board rate... auxheat=1.88746 ;the standard aux. board heater rate from F1 f3rate2= 0.09934*(p(9,1)-p(1,1))/2 ; Scale factor = 0.09934 deg/dn print,f3rate2,format='("Aux Board Heater Rate:",T30,f7.3," deg/min")' print,auxheat,format='("Typical Aux Board Rate:",T30,f7.3," deg/min")' print,100*(f3rate2-auxheat)/auxheat,format='("Difference:",T30,f7.3," %",/)' ;__________________________________________________________________ ;f5... cd,f5+'/DB/HeatTest' f5files=findfile('*', count = nfiles) ; p is an intarr(10,2) containing the strap and aux board data. ;for i=0,nfiles-1 do begin d_read,f5files,h,p ;Date=strmid(d_value(h,15),4,6)+' '+strmid(d_value(h,15),20,4) Date=strmid(d_value(h,15),4,21) print,'f5 - Date: ',date,string(10b) print,'Raw Data:' print,p print,'' f5data=fltarr(10,2) f5data(0,0)=1 & f5data(0,1)=2 f5data(1:9,*)=0.09934*p(1:9,*)+23.973 ;conversion from DN to deg K print,'Converted Data:' print,f5data, format='(i8,9f8.2)' print,'' ;rate1 is the CCD rate... fpheat=4.91733 ;the standard strap heater rate from F1 f5rate1= 0.09934*(p(9,0)-p(1,0))/2 ; Scale factor = 0.09934 deg/dn print,f5rate1,format='("Focal Plane Heater Rate:",T30,f7.3," deg/min")' print,fpheat,format='("Typical FP Heater Rate:",T30,f7.3," deg/min")' print,100*(f5rate1-fpheat)/fpheat,format='("Difference:",T30,f7.3," %",/)' ;rate2 is the aux board rate... auxheat=1.88746 ;the standard aux. board heater rate from F1 f5rate2= 0.09934*(p(9,1)-p(1,1))/2 ; Scale factor = 0.09934 deg/dn print,f5rate2,format='("Aux Board Heater Rate:",T30,f7.3," deg/min")' print,auxheat,format='("Typical Aux Board Rate:",T30,f7.3," deg/min")' print,100*(f5rate2-auxheat)/auxheat,format='("Difference:",T30,f7.3," %",/)' ;______________________________________________________________________________________ ;Plot Data... !P.MULTI=[0,1,2,0,0] ;______________________________________________________________________________________ ;Strap (Focal plane) heater... ;F3/F1 Strap_F3overF1=(f3data(2:9,0)-f3data(1,0))/(f1data(2:9,0)-f1data(1,0)) plot,time(1:8),Strap_F3overF1,title='Strap Heater Ratios', $ xtitle='Minutes', ytitle='Normalized Temp. Rise',xrange=[0,2.5],yrange=[0.90,1.10],linestyle=0 ;legend ypos= 0.935 ;y position of legend oplot,[1.8,2.05],[ypos,ypos],linestyle=0 xyouts,2.1,ypos,'F3/F1' ;legend ;F5/F1 Strap_F5overF1=(f5data(2:9,0)-f5data(1,0))/(f1data(2:9,0)-f1data(1,0)) oplot,time(1:8),Strap_F5overF1,linestyle=1 ;legend ypos= 0.925 ;y position of legend xyouts,2.1,ypos,'F5/F1' ;legend oplot,[1.8,2.05],[ypos,ypos],linestyle=1 ;F5/f3 Strap_F5overF3=(f5data(2:9,0)-f5data(1,0))/(f3data(2:9,0)-f3data(1,0)) oplot,time(1:8),Strap_F5overF3,linestyle=2 ;legend ypos= 0.915 ;y position of legend xyouts,2.1,ypos,'F5/F3' ;legend oplot,[1.8,2.05],[ypos,ypos],linestyle=2 ;______________________________________________________________________________________ ;Aux. Board heater... ;F3/F1 Aux_F3overF1=(f3data(2:9,1)-f3data(1,1))/(f1data(2:9,1)-f1data(1,1)) plot,time(1:8),Aux_F3overF1,title='SH Auxiliary Board Heater Ratios', $ xtitle='Minutes', ytitle='Normalized Temp. Rise',xrange=[0,2.5],yrange=[0.90,1.10],linestyle=0 ;legend ypos= 0.935 ;y position of legend oplot,[1.8,2.05],[ypos,ypos],linestyle=0 xyouts,2.1,ypos,'F3/F1' ;legend ;F5/F1 Aux_F5overF1=(f5data(2:9,1)-f5data(1,1))/(f1data(2:9,1)-f1data(1,1)) oplot,time(1:8),Aux_F5overF1,linestyle=1 ;legend ypos= 0.925 ;y position of legend xyouts,2.1,ypos,'F5/F1' ;legend oplot,[1.8,2.05],[ypos,ypos],linestyle=1 ;F5/f3 Aux_F5overF3=(f5data(2:9,1)-f5data(1,1))/(f3data(2:9,1)-f3data(1,1)) oplot,time(1:8),Aux_F5overF3,linestyle=2 ;legend ypos= 0.915 ;y position of legend xyouts,2.1,ypos,'F5/F3' ;legend oplot,[1.8,2.05],[ypos,ypos],linestyle=2 ;______________________________________________________________________________________ ;Create an output file... close,1 Openw,1,'~/Heatercomp.out' printf,1,'This is the output file from heatercomp.pro' printf,1,'It was created on: ', systime(0) printf,1,'The logs examined are: ',string(10b),F1,string(10b),F3,string(10b),F5,string(10b) Printf,1,'Focal Plane Heater Data...' printf,1,' Time F3/F1 F5/F1 F5/F3' for i=0,7 do printf,1,time(i+1),Strap_F3overF1(i),Strap_F5overF1(i),Strap_F5overF3(i), $ format='(f5.2,3F12.4)' Printf,1,string(10b),'Aux Board Heater Data...' printf,1,' Time F3/F1 F5/F1 F5/F3' for i=0,7 do printf,1,time(i+1),Aux_F3overF1(i),Aux_F5overF1(i),Aux_F5overF3(i), $ format='(f5.2,3F12.4)' close,1 cd,cdir !P.MULTI=[0,1,1,0,0] ;stop end