Robert and I have written code that does what is written in the description of this ticket, i.e., identifies n-sigma outliers from master darks and master flats. This code has been used with AuxTel, TS8, and comCam data. The code is in the following file: FindDefects.py
Documentation can be found in the doctrings inside the code, but the the function that does the main calculations is "find_hot_and_cold_pixels". The code creates a directory ("output_directory") where ASCII files with the defects (as boxes) are printed. In addition, histograms and other diagnose plots are written to PDF files in that same directory.