Practice Reduction on RD300 (14 oct 01) before going on 2002 Nov Run

# Try a rectify-then-proceed-as-usual reduction.

# Before rectifying, overscan-subraction with esibias.cl.
# Note that it seems you gotta load noao.imred.ccdred in
# order for esibias to get at the ccdproc task.  Also, 
# nothing worked until I "unlearned" ccdproc.

task esibias=/kinney/sdawson/keck/02nov/esi_testredux/esibias.cl

ls *fits > list.bias
vi list.bias
:.,$s/.fits//g

esibias @list.bias

# Also before rectifying, you need to apply an x-shift.  I think  (NOT SURE ABOUT THIS)
# Arjun determined this empirically.
# Syntax = imshift esi0023 esi0023s 27 0
# NOTE: this looks like it screwed up the data, but was good for the calibrations.  Weird.

ls *.fits > list.xshift
vi list.xshift
:.,$s/.fits/ 27 0/
:.,$s/esi/imshift esi/g

cl < list.xshift

# Now try rectifying.  Make a shell script.
# Syntax = esishift esi_rect.specs_new esi0023.fits esi0023r.fits

ls *s.fits > list.rect
vi list.rect
:.,$s/esi/esishift esi_rect.specs_new esi/g

source list.rect

# At this point, you might flip it and try running bogus (though
# it might choke on the fact that the spectra are many too many
# pixels long in the dispersion direction.  As an alternative,
# cut out the order of interest and do a longslit-type recution.

imcopy esi0024r[1134:1264,*] esi0024r_o8
imcopy esi0025r[1134:1264,*] esi0025r_o8

background ("esi0024r_o8",
"esi0024r_o8s", axis=1, interactive=no, sample="*", naverage=1,
function="spline3", order=4, low_reject=3., high_reject=2., niterate=3,
grow=0., graphics="stdgraph", cursor="")

background ("esi0025r_o8",
"esi0025r_o8s", axis=1, interactive=no, sample="*", naverage=1,
function="spline3", order=4, low_reject=3., high_reject=2., niterate=3,
grow=0., graphics="stdgraph", cursor="")

display esi0024r_o8s 1
display esi0025r_o8s 2

# Determine shifts for imcombine
# For some reason, there's an x-offset in the splot

splot esi0024r_o8sA
  disp = 1
  nsum = 3000

  x = 153.2

splot esi0025r_o8sA

  x = 136.0

	0	0	
	20	0