a1=readfits('~/telescope_data/Observing/24.03.11/hip31635/60s1.FIT')
a2=readfits('~/telescope_data/Observing/24.03.11/hip31635/60s2.FIT')
a3=readfits('~/telescope_data/Observing/24.03.11/hip31635/60s3.FIT')
a4=readfits('~/telescope_data/Observing/24.03.11/hip31635/60s4.FIT')
a5=readfits('~/telescope_data/Observing/24.03.11/hip31635/60s5.FIT')

flat=readfits('~/telescope_data/Observing/07.03.11/Flats/flat.FIT')
dark=readfits('~/telescope_data/masters/scaledmasterdark.FIT')
dark=dark/30
bias=readfits('~/telescope_data/masters/masterbias.FIT')


ab1=(a1-bias-dark)/flat
ab2=(a2-bias-dark)/flat
ab3=(a3-bias-dark)/flat
ab4=(a4-bias-dark)/flat
ab5=(a5-bias-dark)/flat

;--- files were saved at this point as FITS and then read back in ---
;
;ab1=readfits('~/telescope_data/Observing/24.03.11/hip31635/darkbiasflat1.FIT')
;ab2=readfits('~/telescope_data/Observing/24.03.11/hip31635/darkbiasflat2.FIT')
;ab3=readfits('~/telescope_data/Observing/24.03.11/hip31635/darkbiasflat3.FIT')
;ab4=readfits('~/telescope_data/Observing/24.03.11/hip31635/darkbiasflat4.FIT')
;ab5=readfits('~/telescope_data/Observing/24.03.11/hip31635/darkbiasflat5.FIT')
;
;--- But without this step the program would continue as follows ----

maskhot=readfits('~/telescope_data/masters/sig5mask.FIT')
maskwarm=readfits('~/telescope_data/masters/chmask.FIT')
maskcool=readfits('~/telescope_data/masters/chcmask.FIT')

.compile '~/telescope_data/maskfix/maskfixnospace.pro'

ma1=maskfix(ab1,maskhot) &$
ma2=maskfix(ab2,maskhot) &$
ma3=maskfix(ab3,maskhot) &$
ma4=maskfix(ab4,maskhot) &$
ma5=maskfix(ab5,maskhot)

mb1=maskfix(ab1,maskwarm) &$
mb2=maskfix(ab2,maskwarm) &$
mb3=maskfix(ab3,maskwarm) &$
mb4=maskfix(ab4,maskwarm) &$
mb5=maskfix(ab5,maskwarm)

mc1=maskfix(ab1,maskcool) &$
mc2=maskfix(ab2,maskcool) &$
mc3=maskfix(ab3,maskcool) &$
mc4=maskfix(ab4,maskcool) &$
mc5=maskfix(ab5,maskcool)

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

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

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

find,ma3,x3,y3,flux3,s,r,120,8,[-1.5,1.0],[0.1,1.0]
; This file was difficult. The star was slightly streaked ans so the sharpness
;   and round limits had to be adjusted in order to locate the correct star.

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

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


; The located stars are sorted by flux intensity. The bightest star is the r'

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)))

; The x,y coords and flux of brightest star from each 'find' routine is stored.

info=fltarr(3,5)
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(0)
info(1,3)=y4(0)
info(2,3)=flux4(0)
info(0,4)=x5(0)
info(1,4)=y5(0)
info(2,4)=flux5(0)

