sci_frames = all_frames
outroot = '1131_Ks'
ccdbfile = 'ccmap_final.txt'
skybpm = 'bpm_from_sky.fits'
"""
Do the nonlinearity correction and create the sky file as the first pass
"""
if redpass=='calib_1':
niri.calib_1(sci_frames,rawroot,outroot,bpmfile,rawdir)
""" Calibrate the science frames """
if redpass=='calib_2':
print ""
skyfile = '%s_sky.fits' % outroot
niri.reduce_sci(sci_frames,outroot,skyfile,flatfile)
""" Make the bad pixel mask from the calibrated science frames """
if redpass=='calib_3':
if skybpm is not None:
niri.niri_bpm_from_sky(sci_frames,skybpm,outsky='sky_from_ff.fits')
else:
print ""
print "*** No bad pixel mask requested for this set of input files ***"
print ""
"""
Prepare the files for running SExtractor
The input files will be called ff*fits
The output files will be called fp*fits
"""
fullnames.append('%s_%d' % (t4obsdate,i))
else:
for i in sci_frames:
fullnames.append('%s_%d' % (obsdate,i))
"""
Do the nonlinearity correction and Create the sky file as the first pass
"""
if redpass=='make_sky':
niri.calib_1(sci_frames,rawroot,outroot,bpmfile,rawdir,obsdate)
""" Do the final calibration and make the SExtractor catalogs """
if redpass=='final_calib':
print ""
skyfile = '%s_sky.fits' % outroot
niri.reduce_sci(fullnames,outroot,skyfile,flatfile)
if skybpm is not None:
niri.niri_bpm_from_sky(fullnames,skybpm)
""" Make the SExtractor catalogs """
if redpass=='make_cats':
ff_files = []
for i in fullnames:
ff_files.append('ff%s.fits' % i)
niri.niri_sextractor(ff_files,tilebpm,fixpix=True,catformat='ascii')
"""
Redo the astrometry of all the individual input files onto the full grid
now that we have created it, and then do the final coadd.
"""