msu_xrt_cat,'28-oct-2021 15:00','28-oct-2021 18:00',cat,ofiles
;wdef,0,1000,600
;plot,cat.xcen,cat.ycen,psym=1
;
ss_fl=where(cat.flflg eq 'FLR' and cat.chip_sum ne 8) 
files0=ofiles(ss_fl)
mreadfits,files0,index0,/nodata
ss_384=where(index0.naxis1 eq 384)
mreadfits,files0(ss_384),index,data
stepper,index,safe_log10(data)
;
;***** PREP *****
index=index(0:29) & data=data(*,*,0:29)
xrt_prep,index,data,iout,dout,/float,/norm
;
NEXT:
;***** coalignment *****
jj_ref=10
dout2=dout
print,iout(jj_ref).xcen,iout(jj_ref).ycen
sftx=(iout.xcen-iout(jj_ref).xcen)/iout.platescl
sfty=(iout.ycen-iout(jj_ref).ycen)/iout.platescl
for i=0,n_elements(iout)-1 do $
   dout2(*,*,i)=shift(dout(*,*,i),round(sftx(i)),round(sfty(i)))
;
stepper,iout,safe_log10(dout2)
;
end