raw1=readfits('~/data/m42/5s1.FIT')
raw2=readfits('~/data/m42/5s2.FIT')
raw3=readfits('~/data/m42/5s3.FIT')
raw4=readfits('~/data/m42/5s4.FIT')
raw5=readfits('~/data/m42/5s5.FIT')
raw6=readfits('~/data/m42/5s6.FIT')
raw7=readfits('~/data/m42/5s7.FIT')

dark=readfits('~/telescope_data/masters/scaledmasterdark.FIT')
bias=readfits('~/telescope_data/masters/masterbias.FIT')
mask=readfits('~/telescope_data/masters/sig5mask.FIT')


dark=dark/(360)
a1=raw1-bias-dark
a2=raw2-bias-dark
a3=raw3-bias-dark
a4=raw4-bias-dark
a5=raw5-bias-dark
a6=raw6-bias-dark
a7=raw7-bias-dark

.compile '~/telescope_data/maskfix/maskfixnospace'

a1=maskfix(a1,mask)
a2=maskfix(a2,mask)
a3=maskfix(a3,mask)
a4=maskfix(a4,mask)
a5=maskfix(a5,mask)
a6=maskfix(a6,mask)
a7=maskfix(a7,mask)

sharplim=[0.2,1.0]
roundlim=[-1.0,1.0]

find,a1,x1,y1,flux1,s,r,120,8,roundlim,sharplim

find,a2,x2,y2,flux2,s,r,120,8,roundlim,sharplim

find,a3,x3,y3,flux3,s,r,120,8,roundlim,sharplim

find,a4,x4,y4,flux4,s,r,120,8,roundlim,sharplim

find,a5,x5,y5,flux5,s,r,120,8,roundlim,sharplim

find,a6,x6,y6,flux6,s,r,120,8,roundlim,sharplim

find,a7,x7,y7,flux7,s,r,120,8,roundlim,sharplim


x1=x1(REVERSE(sort(flux1)))
y1=y1(REVERSE(sort(flux1)))
flux1=flux1(REVERSE(sort(flux1)))
x2=x2(REVERSE(sort(flux2)))
y2=y2(REVERSE(sort(flux2)))
flux2=flux2(REVERSE(sort(flux2)))
x3=x3(REVERSE(sort(flux3)))
y3=y3(REVERSE(sort(flux3)))
flux3=flux3(REVERSE(sort(flux3)))
x4=x4(REVERSE(sort(flux4)))
y4=y4(REVERSE(sort(flux4)))
flux4=flux4(REVERSE(sort(flux4)))
x5=x5(REVERSE(sort(flux5)))
y5=y5(REVERSE(sort(flux5)))
flux5=flux5(REVERSE(sort(flux5)))
x6=x6(REVERSE(sort(flux6)))
y6=y6(REVERSE(sort(flux6)))
flux6=flux6(REVERSE(sort(flux6)))
x7=x7(REVERSE(sort(flux7)))
y7=y7(REVERSE(sort(flux7)))
flux7=flux7(REVERSE(sort(flux7)))

info=fltarr(3,7)
info(0,0)=x1(0)
info(1,0)=y1(0)
info(2,0)=flux1(0)
info(0,1)=x2(0)
info(1,1)=y2(0)
info(2,1)=flux2(0)
info(0,2)=x3(0)
info(1,2)=y3(0)
info(2,2)=flux3(0)
info(0,3)=x4(1)
info(1,3)=y4(1)
info(2,3)=flux4(0)
info(0,4)=x5(0)
info(1,4)=y5(0)
info(2,4)=flux5(0)
info(0,5)=x6(0)
info(1,5)=y6(0)
info(2,5)=flux6(0)
info(0,6)=x7(0)
info(1,6)=y7(0)
info(2,6)=flux7(0)

xrange=Round(max(info(0,*)))-round(min(info(0,*)))
yrange=Round(max(info(1,*)))-round(min(info(1,*)))

stack=fltarr(765+xrange,510+yrange,7)

offset=intarr(2,7)
offset(0,0)=Round(info(0,0))-Round(MIN(info(0,*)))
offset(1,0)=Round(info(1,0))-Round(MIN(info(1,*)))
offset(0,1)=Round(info(0,1))-Round(MIN(info(0,*)))
offset(1,1)=Round(info(1,1))-Round(MIN(info(1,*)))
offset(0,2)=Round(info(0,2))-Round(MIN(info(0,*)))
offset(1,2)=Round(info(1,2))-Round(MIN(info(1,*)))
offset(0,3)=Round(info(0,3))-Round(MIN(info(0,*)))
offset(1,3)=Round(info(1,3))-Round(MIN(info(1,*)))
offset(0,4)=Round(info(0,4))-Round(MIN(info(0,*)))
offset(1,4)=Round(info(1,4))-Round(MIN(info(1,*)))
offset(0,5)=Round(info(0,5))-Round(MIN(info(0,*)))
offset(1,5)=Round(info(1,5))-Round(MIN(info(1,*)))
offset(0,6)=Round(info(0,6))-Round(MIN(info(0,*)))
offset(1,6)=Round(info(1,6))-Round(MIN(info(1,*)))

for i=0,764 do begin $
for j=0,509 do begin $
stack(i+offset(0,0),j+offset(1,0),0)=a1(i,j) &$
stack(i+offset(0,1),j+offset(1,1),1)=a2(i,j) &$
stack(i+offset(0,2),j+offset(1,2),2)=a3(i,j) &$
stack(i+offset(0,3),j+offset(1,3),3)=a4(i,j) &$
stack(i+offset(0,4),j+offset(1,4),4)=a5(i,j) &$
stack(i+offset(0,5),j+offset(1,5),5)=a6(i,j) &$
stack(i+offset(0,6),j+offset(1,6),6)=a7(i,j) &$
endfor &$
endfor