header, data_dict = ltools.from_HDF(filepath[0], datatypes)
try:
BRdat = data_dict[datatypes[0]]
except KeyError:
sys.exit("No BR532 data available!")
try:
PRdat = data_dict[datatypes[1]]
except KeyError:
sys.exit("No PR532 data available!")
altrange = np.arange(150,10000,10)
BRdat = ltools.alt_resample(BRdat,altrange)
PRdat = ltools.alt_resample(PRdat,altrange)
altrange = BRdat.columns.values
copolvals = np.hstack(BRdat.values).astype('float32')
depolvals = np.hstack(PRdat.values).astype('float32')
numbins = 100
depolmin = 0.0
depolmax = 0.5
copolmin = 0.0
copolmax = 5.0
copolhist=h2d.fullhist(copolvals,numbins,copolmin,copolmax,-9999.,-8888.)
#double check to make sure the mask files match up with the data files
for d,m in zip(datafiles, maskfiles):
[d_stat,d_date,d_type] = d.split('_')
[m_stat,m_date,m_type] = m.split('_')
print 'Checking mask/data match for %s'%(d_date)
if d_date == m_date and d_stat == m_stat:
print 'Check!'
continue
else:
sys.exit("Error: Mask files don't match data files")
#open, altitude resample, and concatenate data and mask files
for d,m in zip(datafiles, maskfiles):
d_temp, data_prod = LNC.lnc_reader(d)
d_realt = LNC.alt_resample(d_temp,altrange)
try:
d_event = pan.concat([d_event,d_realt])
except NameError:
d_event = d_realt
m_temp, data_prod = LNC.lnc_reader(m)
m_realt = LNC.alt_resample(m_temp,altrange)
try:
m_event = pan.concat([m_event,m_realt])
except NameError:
m_event = m_realt
BRfiles.append(f)
#make sure the files are in a common order of ascending date (assuming they're all
#from the same station
BRfiles.sort()
counter = 0
header = {}
#open, altitude resample, and concatenate data and mask files
dataprods = []
counter = 0
for br in BRfiles:
br_temp, tempprod = LNC.lnc_reader(br)
br_realt = LNC.alt_resample(br_temp,altrange)
if counter == 0:
br_event = br_realt
dataprods.append(BR_type)
else:
br_event = pan.concat([br_event,br_realt])
counter +=1
altrange = br_realt.columns.values
#sort by index to make certain data is in order then set date ranges to match
br_event = br_event.sort_index()
