; K bandHD 1881 NSFCam Data - February 24, 2002 UT
; Will make postscript file.  It will NOT output to xwindow.

; enter the name of the postscript file
psfile = 'hd1881_K_flat.ps'

; Cha
nge plotting scale if needed
min = -100 & max = 10000

; Read in the images?
readin = 1

; Flatten then images?
flatten = 1

; Set the platescale
platescale = 0.3

; NOTE: Data is read in with E to the left and N at the bottom

if readin eq 1 then k104a = readfits('/scrs1/cew/feb24/data0104.a')
if readin eq 1 then k105a = readfits('/scrs1/cew/feb24/data0105.a')
if readin eq 1 then k106a = readfits('/scrs1/cew/feb24/data0106.a')
if readin eq 1 then k107a = readfits('/scrs1/cew/feb24/data0107.a')
if readin eq 1 then k108a = readfits('/scrs1/cew/feb24/data0108.a')

flats,'hd1881_k.list',hd1881_k_sky

; Subtract object and sky
k0o = k104a-hd1881_k_sky
k1o = k105a-hd1881_k_sky
k2o = k106a-hd1881_k_sky
k3o = k107a-hd1881_k_sky
k4o = k108a-hd1881_k_sky



; Flip the images in N-S so N is on top
k0 = reverse(k0o,2)
k1 = reverse(k1o,2)
k2 = reverse(k2o,2)
k3 = reverse(k3o,2)
k4 = reverse(k4o,2)


; Flatten if selected
if flatten eq 1 then restore,'/scrs1/deh/feb23/flats.idl'

if flatten eq 1 then k0 = k0/ksub_flat
if flatten eq 1 then k1 = k1/ksub_flat
if flatten eq 1 then k2 = k2/ksub_flat
if flatten eq 1 then k3 = k3/ksub_flat
if flatten eq 1 then k4 = k4/ksub_flat


; Try and make a badpixel mask
;kmask = fltarr(256,256)

;k0bad = where(k0 gt 1000)
;k1bad = where(k1 gt 1000)
;k2bad = where(k2 gt 1000)
;k3bad = where(k3 gt 1000)
;k4bad = where(k4 gt 1000)
;k5bad = where(k5 gt 1000)

;kmask(k0bad) = 1.0
;kmask(k1bad) = 1.0
;kmask(k2bad) = 1.0
;kmask(k3bad) = 1.0
;kmask(k4bad) = 1.0
;kmask(k5bad) = 1.0

; Set some plot parameters - This looks best for postscript
!x.thick=5
!y.thick=5
!p.thick=4
!p.charsize=1.5
!x.title = '!17R. A. Offset (arcsec)'
!y.title = '!17Dec. Offset (arcsec)'
!x.style=1
!y.style=1
!x.ticklen=0.04
!y.ticklen=0.04

; The center location of the comet
x0= 133. & y0 = 127. & x1 = 108. & y1 = 154. & x2 = 108. & y2 = 100. 
x3 = 163. & y3= 99. & x4 = 165. & y4 = 152. 

; Set pixel scale of array in arcsecond
nband = size(k0)
xa = platescale*findgen(256)
ya = platescale*findgen(256)

; Set up PS parameters
set_plot,'ps'
device,/portrait,/inches,filename=psfile,xoffset=1,yoffset=2,xsize=6.5,ysize=8

; Multiple images on one sheet
!p.multi = [0,2,3]

; START PLOTTING

x = x0*platescale - xa & y = -y0*platescale + ya
dispcont3,bytscl(k0,min=min,max=max),x,y,10000,aspect=1
xyouts,x(0)-2,y(0)+54*platescale,'C/2002 C1 IK 02/24/2002 UT',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+40*platescale,'NSFCam at NASA IRTF',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+27*platescale,'K,0.3"/pix, itime=3s*5,AM=2.415',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+13.5*platescale,'90" N-Nod, center = ',color=255,charsize=0.75
xyouts,x(0)-135*platescale,y(0)+13.5*platescale,fix(x1),color=255,charsize=0.75
xyouts,x(0)-176*platescale,y(0)+13.5*platescale,fix(y1),color=255,charsize=0.75
xyouts,x(0)-(x(0)-x(nband(1)-1))/2.,y(nband(2)-1)-4.,'N',color=255,charsize=0.75
xyouts,x(0)-2.,y(0)+(y(nband(2)-1)-y(0))/2.,'E',color=255,charsize=0.75