; Errors - IGNORE THIS SECTION. MODIFIED VERSION IMPLEMENTED LATER.
;
; error for hot pixel masking
; Set bad pixels to -200 so that they're well out of the data range and
; can be distinguished later with APER.
;erra1=float(ab1)
;erra1(where(maskhot eq 0))=-1
;erra2=float(ab2)
;erra2(where(maskhot eq 0))=-1
;erra3=float(ab3)
;erra3(where(maskhot eq 0))=-1
;erra4=float(ab4)
;erra4(where(maskhot eq 0))=-1
;erra5=float(ab5)
;erra5(where(maskhot eq 0))=-1
;
;errcropa1=fltarr(160,160)
;errcropa2=fltarr(160,160)
;errcropa3=fltarr(160,160)
;errcropa4=fltarr(160,160)
;errcropa5=fltarr(160,160)
;for i=0,159 do begin $
;for j=0,159 do begin $
;errcropa1(i,j)=erra1(ROUND(info(0,0))-80+i,ROUND(info(1,0))-80+j) & $
;errcropa2(i,j)=erra2(ROUND(info(0,1))-80+i,ROUND(info(1,1))-80+j) & $
;errcropa3(i,j)=erra3(ROUND(info(0,2))-80+i,ROUND(info(1,2))-80+j) & $
;errcropa4(i,j)=erra4(ROUND(info(0,3))-80+i,ROUND(info(1,3))-80+j) & $
;errcropa5(i,j)=erra5(ROUND(info(0,4))-80+i,ROUND(info(1,4))-80+j) & $
;endfor & $
;endfor
;
;skyrad=[70,89]
;apr=[20,30,40,50,60,70]
;phpadu=2.3
;badpix=[-200,2000]
;
;APER,errcropa1,80,80,errmaga1,errerra1,errskya1,errskyerra1,phpadu,apr,skyrad,badpix
;APER,errcropa2,80,80,errmaga2,errerra2,errskya2,errskyerra2,phpadu,apr,skyrad,badpix
;APER,errcropa3,80,80,errmaga3,errerra3,errskya3,errskyerra3,phpadu,apr,skyrad,badpix
;APER,errcropa4,80,80,errmaga4,errerra4,errskya4,errskyerra4,phpadu,apr,skyrad,badpix
;APER,errcropa5,80,80,errmaga5,errerra5,errskya5,errskyerra5,phpadu,apr,skyrad,badpix
;
; Magnitude for cropped and fixed data
;
;cropa1=fltarr(160,160)
;cropa2=fltarr(160,160)
;cropa3=fltarr(160,160)
;cropa4=fltarr(160,160)
;cropa5=fltarr(160,160)
;for i=0,159 do begin $
;for j=0,159 do begin $
;cropa1(i,j)=ma1(ROUND(info(0,0))-80+i,ROUND(info(1,0))-80+j) & $
;cropa2(i,j)=ma2(ROUND(info(0,1))-80+i,ROUND(info(1,1))-80+j) & $
;cropa3(i,j)=ma3(ROUND(info(0,2))-80+i,ROUND(info(1,2))-80+j) & $
;cropa4(i,j)=ma4(ROUND(info(0,3))-80+i,ROUND(info(1,3))-80+j) & $
;cropa5(i,j)=ma5(ROUND(info(0,4))-80+i,ROUND(info(1,4))-80+j) & $
;endfor & $
;endfor
;
;skyrad=[60,70]
;apr=[10,20,30,40,50,60]
;phpadu=2.3
;
;APER,cropa1,80,80,maga1,erra1,skya1,skyerra1,phpadu,apr,skyrad
;
;APER,cropa2,80,80,maga2,erra2,skya2,skyerra2,phpadu,apr,skyrad
;
;APER,cropa3,80,80,maga3,erra3,skya3,skyerra3,phpadu,apr,skyrad
;
;APER,cropa4,80,80,maga4,erra4,skya4,skyerra4,phpadu,apr,skyrad
;
;APER,cropa5,80,80,maga5,erra5,skya5,skyerra5,phpadu,apr,skyrad
;






; Magnitude fixed data HOT MASK

; Divide by exposure length?!?
ma1=ma1/60
ma2=ma2/60
ma3=ma3/60
ma4=ma4/60
ma5=ma5/60

skyrad=[100,120]
apr=[2,5,8,10,15,20,40,60,80,100]

APER,ma1,info(0,0),info(1,0),maga1,erra1,skya1,skyerra1,phpadu,apr,skyrad,/FLUX

APER,ma2,info(0,1),info(1,1),maga2,erra2,skya2,skyerra2,phpadu,apr,skyrad,/FLUX

APER,ma3,info(0,2),info(1,2),maga3,erra3,skya3,skyerra3,phpadu,apr,skyrad,/FLUX

APER,ma4,info(0,3),info(1,3),maga4,erra4,skya4,skyerra4,phpadu,apr,skyrad,/FLUX

APER,ma5,info(0,4),info(1,4),maga5,erra5,skya5,skyerra5,phpadu,apr,skyrad,/FLUX







; Magnitude fixed data WARM MASK

skyrad=[100,120]
apr=[2,5,8,10,15,20,40,60,80,100]

