pro quickdark, side=side if n_elements(side) eq 0 then side = 512 wdef,0,side infil = findfile('/yd3/reformat/sfr*') tstart = anytim2ints(!stime,offset=-3600.*28*24) rd_roadmap,infil,rm tt = tim2dset(rm,tstart) rm1 = rm ss = where(gt_filta(rm1) eq 1 and gt_filtb(rm1) eq 3 and $ gt_dpe(rm1) eq 13 and $ gt_percentd(rm1) eq 100 and gt_expmode(rm1) eq 1) rd_xda,infil,ss,index,data tots = reform(total(total(sxt_decomp(data),1),1))/(256.^2)/4 utplot,index,tots,psym=-2,ytit='DN/sec/HRpix',yr=[28,36],tit='Recent dark level' fmt_timer,index,t1,t2 outplot,[t1,t2],[31.0894,31.0894],lines=1 saa = where(saa_in(index)) evt_grid,index(saa) ;print,stdev(tots,mean),mean end