# Try iqpkg try: import iqpkg import iqutils except: print "No iqpkg, iqutils. Can't run lmi_stats." # Necessary packages iraf.images() iraf.immatch() # iraf.imfilter() iraf.noao() iraf.imred() iraf.ccdred() iraf.stsdas() iraf.hst_calib() iraf.nicmos() iraf.imutil() yes = iraf.yes no = iraf.no INDEF = iraf.INDEF globclob = yes globver = yes LMIFILTS = ["U", "B", "V", "R", "I", "SDSS-G", "SDSS-R", "SDSS-I", "SDSS-Z", "SDSS-U"] LMIPIXSCALE = 0.240 ASTROMCMD = "/Users/scenko/python/astrometry/vlt_autoastrometry.py" FDICT = {"SDSS-U": "UMAG", "SDSS-G": "GMAG", "SDSS-R": "RMAG", "SDSS-I": "IMAG", "SDSS-Z": "ZMAG"} SEEDICT = {"SDSS-U": "co", "SDSS-G": "go", "SDSS-R": "ro", "SDSS-I": "bo", "SDSS-Z": "ko"} LMIGAIN = 2.95
# These statements import python modules that will be used below... import glob import os import string from pyraf import iraf from shutil import copyfile # Load necessary IRAF (pyraf) modules... iraf.stsdas() iraf.hst_calib() iraf.synph # cd to current working directory iraf.cd(os.getcwd()) # This creates a list of filter file names... print "Here is the list of CCDs:" filterFileNameList = glob.glob('des_band_20150601a_syn.average.list') filterFileNameList += glob.glob('des_band_20150601a_syn.ccd*.list') filterFileNameList += glob.glob('des_band_20150601a_syn.ccd*_amp1.list') filterFileNameList += glob.glob('des_band_20150601a_syn.ccd*_amp2.list') filterFileNameList.sort() print filterFileNameList # This is the file containing the list of calspec files, # which you created previously using the pyraf "files" # command (see top of this script)... calspecListFileName = 'pickles.list' # Open and the calspec list file... fd = open(calspecListFileName, 'r')