mb1=mb1/60
mb2=mb2/60
mb3=mb3/60
mb4=mb4/60
mb5=mb5/60

APER,mb1,info(0,0),info(1,0),magb1,errb1,skyb1,skyerrb1,phpadu,apr,skyrad,/FLUX

APER,mb2,info(0,1),info(1,1),magb2,errb2,skyb2,skyerrb2,phpadu,apr,skyrad,/FLUX

APER,mb3,info(0,2),info(1,2),magb3,errb3,skyb3,skyerrb3,phpadu,apr,skyrad,/FLUX

APER,mb4,info(0,3),info(1,3),magb4,errb4,skyb4,skyerrb4,phpadu,apr,skyrad,/FLUX

APER,mb5,info(0,4),info(1,4),magb5,errb5,skyb5,skyerrb5,phpadu,apr,skyrad,/FLUX






; Magnitude fixed data COLD MASK

skyrad=[100,120]
apr=[2,5,8,10,15,20,40,60,80,100]

mc1=mc1/60
mc2=mc2/60
mc3=mc3/60
mc4=mc4/60
mc5=mc5/60

APER,mc1,info(0,0),info(1,0),magc1,errc1,skyc1,skyerrc1,phpadu,apr,skyrad,/FLUX

APER,mc2,info(0,1),info(1,1),magc2,errc2,skyc2,skyerrc2,phpadu,apr,skyrad,/FLUX

APER,mc3,info(0,2),info(1,2),magc3,errc3,skyc3,skyerrc3,phpadu,apr,skyrad,/FLUX

APER,mc4,info(0,3),info(1,3),magc4,errc4,skyc4,skyerrc4,phpadu,apr,skyrad,/FLUX

APER,mc5,info(0,4),info(1,4),magc5,errc5,skyc5,skyerrc5,phpadu,apr,skyrad,/FLUX





; Magnitudes for RAW data

skyrad=[100,120]
apr=[2,5,8,10,15,20,40,60,80,100]

reduced1=ab1/60
reduced2=ab2/60
reduced3=ab3/60
reduced4=ab4/60
reduced5=ab5/60

APER,reduced1,info(0,0),info(1,0),magr1,errr1,skyr1,skyerrr1,phpadu,apr,skyrad,/FLUX

APER,reduced2,info(0,1),info(1,1),magr2,errr2,skyr2,skyerrr2,phpadu,apr,skyrad,/FLUX

APER,reduced3,info(0,2),info(1,2),magr3,errr3,skyr3,skyerrr3,phpadu,apr,skyrad,/FLUX

APER,reduced4,info(0,3),info(1,3),magr4,errr4,skyr4,skyerrr4,phpadu,apr,skyrad,/FLUX

APER,reduced5,info(0,4),info(1,4),magr5,errr5,skyr5,skyerrr5,phpadu,apr,skyrad,/FLUX




; Multiply the image by the mask and perform Aper on that
; Does this give a better error estimation?

hotmasked1=ab1*maskhot/60
hotmasked2=ab2*maskhot/60
hotmasked3=ab3*maskhot/60
hotmasked4=ab4*maskhot/60
hotmasked5=ab5*maskhot/60

warmmasked1=ab1*maskwarm/60
warmmasked2=ab2*maskwarm/60
warmmasked3=ab3*maskwarm/60
warmmasked4=ab4*maskwarm/60
warmmasked5=ab5*maskwarm/60

coolmasked1=ab1*maskcool/60
coolmasked2=ab2*maskcool/60
coolmasked3=ab3*maskcool/60
coolmasked4=ab4*maskcool/60
coolmasked5=ab5*maskcool/60


;--- Dist_circle method ---

dist_circle,circle1,[765,510],info(0,0),info(1,0)
dist_circle,circle2,[765,510],info(0,1),info(1,1)
dist_circle,circle3,[765,510],info(0,2),info(1,2)
dist_circle,circle4,[765,510],info(0,3),info(1,3)
dist_circle,circle5,[765,510],info(0,4),info(1,4)

; Hot masking

