#basic image reductions
infilelist = sorted(glob.glob('../raw/P*fits'))
imred(infilelist,
'./',
datadir + 'bpm_rss_11.fits',
crthresh=False,
cleanup=True)
#basic polarimetric reductions
logfile = 'specpol' + obsdate + '.log'
#wavelength map
infilelist = sorted(glob.glob('m*fits'))
linelistlib = ""
specpolwavmap(infilelist, linelistlib=linelistlib, logfile=logfile)
#background subtraction and extraction
infilelist = sorted(glob.glob('wm*fits'))
dyasec = 5. # star +/-5, bkg= +/-(25-35)arcsec: 2nd order is 9-20 arcsec away
yoffo = 0. # optional offset of target (bins) from brightest in O (bottom) image
specpolextract_sc(infilelist, yoffo, dyasec, logfile=logfile)
#raw stokes
infilelist = sorted(glob.glob('e*fits'))
specpolrawstokes(infilelist, logfile=logfile)
#final stokes
infilelist = sorted(glob.glob('*_h*.fits'))
specpolfinalstokes(infilelist, logfile=logfile)
#basic polarimetric reductions
# debug=True
debug=False
logfile='specpol'+obsdate+'.log'
#target and wavelength map
infile_list = sorted(glob.glob('m*fits'))
linelistlib=""
specpolwavmap(infile_list, linelistlib=linelistlib, inter=True, logfile=logfile)
#background subtraction and extraction
#infile_list = sorted(glob.glob('wm*099.fits')+glob.glob('wm*10[0-6].fits')) # optional subselection
infile_list = sorted(glob.glob('wm*fits'))
specpolextract(infile_list, logfile=logfile, debug=debug)
#raw stokes
#infile_list = sorted(glob.glob('e*06[8-9].fits')+glob.glob('e*07[0-1].fits')+glob.glob('e*07[4-7].fits')) # optional subselection
infile_list = sorted(glob.glob('e*fits'))
specpolrawstokes(infile_list, logfile=logfile)
#final stokes
#polcal = 'polcal0.txt' # null calibration
#infile_list = sorted(glob.glob('*_h[0,2]*.fits')) # optional subselection
polcal = 'polcal.txt'
infile_list = sorted(glob.glob('*_h*.fits'))
specpolfinalstokes(infile_list, polcal=polcal, logfile=logfile, debug=debug)