#Dome Flat
imdel tmp0,tmp1,tmp2,tmp3,mask,mask2,irs1_cont_BANDm.fits
imarith irs1_cont_BAND.fits / ../flat tmp2
fit1d tmp2 tmp3 fit axis=2 inter- function=legendre order=1 high_reject=2 low_reject=2 niterate=5 grow=0
imrepl tmp3 2000 lower=INDEF upper=2000
median tmp3 tmp3 300 1
imarith tmp2 / tmp3 tmp3
imarith tmp3 * 10000. tmp0.fits
displ tmp0 1 zs+
median tmp0 tmp1 3 3
imarith tmp0 - tmp1 tmp1
imarith tmp1 * ../badpixr tmp1
imrepl tmp1 0. lower=INDEF upper=0.
imrepl tmp1 65000. upper=INDEF lower=65000.
!rm tmp.tbl tmp1.tbl
imhist tmp1 z1=-50 z2=65050 binwidt=100 nbins=651 hist_type=cumulative listout+ > tmp.tbl
!../thresh tmp.tbl 1300 tmp1 > tmp1.tbl
cl < tmp1.tbl
imstati tmp1 lower=INDEF upper=0.5
imarith ../badpixr * tmp1 mask.fits
convolve mask.fits mask.fits "0. 0.2 0.; 0.2 0.2 0.2; 0. 0.2 0."
imrepl mask 0. lower=INDEF upper=0.9
displ mask 2 zs- zr- z1=1 z2=0
imstati mask lower=INDEF upper=0.5
imcopy ../badpix2r mask2.fits
!../fixbadpix tmp0.fits irs1_cont_BANDm.fits
displ irs1_cont_BANDm 3 zs+

imdel tmp1,irs1_cont_BANDe.fits,irs1_cont_BANDs.fits,irs1_cont_BANDa.fits,check,bkg,tmp1_bpm
!echo "0.5 1 1 1 1 1 1 1 0.5" > extract.dat
#imsur irs1_cont_BANDm tmp1.fits 5 1 type_ou=response functio=spline3 
imcopy irs1_cont_BANDm tmp1.fits
!../extract tmp1.fits irs1_cont_BANDe.fits > irs1_cont_BANDe.dat
displ check 4 zs+
displ irs1_cont_BANDe 5 zs+
median bkg bkg 27 9
!../extract2 bkg.fits bkg.fits
magnify bkg bkg 1. 25. y1=1.01 y2=9.0 interpo=linear boundar=nearest fluxcon-
displ bkg 6 zs+
imarith irs1_cont_BANDe.fits - bkg irs1_cont_BANDs.fits
displ irs1_cont_BANDs 7 zs+
!../align irs1_cont_BANDs.fits irs1_cont_BANDa.fits > irs1_cont_BANDa.dat
#!../align irs1_cont_BANDe.fits irs1_cont_BANDa.fits > irs1_cont_BANDa.dat
imdel irs1_cont_BANDa_bak
imrename irs1_cont_BANDa.fits irs1_cont_BANDa_bak.fits
displ irs1_cont_BANDa_bak 8 zs+
!rm tmp.tbl
imstati irs1_cont_BANDa_bak.fits fields=mean format- > tmp.tbl
imarith irs1_cont_BANDa_bak.fits / @tmp.tbl irs1_cont_BANDa.fits
median irs1_cont_BANDa.fits irs1_cont_BANDa.fits 400 1 zloreje=0.1
displ irs1_cont_BANDa 9 zs+