skyhotmasked1=Total(hotmasked1(where(circle1 gt 100 AND circle1 lt 120)))/Total(maskhot(where(circle1 gt 100 AND circle1 lt 120)))
skyhotmasked2=Total(hotmasked2(where(circle2 gt 100 AND circle2 lt 120)))/Total(maskhot(where(circle2 gt 100 AND circle2 lt 120)))
skyhotmasked3=Total(hotmasked3(where(circle3 gt 100 AND circle3 lt 120)))/Total(maskhot(where(circle3 gt 100 AND circle3 lt 120)))
skyhotmasked4=Total(hotmasked4(where(circle4 gt 100 AND circle4 lt 120)))/Total(maskhot(where(circle4 gt 100 AND circle4 lt 120)))
skyhotmasked5=Total(hotmasked5(where(circle5 gt 100 AND circle5 lt 120)))/Total(maskhot(where(circle5 gt 100 AND circle5 lt 120)))

hotmasked1=hotmasked1-skyhotmasked1
hotmasked2=hotmasked2-skyhotmasked2
hotmasked3=hotmasked3-skyhotmasked3
hotmasked4=hotmasked4-skyhotmasked4
hotmasked5=hotmasked5-skyhotmasked5
;hotmasked1(where(hotmasked1 lt 0))=0
;hotmasked2(where(hotmasked2 lt 0))=0
;hotmasked3(where(hotmasked3 lt 0))=0
;hotmasked4(where(hotmasked4 lt 0))=0
;hotmasked5(where(hotmasked5 lt 0))=0

;Just turn the bad pixels back to zero. Now a good sky pixel will average about zero (hopefully)
hotmasked1=hotmasked1*maskhot
hotmasked2=hotmasked2*maskhot
hotmasked3=hotmasked3*maskhot
hotmasked4=hotmasked4*maskhot
hotmasked5=hotmasked5*maskhot


print,total(hotmasked1(where(circle1 le 2))),$
total(hotmasked1(where(circle1 le 5))),$
total(hotmasked1(where(circle1 le 8))),$
total(hotmasked1(where(circle1 le 10))),$
total(hotmasked1(where(circle1 le 15))),$
total(hotmasked1(where(circle1 le 20))),$
total(hotmasked1(where(circle1 le 40))),$
total(hotmasked1(where(circle1 le 60))),$
total(hotmasked1(where(circle1 le 80))),$
total(hotmasked1(where(circle1 le 100))) &$
$
print,total(hotmasked2(where(circle2 le 2))),$
total(hotmasked2(where(circle2 le 5))),$
total(hotmasked2(where(circle2 le 8))),$
total(hotmasked2(where(circle2 le 10))),$
total(hotmasked2(where(circle2 le 15))),$
total(hotmasked2(where(circle2 le 20))),$
total(hotmasked2(where(circle2 le 40))),$
total(hotmasked2(where(circle2 le 60))),$
total(hotmasked2(where(circle2 le 80))),$
total(hotmasked2(where(circle2 le 100))) &$
$
print,total(hotmasked3(where(circle3 le 2))),$
total(hotmasked3(where(circle3 le 5))),$
total(hotmasked3(where(circle3 le 8))),$
total(hotmasked3(where(circle3 le 10))),$
total(hotmasked3(where(circle3 le 15))),$
total(hotmasked3(where(circle3 le 20))),$
total(hotmasked3(where(circle3 le 40))),$
total(hotmasked3(where(circle3 le 60))),$
total(hotmasked3(where(circle3 le 80))),$
total(hotmasked3(where(circle3 le 100))) &$
$
print,total(hotmasked4(where(circle4 le 2))),$
total(hotmasked4(where(circle4 le 5))),$
total(hotmasked4(where(circle4 le 8))),$
total(hotmasked4(where(circle4 le 10))),$
total(hotmasked4(where(circle4 le 15))),$
total(hotmasked4(where(circle4 le 20))),$
total(hotmasked4(where(circle4 le 40))),$
total(hotmasked4(where(circle4 le 60))),$
total(hotmasked4(where(circle4 le 80))),$
total(hotmasked4(where(circle4 le 100))) &$
$
print,total(hotmasked5(where(circle5 le 2))),$
total(hotmasked5(where(circle5 le 5))),$
total(hotmasked5(where(circle5 le 8))),$
total(hotmasked5(where(circle5 le 10))),$
total(hotmasked5(where(circle5 le 15))),$
total(hotmasked5(where(circle5 le 20))),$
total(hotmasked5(where(circle5 le 40))),$
total(hotmasked5(where(circle5 le 60))),$
total(hotmasked5(where(circle5 le 80))),$
total(hotmasked5(where(circle5 le 100)))


