# Modified version of Dave Cook's recipe on his website
#
# Frames are binned 2x2 resulting in 2048x2048 images
# Can't really use flats with exposures 2 sec or less -- too much closing-shutter effect

# 1. Update headers to indicate what type of image each fits file represents.

hedit @bias.list ccdtype bias add+ verify-
hedit @bias.list imagetype zero add+ verify-
hedit @object.list ccdtype object add+ verify-
hedit @object.list imagetype object add+ verify-
hedit @flat.list ccdtype flat add+ verify-
hedit @flat.list imagetype flat add+ verify-

# 2. Combine bias images using zerocombine.  

task prevu=/d/crow1/tools/iraf/galphot/sep00/prevu.cl
zerocombine input=@bias.list output=masterbias.fit combine=average reject=minmax ccdtype=zero process=no delete=no clobber=no

# 3. Determine the overscan and trim sections.

displ a021 1
implot a021
:l 1000 1100
:c 1000 1100
[ typing 'e' twice will define a zoomed-in box and 'r' replots; trimsec=[1:2048,1:2048] -- the entire image -- and I didn't find any apparent overscan region ]

# 4. Trim the flats and objects. 

ccdproc @flat.list ccdtype=flat fixpix=no overscan=no trim=yes zerocor=yes darkcor=no flatcor=no trimsec=[1:2048,1:2048] zero=masterbias
ccdproc @object.list ccdtype=object fixpix=no overscan=no trim=yes zerocor=yes darkcor=no flatcor=no trimsec=[1:2048,1:2048] zero=masterbias

# 5. Combine flats. MAKE SURE TO COMBINE SEPARATE FILTERS SEPARATELY. 

flatcombine input=@flat.list output=master_flat.fit combine=median reject=avsigclip ccdtype=flat subsets=yes process=no delete=no clobber=no scale=mode
flatcombine input=@flat.subset.list output=master_flat.subset.fit combine=median reject=avsigclip ccdtype=flat subsets=yes process=no delete=no clobber=no scale=mode

# 6. Flat-field your object frames for flats

ccdproc @object.list ccdtype=object fixpix=no overscan=no trim=no zerocor=no darkcor=no flatcor=yes flat=master_flat.subset
      
# 7. Align and stack images; adapted from /crow1/wysh/stack/combine.doc
      
######################################################
Do the photometry on the images.  Use prevu.cl to determine the best center for all images in both filters
phot a023 output=a023.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a024 output=a024.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a025 output=a025.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a026 output=a026.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a027 output=a027.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a028 output=a028.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a029 output=a029.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a030 output=a030.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a031 output=a031.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a032 output=a032.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
phot a033 output=a033.phot radplots=yes centerpars.calgori=centroid centerpars.cbox=10 fitskypars.annulus=15 fitskypars.dannulus=5 photpars.apertur=15
########################################
# create .txt version of .phot by putting in star coords in columns 3&4 (reference star coords in columns 1&2
ls a*txt > ToCoordinates.list
geomap @ToCoordinates.list Coordinates.list 1 2048 1 2048               # hit 'q' in Textronic for each frame
cp object.subset.list Aligned.list				
n Aligned.list&								# add '.aligned' to each name
geotran @object.subset.list @Aligned.list Coordinates.list @ToCoordinates.list
prevu @Aligned.list
imcombine @Aligned.list stack.fit combine=average reject=avsigclip      # weight=@weighting.bin
imstat @Aligned.list
imstat stack[1610:1830,630:960]
fitsky stack fitskypars.annulus=0 fitskypars.dannulus=10 
######## add wcs to header #####################################################
gasp
n a*.txt&								# pop up coordinates from 'master' and 'slave' .mag.1 files
displ ngc4214_dss[0] 2 					                
displ stack 1
imexa rimexam.radius=10 logfile=dss.txt keep+						# get coordinates with 'a' for the same 5 stars
xyeq iminfo=yes image=ngc4214_dss[0] xy=dss.txt pix_center=iraf xcol=1 ycol=2 nskip=1 orig=yes new=no ra_h=yes ra_f=%h dec_f=%h > dss.coords
# copy and paste into ccmap.in the stack pix coords using the first two columns from any *.txt file; cut and paste ra+dec from dss.coords into 3rd & 4th columns
n dss.coords&
n ccmap.in& 
ccmap input=        ccmap.in database=ccmap.database images=stack.fit lngunits=hours latunits=degrees results=ccmap.out update=yes interactive=no 

############ bootstrap calibration based on existing calibrated imaging ##################################
############ LVL BVRI at; SDSS at pawnee:/pawnee3/davec/LVLbroadband/Data/sdss/v4_sdss_images_unzip ################################

i.e., my neighbor has a 6cm-long mouse, and mine is half that length.  How long is my mouse?