# 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')