; Warm masking

skywarmmasked1=Total(warmmasked1(where(circle1 gt 100 AND circle1 lt 120)))/Total(maskwarm(where(circle1 gt 100 AND circle1 lt 120)))
skywarmmasked2=Total(warmmasked2(where(circle2 gt 100 AND circle2 lt 120)))/Total(maskwarm(where(circle2 gt 100 AND circle2 lt 120)))
skywarmmasked3=Total(warmmasked3(where(circle3 gt 100 AND circle3 lt 120)))/Total(maskwarm(where(circle3 gt 100 AND circle3 lt 120)))
skywarmmasked4=Total(warmmasked4(where(circle4 gt 100 AND circle4 lt 120)))/Total(maskwarm(where(circle4 gt 100 AND circle4 lt 120)))
skywarmmasked5=Total(warmmasked5(where(circle5 gt 100 AND circle5 lt 120)))/Total(maskwarm(where(circle5 gt 100 AND circle5 lt 120)))

warmmasked1=warmmasked1-skywarmmasked1
warmmasked2=warmmasked2-skywarmmasked2
warmmasked3=warmmasked3-skywarmmasked3
warmmasked4=warmmasked4-skywarmmasked4
warmmasked5=warmmasked5-skywarmmasked5
;warmmasked1(where(warmmasked1 lt 0))=0
;warmmasked2(where(warmmasked2 lt 0))=0
;warmmasked3(where(warmmasked3 lt 0))=0
;warmmasked4(where(warmmasked4 lt 0))=0
;warmmasked5(where(warmmasked5 lt 0))=0
warmmasked1=warmmasked1*maskwarm
warmmasked2=warmmasked2*maskwarm
warmmasked3=warmmasked3*maskwarm
warmmasked4=warmmasked4*maskwarm
warmmasked5=warmmasked5*maskwarm

print,total(warmmasked1(where(circle1 le 2))),$
total(warmmasked1(where(circle1 le 5))),$
total(warmmasked1(where(circle1 le 8))),$
total(warmmasked1(where(circle1 le 10))),$
total(warmmasked1(where(circle1 le 15))),$
total(warmmasked1(where(circle1 le 20))),$
total(warmmasked1(where(circle1 le 40))),$
total(warmmasked1(where(circle1 le 60))),$
total(warmmasked1(where(circle1 le 80))),$
total(warmmasked1(where(circle1 le 100))) &$
$
print,total(warmmasked2(where(circle2 le 2))),$
total(warmmasked2(where(circle2 le 5))),$
total(warmmasked2(where(circle2 le 8))),$
total(warmmasked2(where(circle2 le 10))),$
total(warmmasked2(where(circle2 le 15))),$
total(warmmasked2(where(circle2 le 20))),$
total(warmmasked2(where(circle2 le 40))),$
total(warmmasked2(where(circle2 le 60))),$
total(warmmasked2(where(circle2 le 80))),$
total(warmmasked2(where(circle2 le 100))) &$
$
print,total(warmmasked3(where(circle3 le 2))),$
total(warmmasked3(where(circle3 le 5))),$
total(warmmasked3(where(circle3 le 8))),$
total(warmmasked3(where(circle3 le 10))),$
total(warmmasked3(where(circle3 le 15))),$
total(warmmasked3(where(circle3 le 20))),$
total(warmmasked3(where(circle3 le 40))),$
total(warmmasked3(where(circle3 le 60))),$
total(warmmasked3(where(circle3 le 80))),$
total(warmmasked3(where(circle3 le 100))) &$
$
print,total(warmmasked4(where(circle4 le 2))),$
total(warmmasked4(where(circle4 le 5))),$
total(warmmasked4(where(circle4 le 8))),$
total(warmmasked4(where(circle4 le 10))),$
total(warmmasked4(where(circle4 le 15))),$
total(warmmasked4(where(circle4 le 20))),$
total(warmmasked4(where(circle4 le 40))),$
total(warmmasked4(where(circle4 le 60))),$
total(warmmasked4(where(circle4 le 80))),$
total(warmmasked4(where(circle4 le 100))) &$
$
print,total(warmmasked5(where(circle5 le 2))),$
total(warmmasked5(where(circle5 le 5))),$
total(warmmasked5(where(circle5 le 8))),$
total(warmmasked5(where(circle5 le 10))),$
total(warmmasked5(where(circle5 le 15))),$
total(warmmasked5(where(circle5 le 20))),$
total(warmmasked5(where(circle5 le 40))),$
total(warmmasked5(where(circle5 le 60))),$
total(warmmasked5(where(circle5 le 80))),$
total(warmmasked5(where(circle5 le 100)))