x = x1*platescale - xa & y = -y1*platescale + ya
dispcont3,bytscl(k1,min=min,max=max),x,y,10000,aspect=1
xyouts,x(0)-2,y(0)+54*platescale,'C/2002 C1 IK 02/24/2002 UT',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+40*platescale,'NSFCam at NASA IRTF',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+27*platescale,'K,0.3"/pix, itime=3s*5,AM=2.415',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+13.5*platescale,'90" N-Nod, center = ',color=255,charsize=0.75
xyouts,x(0)-135*platescale,y(0)+13.5*platescale,fix(x1),color=255,charsize=0.75
xyouts,x(0)-176*platescale,y(0)+13.5*platescale,fix(y1),color=255,charsize=0.75
xyouts,x(0)-(x(0)-x(nband(1)-1))/2.,y(nband(2)-1)-4.,'N',color=255,charsize=0.75
xyouts,x(0)-2.,y(0)+(y(nband(2)-1)-y(0))/2.,'E',color=255,charsize=0.75

x = x2*platescale - xa & y = -y2*platescale + ya
dispcont3,bytscl(k2,min=min,max=max),x,y,10000,aspect=1
xyouts,x(0)-2,y(0)+54*platescale,'C/2002 C1 IK 02/24/2002 UT',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+40*platescale,'NSFCam at NASA IRTF',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+27*platescale,'K,0.3"/pix, itime=3s*5,AM=2.415',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+13.5*platescale,'90" N-Nod, center = ',color=255,charsize=0.75
xyouts,x(0)-135*platescale,y(0)+13.5*platescale,fix(x1),color=255,charsize=0.75
xyouts,x(0)-176*platescale,y(0)+13.5*platescale,fix(y1),color=255,charsize=0.75
xyouts,x(0)-(x(0)-x(nband(1)-1))/2.,y(nband(2)-1)-4.,'N',color=255,charsize=0.75
xyouts,x(0)-2.,y(0)+(y(nband(2)-1)-y(0))/2.,'E',color=255,charsize=0.75

x = x3*platescale - xa & y = -y3*platescale + ya
dispcont3,bytscl(k3,min=min,max=max),x,y,10000,aspect=1
xyouts,x(0)-2,y(0)+54*platescale,'C/2002 C1 IK 02/24/2002 UT',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+40*platescale,'NSFCam at NASA IRTF',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+27*platescale,'H,0.3"/pix, itime=3s*5,AM=2.415',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+13.5*platescale,'90" N-Nod, center = ',color=255,charsize=0.75
xyouts,x(0)-135*platescale,y(0)+13.5*platescale,fix(x1),color=255,charsize=0.75
xyouts,x(0)-176*platescale,y(0)+13.5*platescale,fix(y1),color=255,charsize=0.75
xyouts,x(0)-(x(0)-x(nband(1)-1))/2.,y(nband(2)-1)-4.,'N',color=255,charsize=0.75
xyouts,x(0)-2.,y(0)+(y(nband(2)-1)-y(0))/2.,'E',color=255,charsize=0.75

x = x4*platescale - xa & y = -y4*platescale + ya
dispcont3,bytscl(k4,min=min,max=max),x,y,10000,aspect=1
xyouts,x(0)-2,y(0)+54*platescale,'C/2002 C1 IK 02/24/2002 UT',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+40*platescale,'NSFCam at NASA IRTF',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+27*platescale,'K, 0.3"/pix, itime=3s*5,AM=2.415',color=255,charsize=0.75
xyouts,x(0)-2,y(0)+13.5*platescale,'90" N-Nod, center = ',color=255,charsize=0.75
xyouts,x(0)-135*platescale,y(0)+13.5*platescale,fix(x1),color=255,charsize=0.75
xyouts,x(0)-176*platescale,y(0)+13.5*platescale,fix(y1),color=255,charsize=0.75
xyouts,x(0)-(x(0)-x(nband(1)-1))/2.,y(nband(2)-1)-4.,'N',color=255,charsize=0.75
xyouts,x(0)-2.,y(0)+(y(nband(2)-1)-y(0))/2.,'E',color=255,charsize=0.75

device,/close

set_plot,'x'

!p.multi=0

; To find center

window,0,xsize=500,ysize=500
disp,k0,2,min=min,max=1000,/c
disp,k1,2,min=min,max=1000,/c
disp,k2,2,min=min,max=1000,/c
disp,k3,2,min=min,max=1000,/c
disp,k4,2,min=min,max=1000,/c

end