; Cool masking

skycoolmasked1=Total(coolmasked1(where(circle1 gt 100 AND circle1 lt 120)))/Total(maskcool(where(circle1 gt 100 AND circle1 lt 120)))
skycoolmasked2=Total(coolmasked2(where(circle2 gt 100 AND circle2 lt 120)))/Total(maskcool(where(circle2 gt 100 AND circle2 lt 120)))
skycoolmasked3=Total(coolmasked3(where(circle3 gt 100 AND circle3 lt 120)))/Total(maskcool(where(circle3 gt 100 AND circle3 lt 120)))
skycoolmasked4=Total(coolmasked4(where(circle4 gt 100 AND circle4 lt 120)))/Total(maskcool(where(circle4 gt 100 AND circle4 lt 120)))
skycoolmasked5=Total(coolmasked5(where(circle5 gt 100 AND circle5 lt 120)))/Total(maskcool(where(circle5 gt 100 AND circle5 lt 120)))

coolmasked1=coolmasked1-skycoolmasked1
coolmasked2=coolmasked2-skycoolmasked2
coolmasked3=coolmasked3-skycoolmasked3
coolmasked4=coolmasked4-skycoolmasked4
coolmasked5=coolmasked5-skycoolmasked5
;coolmasked1(where(coolmasked1 lt 0))=0
;coolmasked2(where(coolmasked2 lt 0))=0
;coolmasked3(where(coolmasked3 lt 0))=0
;coolmasked4(where(coolmasked4 lt 0))=0
;coolmasked5(where(coolmasked5 lt 0))=0
coolmasked1=coolmasked1*maskcool
coolmasked2=coolmasked2*maskcool
coolmasked3=coolmasked3*maskcool
coolmasked4=coolmasked4*maskcool
coolmasked5=coolmasked5*maskcool

print,total(coolmasked1(where(circle1 le 2))),$
total(coolmasked1(where(circle1 le 5))),$
total(coolmasked1(where(circle1 le 8))),$
total(coolmasked1(where(circle1 le 10))),$
total(coolmasked1(where(circle1 le 15))),$
total(coolmasked1(where(circle1 le 20))),$
total(coolmasked1(where(circle1 le 40))),$
total(coolmasked1(where(circle1 le 60))),$
total(coolmasked1(where(circle1 le 80))),$
total(coolmasked1(where(circle1 le 100))) &$
$
print,total(coolmasked2(where(circle2 le 2))),$
total(coolmasked2(where(circle2 le 5))),$
total(coolmasked2(where(circle2 le 8))),$
total(coolmasked2(where(circle2 le 10))),$
total(coolmasked2(where(circle2 le 15))),$
total(coolmasked2(where(circle2 le 20))),$
total(coolmasked2(where(circle2 le 40))),$
total(coolmasked2(where(circle2 le 60))),$
total(coolmasked2(where(circle2 le 80))),$
total(coolmasked2(where(circle2 le 100))) &$
$
print,total(coolmasked3(where(circle3 le 2))),$
total(coolmasked3(where(circle3 le 5))),$
total(coolmasked3(where(circle3 le 8))),$
total(coolmasked3(where(circle3 le 10))),$
total(coolmasked3(where(circle3 le 15))),$
total(coolmasked3(where(circle3 le 20))),$
total(coolmasked3(where(circle3 le 40))),$
total(coolmasked3(where(circle3 le 60))),$
total(coolmasked3(where(circle3 le 80))),$
total(coolmasked3(where(circle3 le 100))) &$
$
print,total(coolmasked4(where(circle4 le 2))),$
total(coolmasked4(where(circle4 le 5))),$
total(coolmasked4(where(circle4 le 8))),$
total(coolmasked4(where(circle4 le 10))),$
total(coolmasked4(where(circle4 le 15))),$
total(coolmasked4(where(circle4 le 20))),$
total(coolmasked4(where(circle4 le 40))),$
total(coolmasked4(where(circle4 le 60))),$
total(coolmasked4(where(circle4 le 80))),$
total(coolmasked4(where(circle4 le 100))) &$
$
print,total(coolmasked5(where(circle5 le 2))),$
total(coolmasked5(where(circle5 le 5))),$
total(coolmasked5(where(circle5 le 8))),$
total(coolmasked5(where(circle5 le 10))),$
total(coolmasked5(where(circle5 le 15))),$
total(coolmasked5(where(circle5 le 20))),$
total(coolmasked5(where(circle5 le 40))),$
total(coolmasked5(where(circle5 le 60))),$
total(coolmasked5(where(circle5 le 80))),$
total(coolmasked5(where(circle5 le 100)))

; -------
;
;
;skyrad=[100,120]
;apr=[2,5,8,10,15,20,40,60,80,100]
;
;APER,hotmasked1,info(0,0),info(1,0),magea1,errea1,skyea1,skyerrea1,phpadu,apr,skyrad,/FLUX
;
;APER,hotmasked2,info(0,1),info(1,1),magea2,errea2,skyea2,skyerrea2,phpadu,apr,skyrad,/FLUX
;
;APER,hotmasked3,info(0,2),info(1,2),magea3,errea3,skyea3,skyerrea3,phpadu,apr,skyrad,/FLUX
;
;APER,hotmasked4,info(0,3),info(1,3),magea4,errea4,skyea4,skyerrea4,phpadu,apr,skyrad,/FLUX
;
;APER,hotmasked5,info(0,4),info(1,4),magea5,errea5,skyea5,skyerrea5,phpadu,apr,skyrad,/FLUX
;
;
;APER,warmmasked1,info(0,0),info(1,0),mageb1,erreb1,skyeb1,skyerreb1,phpadu,apr,skyrad,/FLUX
;
;APER,warmmasked2,info(0,1),info(1,1),mageb2,erreb2,skyeb2,skyerreb2,phpadu,apr,skyrad,/FLUX
;
;APER,warmmasked3,info(0,2),info(1,2),mageb3,erreb3,skyeb3,skyerreb3,phpadu,apr,skyrad,/FLUX
;
;APER,warmmasked4,info(0,3),info(1,3),mageb4,erreb4,skyeb4,skyerreb4,phpadu,apr,skyrad,/FLUX
;
;APER,warmmasked5,info(0,4),info(1,4),mageb5,erreb5,skyeb5,skyerreb5,phpadu,apr,skyrad,/FLUX
;
;
;APER,coolmasked1,info(0,0),info(1,0),magec1,errec1,skyec1,skyerrec1,phpadu,apr,skyrad,/FLUX
;
;APER,coolmasked2,info(0,1),info(1,1),magec2,errec2,skyec2,skyerrec2,phpadu,apr,skyrad,/FLUX
;
;APER,coolmasked3,info(0,2),info(1,2),magec3,errec3,skyec3,skyerrec3,phpadu,apr,skyrad,/FLUX
;
;APER,coolmasked4,info(0,3),info(1,3),magec4,errec4,skyec4,skyerrec4,phpadu,apr,skyrad,/FLUX
;
;APER,coolmasked5,info(0,4),info(1,4),magec5,errec5,skyec5,skyerrec5,phpadu,apr,skyrad,/FLUX
;



cropa1=fltarr(200,200)
cropa2=fltarr(200,200)
cropa3=fltarr(200,200)
cropa4=fltarr(200,200)
cropa5=fltarr(200,200)
for i=0,199 do begin $
for j=0,199 do begin $
cropa1(i,j)=hotmasked1(ROUND(info(0,0))-100+i,ROUND(info(1,0))-100+j) & $
cropa2(i,j)=hotmasked2(ROUND(info(0,1))-100+i,ROUND(info(1,1))-100+j) & $
cropa3(i,j)=hotmasked3(ROUND(info(0,2))-100+i,ROUND(info(1,2))-100+j) & $
cropa4(i,j)=hotmasked4(ROUND(info(0,3))-100+i,ROUND(info(1,3))-100+j) & $
cropa5(i,j)=hotmasked5(ROUND(info(0,4))-100+i,ROUND(info(1,4))-100+j) & $
endfor & $
endfor